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A review is given of attempts to bridge the gap between everyday par- 
ticle and nuclear physics — involving many quarks — and the basic 
underlying theory of QCD that can only be evaluated exactly for few 
quark systems. Even the latter requires the original theory of QCD to 
be discretised to give Lattice QCD — but this modification can still yield 
exact results for the original theory. These LQCD results can then be 
considered on a similar footing to experimental data — namely as cor- 
nerstones that must be fitted by phenomenological models. In this way, 
the hope is that "QCD inspired" models can become more and more 
"QCD based" models, by fixing — in the few-quark case where LQCD 
can be carried out — the form of these models in such a way that they 
can be extended to multi-quark systems. 



1.1 Introduction 

Even though for over 30 years QCD has been thought to be the theory 
of strong interactions, it has had a rather limited impact on most other 
branches of physics — except for few-quark hadron physics. Of course, 
the reason is well known — to write down the Lagrangian that describes 
exactly the quark and gluon interactions is easy, but actually performing 
calculations directly with this Lagrangian has turned out to be extremely 
difficult. The one exception to this last statement is at high energies, where 
— due to asymptotic freedom — the interactions become sufficiently weak 
for perturbation theory to be applicable and this has had much success 
pp. However, most of "everyday" physics is far from this limit. Further- 
more, some quantities such as masses depend on the interaction strength 
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2 Bridges from Lattice QCD to Nuclear Physics 

g as exp(— which immediately rules out a perturbation expansion in 
powers of g. 

This inability to treat the QCD Lagrangian directly has led to several 
different types of approximation being made. Essentially these fall into two 
broad categories: the numerical and the effective theory approaches. Unfor- 
tunately, the latter — especially the Effective Potential Theories (EPTs), 
which are the main subject of this chapter — often have little overlap with 
the numerical approaches. Those approaches that use the direct numerical 
way concentrate mainly on the description of single mesons or baryons i. e. 
qq or qqq states, whereas most of the EPTs tend to concentrate more on 
multiquark systems. However, even though these EPTs are often adver- 
tized as being "QCD motivated", in most cases they are simply based on 
many-body ideas and techniques that are well founded in nuclear physics 
— but are not necessarily justified for the description of multiquark states. 
The main purpose of this chapter is to see to what extent these "nuclear 
physics motivated" methods are justified. But first a few general words 
should be said about these two approaches. 



1.1.1 Numerical treatment of QCD 

The numerical treatment of the QCD Lagrangian — Lattice QCD — has 
been the main subject of this volume and so, to avoid too much repeti- 
tion, only points relevant to this chapter will be mentioned. In Lattice 
QCD the original exact Lagrangian is replaced by an approximate form 
that is discretized on a 4-dimensional lattice with links of length a. This 
discretization can be done in many ways, but in all cases the original QCD 
Lagrangian must be recovered as a — s- 0. The discretization also reduces the 
subject to being mainly numerical. However, as emphasized by Liischer j2|: 
" In general, numerical simulations have the reputation of being an approx- 
imate method that mainly serves to obtain qualitative information on the 
behaviour of complex systems. This is, however, not so in lattice QCD, 
where the simulations produce results that are exact (on a given lattice) up 
to statistical errors. The systematic uncertainties related to the non-zero 
lattice spacing and the finite lattice volume then still need to be investi- 
gated, but these effects are theoretically well understood and can usually 
be brought under control." Therefore, in order to recover results that are 
appropriate to the original continuum Lagrangian, two main limits need to 
be studied: 
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Introduction 3 

• Limit 1. Are the results stable as a 0? 

Of course, this limit must be approaclied with consideration of the 
number of spatial sites in the lattice N^, since the volume {V) of 
the system being studied, e.g. 2-, 3- or 4-quarks, should be much 
smaller than the physical lattice size = (aN)^. Furthermore, 
the Euclidean time — needed to extract observables — should be 
much greater than L. In practice, usually T ^ 2L suffices. There- 
fore, the two inequalities that need to be satisfied can be combined 
as <^L<^T. This can instantly lead to problems for a meson- 
meson system {q^q^) with mesons of size w 1 fm, since this would 
need V^^^ > 2 fm, if the gradual separation of the mesons is of 
interest. In this lattice spacing of « 0.1 fm would require 

iV > 20 and so T > 40. Such large lattices are used by some groups 
e.g. 32"^ x 60 in the study of scattering lengths in Ref. |3] — with 
even larger lattices 64'^ x 128 now becoming feasible |2. In fact, 
the progress in computer technology allows the lattice extents to 
be doubled in all directions roughly every 8 years [^J. Unfortu- 
nately, at present many of us have to be satisfied with sizes more 
like 16'^ X 24, which rules out the study of completely separated 
mesons. One way of partially overcoming this problem is to use so 
called Improved Actions. These incorporate modifications to the 
standard lattice Lagrangian in order to remove the lowest order 
dependences on a. This enables coarser lattices to be used, so that, 
in some cases, a ~ 0.5 fm suffices compared with more usual values 
of a « 0.1 fm — see Refs. [J. In this way, the inequality V^^"^ <C L 
is satisfied by increasing a and not N. Of course, there is often a 
price to pay. Since improved actions are more complicated, they 
take more computer time to implement compared with the stan- 
dard actions. Sometimes this is sufficient to remove the advantage 
of using a smaller lattice. Having said that, there are some im- 
proved actions that seem to be always advantageous. Probably the 
most common of these is the so-called clover action for improving 
the quark part of the QCD Lagrangian |5I ()| — see Subsec. 11.1.2.11 
and also Appendix B.5 of Chapter 4. 

Often the a ^ limit is checked in three steps: 



— i) First a benchmark calculation is performed with, say, a 
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4 Bridges from Lattice QCD to Nuclear Physics 

lattice 16'^ x 24 and a — 0.12 fm. This is the least time and 
storage consuming of these three steps. 

— ii) Then the finite size effect V-^^^ <^ L is checked with a 
larger lattice, say, 24'^ x 32 but with the same a = 0.12 fm. 
This is much more time and storage consuming by a factor 
of about 4, but will show whether or not the process being 
treated "fits" into the lattice. It is concluded that the finite 
size effect is no problem, if these latter results are directly the 
same, within error bars, as the benchmark results. 

— iii) Finally, the calculation is repeated with the above larger 
lattice but using, say, a' = 0.08 fm and the results compared 
with the above by now including appropriate factors of a' /a. 
These final scaled results should now agree, within error bars, 
with the results from the above two a = 0.12 fm lattices. 

A specific example of this scaling procedure appears later in Sub- 
sec. 11.4.3.31 This raises the question: Given a fixed number of 
"computer units" for a calculation, then what is the most efficient 
way of spending these units? As pointed out by Kronfeld |7j , it is 
much more efficient to run at several lattice spacings than to put 
all the resources onto the finest conceivable lattice. Kronfeld gives 
the example of a computer budget of 100 units and suggests using 
65, 25 and 10 units on a series of coarser lattices with spacings ag, 
oi = and a2 = The time needed to create statisti- 

cally independent lattice gauge fields grows as Tg cx a~("'+^\ where 
the 4 in the exponent arises because the number of variables to 
process grows as in a 3-1-1 dimensional world and z « 1 — 2 de- 
pending on the algorithm for updating the lattice. In this case the 
three sets of lattices would have comparable error bars that would 
only be shghtly larger — by a factor of about 1.25 (i.e. l/VO.65) 
— than the case of using all 100 units on the finest lattice with 
ao- However, using all three sets enables an estimate to be made 
of the discretization effect. As Kronfeld says "The slightly larger 
statistical error seems a small price to pay" . 

• Limit 2. The mass of the Hght quarks should be realistic. 

This seems to be a much more difficult limit to achieve. In practice, 
light quarks (u, d) with a mass ^ 100 MeV are often used instead 
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of the true values of less than 10 MeV. This is reflected in the 
computed ratio of the tt— and p— masses R^p = niTr/mp. The ex- 
perimental masses give Rj^p = 0.2, whereas Rt^p ~ 0.7, if rriu.d ~ n^s 
— the strange quark mass (a value used in many works). Unfortu- 
nately, extrapolating from results using different bare quark masses 
to get the observed i?,rp = 0.2 is not straightforward. At sufficiently 
low values of the bare quark masses the effective field theory of 
Chiral Perturbation Theory becomes applicable and shows that in 
these quark mass extrapolations logarithmic terms arise in addition 
to simple power-law behaviour — see Subsec. 11.1.2.11 

The above limits are discussed in more detail in Ref. |H]. There it is 
pointed out that QCD is a multiscale problem. Not only is there a charac- 
teristic scale of QCD (Aqcd ~ 200 — 250 MeV 0) but also a wide range 
of quark masses with light quark masses « 10 MeV up to toq w 5 GeV 
for the 6-quark, leading to the hierarchy 

mg < Aqcd < mq. (1.1) 

But two more scales are needed before QCD can be put on a lattice. Firstly, 
for light quarks the lattice size (L) must be larger than the size of a light 
quark i.e. <C L. Secondly, for heavy quarks the lattice spacing (a) 

must be finer than the size of such quarks i.e. a <C mg^- The hierachy in 
Eg. II. II then becomes 

< < Aqc_d < mg < (1.2) 

However, for numerical reasons it is not possible to satisfy all these con- 
ditions. So that, in practice, finite computer resources force the hierachy 



L 1 < < Aqcd < mg - a ^ (1.3) 

instead of the idealized one. 

In Ref. ^ the dependence on a and the appropriate value of TOtt for the 
work required to obtain a "new" configuration are combined into the single 
approximate expression 

- 0.157 (iV"n'"m"', 



config J \ o 
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6 Bridges from Lattice QCD to Nuclear Physics 

where 1 Gflop is 10^ computer operations per second. This shows the 
~ dependence that is the main numerical problem for Lattice QCD 
to overcome and is the reason for being interested in using coarser lattices 
with improved actions to be discussed in the next Subsection. Estimates 
similar to Eg. 11.41 are also made in Refs. 2 10 . 

For those readers who would like a detailed development of Lattice QCD 
the text books by Creutz Montvay and Miinster j^, and Rothe [T^ 
are recommended. Also there are many review articles and summer school 
lecture series — see Ref. for a partial listing of these. At the time of 
writing, some of the most recent reviews for a general audience are listed 
in Ref. 

1.1.2 Effective Field/ Potential Theories 

Effective theory formulations describing quark-gluon systems fall into dis- 
tinct categories. On the one extreme are the Effective Field Theories 
(EFTs) that have a rigorous basis, whereas at the other extreme we have 
the Effective Potential Theories (EFTs), which are essentially phenomeno- 
logical being based on models with potentials in differential equations. It is 
important to discuss these two types of theory separately, since they play 
very different roles in the present chapter and should not be confused with 
each other. 

1.1.2.1 Effective Field Theories (EFTs) 

Effective field theories play a crucial part in extracting continuum (physical) 
results from the purely numerical lattice techniques. A review of this topic 
has been given by Kronfeld ^7,. 

In these theories an energy scale (A) is introduced. This essentially 
separates the short distance effects {i.e. less than 1/A), which are lumped 
into the coupling constants of the theory, from long distance effects, which 
are described explicitly by the operators of the theory. Such theories are 
then only applicable to processes involving energies less than A. In QCD 
the energy scale characteristic of non-perturbative effects is ~ 1 GeV. There 
are several theories that fall into this category, examples of which are: 



*A popular level review in the February 2004 edition of Physics Today 14. has not been 
well received by everyone 1151 . 
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Introduction 7 

Symanzik effective field theory 

The most obvious difference between Lattice QCD and the real life 
situation of continuum QCD is the presence of the lattice with spac- 
ing a. Only the a limit has a physical meaning. A systematic 
way of studying this lattice artifact was developed by Symanzik \W\ . 
He showed how 0{a) effects could be removed in a systematic way 
from lattice results by assuming Aa is small and treating lattice 
artifacts as perturbations. He achieves this through creating an 
effective field theory by adding terms with increasing powers of a 
(and containing parameters Ci) to the basic lattice QCD Lagrangian 
of Wilson. The Ci parameters are then adjusted (tuned) to kill off 
the offending 0{a) effects. This has now been developed into an 
industry for generating improved actions that only contain O(a^) 
lattice spacing corrections — see Subsec. 4.1.7 and Appendix B of 
Chapter 4 for a more detailed description and for references. 
A simple example of this is the quark-gluon coupling T^{p,p'), 
where p, p' are the 4-momenta of the initial and final quark. The 
above strategy is to first replace the continuum coupling T'j^{p,p') 
by its lattice counterpart T^{p,p') and then to expand the latter 
in powers of a i.e. 

r^(p,p') = c7M-^ 

^jliPiP') = c{7^cos[i(p + p')^a] - ism[^{p + p')^a]} 
~^c{j^,~ ^a{p + p')^ +0{a^)}. 

To remove the 0{a) term, Sheikholeslami and Wohlert [5] sug- 
gested adding a lattice form of (J^h^F^" to the Wilson action so 
that r^(p,p') became 

(P, P') = c{7/. cos[i (p + p')f,a] - i sin[i {p + p')^,a] + 

]j^cswc^tLu cos[ifc^a] sin[ifc^.a]}. (1.5) 

Expanding Tf^{p,p') — >c{'y^- ^a[{p + p')^ + cswi'^f,^k''] + 0{a'^)}, 
where k = p' — p and ap,i, = i[jfj_,j^]/2. On the mass shell, if 
csw is now "tuned" to unity, then the two terms of 0{a) cancel to 
leave corrections of only O(a^). When this procedure is applied to 
Wilson's fermion action the outcome is usually referred to as the 
clover action. 
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8 Bridges from Lattice QCD to Nuclear Physics 

(2) Chiral Perturbation theory of Gasser and Leutwyler |17j 

Light quarks have a mass mi « 10 MeV i.e. mi <C A. This makes it 
numerically impractical to perform lattice QCD calculations with 
such masses since the algorithms for computing the quark propaga- 
tors become slower and slower — as seen from Eq. 11.41 Therefore, 
the procedure to reach this physical region is to first carry out the 
lattice calculations with a sequence of masses (wg) in a range of, 
say, 0.2ms < mg < ms, where ms 100 MeV - appropriate for 
the strange quark. Given this sequence of results (masses or matrix 
elements) , the task is then to use a reliable method for extrapolat- 
ing these results to quark masses appropriate for the light quarks. 
By far the most successful method for this extrapolation is based 
on Chiral Perturbation Theory (^PT), which can be viewed as an 
expansion in mq/A. The numerical data, with mq in the range 
0.2ms < mq < ms, can then be tested against the leading order 
(next-to-lcading order or next-to-next-to-leading order) prediction 
of xPT. If this is successful, then it gives confidence in extrapolating 
mq to the light quark masses |18J. 

(3) Heavy-quark efTective theory and Non-relativistic QCD 
For a brief review of Heavy-quark effective theory (HQET) and 
Non-relativistic QCD (NRQCD) see Ref. ^^l- In situations involv- 
ing heavy quarks — such as i?-physics, where some of the quarks 
have a mass of mg « 5 GeV and are non-relativistic — it is appro- 
priate to make expansions in terms of A/mg or ly, the relative ve- 
locity between the quark and the antiquark in the i3-meson. These 
two expansions are usually referred to, respectively, as HQET — 
for systems containing a single heavy-quark — and NRQCD for a 
heavy-quark heavy- antiquark system. One way of deriving these 
effective theories is to write down the heavy-quark theory as an 
expansion in terms of A/mg or ly in the continuum and then re- 
place the derivatives that arise by their lattice counterparts to give 
Lattice HQET and Lattice NRQCD. These ideas are still being de- 
veloped. For example in Ref. [SOj some of the irrelevant degrees 
of freedom in NRQCD are integrated out to yield a theory called 
potential NRQCD (or pNRQCD), which is much simpler to treat. 

EFTs are not only used with few-quark systems but also have a long his- 
tory in few- and many-nucleon systems — see the works of van Kolck et 
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al. |21l 22\ . This approach was first advocated by Weinberg [23], who il- 
lustrated how the nucleon-nucleon and many-nucleon potentials could be 
qualitatively understood. For example, these arguments show that, if the 
strength of the NN-potential is 10 MeV, then those of the NNN- and 
NNNN-potentials are ^ 0.5 and ^ 0.02 MeV respectively — numbers that 
are in accord with detailed few-nucleon phenomenology based on realis- 
tic potentials such as that of Argonne [21] • However, it is not clear that 
this approach — in spite of its impressive qualitative results — could 
ever compete quantitatively with standard meson exchange models for 
describing the NN-potential, where baryon resonances such as the A(1236) 
and iV*(1535) are included explicitly. As pointed out by van Kolck himself, 
even the inclusion of the p- and cj-mesons give rise to interactions that are 
"at present an insurmountable obstacle for a systematic approach" . 

In multi-particle systems EFTs are often converted into Mean Field 
Theories (MFTs), in which the emphasis is on single particle properties 
with all other particles being treated "on the average" . Unfortunately, in 
some applications — such as the equation of state of high density nuclear 
matter, as encountered in relativistic heavy ion collisions or in neutron stars 
— there are serious questions concerning their validity. Even advocates of 
the MFT approach (see, for example, Glendenning on pages 127 and 287 
in Ref. j25j ) express reservations by writing "In many ways it is not as 
good a theory as the Schrodinger-based theory of nuclear physics" and 
"The status of an exotic solution of an effective theory is more tenuous 
than from a fundamental theory" . Others are even more critical — see 
for example Ref. |21j: "The Relativistic Mean Field approximation is very 
elegant and pedagogically useful, but is not valid in the context of what is 
known about nuclear forces • • • . It requires /i (r) <^ 1, where (r) is the 
average interparticlc distance and l//i a meson range. However, for pions 
Ptt {r) ~ 0.8 — 1.4 but for vector mesons /i^ (r) ~ 4.7 — 5.8." 



1.1.2.2 Effective Potential Theories (EPTs) 

The reason for briefly describing the above Effective Field Theories (EFTs) 
is to emphasize their difference from Effective Potential Theories (EPTs), 
which are the main interest in most of this chapter. The above EFTs are 
an integral part of the development of Lattice QCD and play a crucial role 
in extracting precise continuum results for few (2 or 3) quark systems i.e. 
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10 Bridges from Lattice QCD to Nuclear Physics 

Lattice QCD+EFTs — > Continuum results for few quark systems. 

On the other hand, the EPTs attempt to understand (interpret) these con- 
tinuum results for few quark systems in such a way that the theory can be 
extended to the multiquark systems of interest to nuclear physics i.e. 

Continuum results for few quark systems + EPTs — > 
Descriptions of multiquark systems (Hadron Physics). 

This step is here referred to as "Bridges from Lattice QCD to Nuclear 
Physics" — the title of this chapter. In most cases this step is mainly 
phenomenological. 

The EFTs mainly concentrate on the properties of a single particle, so 
that, for example, the energy of a multiparticle system is expressed as the 
sum of the effective masses of the separate particles with the effect of all 
other particles being treated in an average manner. In this way symmetries 
of the fundamental Lagrangian can be preserved. However, as said above, 
the extension of these ideas to more complicated Lagrangians or many-body 
systems presents problems. On the other hand, in the less ambitious and 
more phenomenological approach of Effective Potential Theories (EPTs), 
the emphasis is first on the two-body system. In this case, two-body poten- 
tials are the main ingredient. 

For multi-nucleon systems the NN-potential can be mainly phenomeno- 
logical or based on EFTs as with the Argonne and Bonn potentials respec- 
tively | 24l 27) . Those based on EFTs can be generated with varying degrees 
of ability to describe the two-body system. We have already mentioned the 
works of Weinberg [23] and van Kolck , which follow the procedure — 
referred to by van Kolck as Weinberg's "theorem" : 

(1) Identify the relevant degrees of freedom and symmetries involved. 

(2) Construct the most general Lagrangian consistent with item (1). 

(3) Do standard quantum field theory with this Lagrangian. 

The outcome is qualitatively correct being within « 10% of the two-nucleon 
data — but this is an accuracy that is often insufficient for understanding 
nuclear phenomena. At the other extreme we have the one-boson-exchange 
(OBE) potentials in which various meson-baryon couplings are tuned to 
ensure a good fit to the NN experimental data. These OBE potentials, 
even though they are based on EFT- like Lagrangians, often incorporate 
couplings such as A^A(1236)p that can not be treated systematically by 
Weinberg's "theorem" . 
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In contrast, to implement EPTs for multiquark systems three ingredi- 
ents are necessary — a wave equation (differential or integrodifferential) , 
an interquark potential and effective quark masses: 

A wave equation. 

Since EPTs are not derived from more basic principles, the forms of the 
wave equations are not predetermined and can vary considerably. Even for 
two-quark systems there is a choice. 

• A non-relativistic Schrddinger equation. This is suitable for heavy- 
quark mesons (QQ) with mg 3> 1 GeV and was the form used in 
Refs. |28l 29j for extracting the QQ-potentials. The best cases for 
this are the Bottomonium mesons such as the T(66, 9.5GeV) since 
mfe ~ 5 GeV. 

• The Dirac equation. Once the quarks are light, i.e. mq <C 1 GeV, 
relativistic effects become important and we, therefore, enter the 
realm of large/small wave function components, pair-creation and 
quantum field theory. This means that the use of the Dirac equation 
is less "clean" phenomenologically than the Schrodinger equation, 
since it deals explicitly with large/small wave function components 
but not with the related effect of pair-creation. In spite of this, 
for Qq systems such as the B-meson, where one quark is heavy, 
the Dirac equation is essentially a one-body equation and so the 
full complications of the two-body relativistic problem are avoided. 
There are many references where this one-body Dirac equation has 
been applied, e.g. [HOI :i II 321 3^. 

• The Bethe-Salpeter equation. When the two quarks are both light 
the correct relativistic scattering equation is the full Bethe-Salpeter 
equation. Unfortunately, direct use of this equation — for physi- 
cally interesting cases — presents severe problems (see Ref. for 
a recent discussion). Therefore, it is usually reduced from a four- 
to a three-dimensional scattering equation by inserting appropri- 
ate (5-functions of the energy. This can be carried out in several 
ways and leads to a number of different two-body equations that 
are covariant and satisfy relativistic unitarity. Examples are the 
equations of Blankenbecler-Sugar, Gross, Kadyshevsky, Thomp- 
son, Erkelenz-Holinde, ...see Chapter 6 in Ref. [33]. Depending 
on the problem, these alternatives have their separate advantages. 
For example, the Gross form — unlike those of Blankenbecler- 
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Sugar and Thompson — treats the two quarks asymmetrically and 
has the feature that it reduces to the Dirac equation when one of 
the quarks becomes infinitely heavy. This suggests that this form 
is perhaps more appropriate for describing the i3-meson. On the 
other hand, when the quarks are equal in mass, as in the J /tp{cc) 
system, the Blankenbecler-Sugar or Thompson equations are prob- 
ably preferable. 

• Multi-quark wave equations. For multi-quark systems the choice 
of wave equation is very limited with the Resonating Group ap- 
proach being the most usual — see the review by Oka and Yazaki 
in Chapter 6 of Ref. and also more recently Refs. \X7l 38j . This 
reduces the interaction between quark clusters {q^ + q^ for the 
NN-potential and qq -\- qq for meson-meson interactions) to a non- 
local Schrodinger-like equation involving the relative distance be- 
tween the clusters. This is achieved by integrating out the explicit 
quark degrees of freedom, which usually requires the introduction 
of gaussian radial factors in order to carry out the multiple integrals 
involved. At first sight such factors may seem to be unrealistic with 
exponential or Yukawa forms being more physical. However, it will 
be seen in Subsec. 11.9.3.21 that the one-body Dirac equation with 
a linear confining potential can lead to gaussian forms asymptot- 
ically. Unfortunately, the effective mass needed for the quarks is 
« 300-400 MeV and so makes this non-relativistic approach some- 
what questionable. On the other hand, the role of relativity in 
many-body systems is still an open question — a recent summary 
being given by Coester EH]. 

An interquark potential (Vqq) — the second EPT ingredient. 

For two-quark systems the interquark potential is often taken to have the 
form suggested by the static limit of infinitely heavy quarks, namely 



where the first term is that expected at short distances due to one-gluon 
exchange and the second term that expected from quark confinement — 
c simply being an additive constant. However, it should be emphasized that 
the form in Eg. II. (il ls, strictly speaking, only appropriate for the interaction 
between static quarks. At the present time, there seems to be no really con- 



(1.6) 
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vincing evidence for a significant one-gluon exchange interaction in systems 
where only hght constituent quarks (up to ~40Q MeV) are involved. This 
becomes evident in the N and A spectra. There a one-gluon interaction is 
unable to describe the empirical ordering of the positive and negative par- 
ity states. However, this ordering can be accomplished by Goldstone-boson 
exchange mechanisms 00]; even though that model also has bad features — 
some states in the A and A spectra are poorly reproduced. The reason why 
the one-gluon exchange becomes ineffective with light quarks is because the 
use of minimal relativity — in, say, the Blankenbecler-Sugar equation for 
describing the Qq interaction — introduces relativistic square root factors. 
Indeed, if both quarks are light, then the effective one-gluon interaction 
is essentially flat and very weak for short distances. Since this damping 
of one-gluon exchange also enters in meson spectra, another mechanism is 
needed for the necessary short range attraction. A possible candidate for 
this is an effective interaction generated by instantons, which is usually 
expressed in terms of an attractive J-function in r 41 . 

More details on VQQ{r) can be found in Chapter 3, which is devoted to 
this topic. The above is a two-body potential. However, for multi-quark 
systems there are strong indications that multi-quark potentials and/or 
potentials involving excited gluon states could also play a major role — see 
Sees. O and O 

For light quarks and for nucleons in dense matter, relativistic effects 
enter and in some cases these can be expressed as corrections to the non- 
relativistic potential involved. An example of this is the so-called Relativis- 
tic Boost Correction shown to be important in high density nuclear matter 

m 

EfTective quark masses — the third ingredient for an EPT. 

The masses of the quarks involved in this "Wave Equation -I- Potential" 
approach are not those of bare quarks i.e. not « 10 MeV for the u and d 
quarks. They are essentially free parameters, which are often taken to be 
~ ^^nucicon/3 ~ 300 McV. Howevcr, it has been suggested 03j that a more 
natural choice would be w A/a/S ~ 400 MeV. 

So we see that the idea of EPTs covers an enormous number of theories, 
models and approaches that are frequently used in nuclear physics. The 
subject of the next section is to see how these EPT ideas can possibly be 
utilized in the understanding of QCD. The magnitude of this step should 
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not be underestimated, since as the authors of Ref. say: 

"On the one hand, the concensus of the majority of the nuclear physics 
community holds that in nuclei 

• nucleons are non-relativistic 

• they interact via essentially two-body forces, with smaller contri- 
butions from many-body forces 

• the two-nucleon interaction generally possesses a high degree of 
isospin symmetry 

• external probes usually interact with mainly one nucleon at a time. 

By contrast, in QCD 

• the u, d and even s quarks are relativistic 

• the interaction is manifestly multi-body, involving exchange of mul- 
tiple gluons 

• there is no obvious isospin symmetry 

• external probes can, and often do, interact with many quarks at 
once. 

It should not be surprising, then, that some new ideas are required to merge 
these two extraordinarily different bodies of theory." Others might simply 
say that this is a case of "Mission impossible" . 

1.2 What Is Meant by "A Bridge"? 

The main theme of this article is the study of bridges between lattice QCD 
and nuclear physics. Unfortunately, what constitutes a possible bridge is 
rather subjective, since the basic idea is to compare some quantity that can 
be measured by lattice QCD with a "corresponding quantity" that arises 
in more conventional physics as the outcome of some EPT or is directly 
connected with experiment. Of course, the question instantly arises as to 
whether these two quantities are indeed comparable i.e. to what extent are 
we confident that we indeed have "corresponding quantities" . In this chap- 
ter the main quantities to be related will be energies or radial correlations. 
This is probably best illustrated by the following simple example. 
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1.2.1 A simple example of a bridge 

Since it is the main goal in this chapter, it may prove useful to the reader 
to first see a simple example of what is meant by "Bridges from LQCD to 
Nuclear Physics" . 

The results of any lattice QCD calculation are quantities expressed as 
dimensionless numbers. Therefore, to be able to make a connection with 
"real life" , one — or more — of these numbers must be compared with its 
continuum counterpart that can actually be measured experimentally. This 
then sets the physical scale for lattice QCD. 

1.2.1.1 Setting the scale from the string tension 

For many years a quantity frequently used for this comparison was the 
string tension (&s), which — as its name implies - is simply the energy/unit 
length of the flux-tube (i.e. string) connecting two quarks and appears 
in Eq. 11.61 Experimentally, estimates can be made of this string tension 
from the spectra of mesons and baryons with increasing orbital angular 
momentum {L) — a series of energies (E) that depend crucially on the 
string increasing in length. This can be carried out with varying degrees 
of sophistication. By simply plotting L versus E"^ this so-called Regge 
trajectory is found to be linear for both mesons and baryons — the slope 
(a) in each case being about 0.9 GeV~^ — see, for example, Figs. 7.33 and 
7.34 in Ref. 01]. As shown in Ref. @ni using a simple classical model this 
slope is directly related to the string tension by the expression a « l/{8bs)- 
This results in a value of a/EI « 380 MeV — a number that is somewhat 
smaller than the accepted value of ~ 440 MeV, which is more in line with 
estimates from string models that give a ~ l/(27r&s) |46j . 

A less direct, but more precise, way to extract the string tension is 
to first find an effective quark-antiquark potential (Vqq) that describes 
— by way of a non-relativistic Schrodinger equation — the above meson 
energy spectra. Naturally, for this non-relativistic approach to be realistic 
the mesons must be constructed from quarks that are much heavier than 
the proton. This, therefore, restricts the analysis to the Bottomonium bb 
mesons, where the b quark has a mass of about 5 GeV and possibly the be 
and cc mesons, where the c quark has a mass of about 1.5 GeV. In fact, these 
spectra can be described by a single effective potential of the form given in 
Eq. 11.61 so that a value for the string tension of y/b^ « 440 MeV results — 
equivalent to 6^ « 1 GeV/fm or « 5 fm~^ in other units. The potentials 
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most frequently quoted are those of Richardson [2H| and Cornell ■ This 
value of bs is now the experimental number to which the lattice estimate of 
VQQ{r) must be matched. Usually the latter is extracted by measuring a 
rectangular Wilson loop W{1, t) of area a^lt for two infinitely heavy quarks 
a distance r = al apart on a lattice and propagating a Euclidean time at — 
a being the lattice spacing to be determined. Wilson loops will be discussed 
in more detail in Sec. 11.4.1.21 A key observation, first made by Wilson 02] 
in 1974, was that 

W{l,t) e^p[~tVQQ{l)] ast->oo. (1.7) 

Therefore, for sufficiently large I, Vqq{1) — > b'J, where 6^ is the dimension- 
less counterpart to the experimental string tension bg defined in Eq. 11.61 
The two are then matched by way of the "bridging equation" b'^ = a?bs to 
give a. A typical number for b'^ is « 0.05 giving a w -^/O.OS/S « 0.1 fm. 

1.2.1.2 Sommer's prescription for setting the scale 

In order to set the scale, the above has compared the lattice result with 
experiment for the most simple of quantities — the string tension. However, 
the experimental data is mainly probing distances of r « 0.2 fm to r « 1 fm 
and not r —> oo, since the rms-radii of the bb mesons cover the range from 
about 0.2 to 0.7 fm and the cc mesons the range from about 0.4 to 1 fm. 
Therefore, the experimental data encoded in the potential Vqq{1) is not 
optimal for studying the string tension. Furthermore, lattice calculations 
of the Wilson loop W{l,t) require t ^ I, so that those W{l,t) dominated 
by the string tension need to be evaluated for large values t. Unfortunately, 
as t increases the Signal/Noise ratio on W{1, t) also increases — eventually 
making measurements for large r = al meaningless. Therefore, on both 
the experimental and lattice sides there are problems for making a reliable 
estimate of a from the string tension. 

In an attempt to overcome this problem, Sommer j48| proposed com- 
paring the potential Vqq{1) as extracted from the lattice with that from 
experiment. However, to use the words of Sommer, "We must remember 
that the relationship between the static QCD potential and the effective 
potential used in phenomenology is not well understood." In spite of this, 
he suggests that the comparison be made at some value of r = in the 
optimal experimental range of r « 0.5 fm. Also at these values of r, Vqq{1) 
can be more reliably extracted on a lattice. In practice, it is the force. 
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defined essentially as 



F{1) 



(1.8) 



a 



that is compared through the expression 



r^i^(r')|r.=fi(c) = c. 



(1.9) 



Sommer chose the dimensionless parameter c = 1.65, since — for the exper- 
imental potentials — this corresponds to a distance i?(1.65) = i?o ~ 0-5 fm. 
Using the lattice forms of F{r) that fit the lattice potential, Ea. ll.9l can be 
solved for r in lattice units a. Comparing this r with Rq then gives a — see 
Ref. 0H1 for more technical details. However, it should be added that the 
value of c is somewhat uncertain with some authors |49| preferring c = 2.44, 
which corresponds to Rq « 0.66 fm. This second choice of c gives values 
of a that are a few percent larger than before and also in better agreement 
with the string tension estimate. 

The reason for this rather lengthy description for extracting the scale 
a is to show that "A Bridge from Lattice QCD to Nuclear Physics" has 
existed for many years. It should be added that another way of extracting 
a, when dealing with light quarks, is to use directly the mass of the p-meson 
as a corner stone and simply compare the experimental mass nip with the 
outcome of the lattice QCD calculation — the dimensionless combination 
airip 50 . However, is known to be a sensitive indicator of scaling 
violations i.e. how the lattice results depend on a as a ^ 0. It is, therefore, 
sometimes reserved for this purpose with the above method of Sommer 
being used to extract actual values of a "ST. The reason for choosing mp 
and not is because the 7r-meson, being so light, is more difficult to treat 
on a lattice. 

1.2.2 Are there bridges other than Vqq ? 

The above simple example showed how the QQ potential Vqq could be 
related to its lattice QCD counterpart and serve as a means for extracting 
the lattice spacing a. The question then arises concerning the possibility 
of there being other quantities that could be compared. However, it must 
be noted that Vqq and the related force F{r) are somewhat special and 
that it is still true what Sommer wrote in 1993: "As to today's knowledge, 
the force F(r) between two static quarks is the quantity which can be 
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Table 1.1 Possible bridges between lattice QCD and nuclear physics 



System 


Quantity matched 


Model 


Refs. 


(QQ) 


String Tension 


Regge Trajectory 
Schrodinger Equation 


m 

EHI 


[{QQ){QQ)] 


Energies 
Flux tube structure 


Matrix diagonalistion 

Discretized String 
Dual Potential Model 


|h;>|h8| 

m 


(Qq) 


Energies 
Density distributions 


Dirac Equation 
Dirac Equation 


m-m 


[{Qq)iQq)] 


Energies 


Variational 


m 


(QQq) 


Energies 
Density distributions 




m 
m 



calculated most precisely" . There are several reasons for this and they 
should be kept in mind in the following discussion. Firstly, both the lattice 
and experimental determination of Vqq can be done with good statistical 
precision and, secondly, the two are what we think they are. In contrast, the 
string tension can only be extracted at values of r that are not necessarily 
sufficiently asymptotic and where there could be corrections from model 
dependent sub-leading terms. So for the purposes of setting a scale the use 
of Vqq is still the best. 

Possibilities for bridges, in addition to the use of Vqq, are listed in 
Table 11.11 This is essentially a "Table of Contents" for the rest of this 
chapter. The Lattice QCD Nuclear Physics relationship changes as we 
go through this list. The first two rows for the QQ static quark system 
have been discussed above. Here the roles of the string tension and Vqq 
are to set a scale for QCD — a necessary step in order for lattice QCD to be 
compared with experiment i.e. the flow of information is Nuclear Physics 
Lattice QCD. However, once the results of lattice QCD can be expressed 
reliably with physical dimensions then the information flow is completely 
Lattice QCD Nuclear Physics. We can now consider the results of lattice 
QCD on the same footing as experimental data — assuming that the a — > 
limit is under control and that the quarks are sufficiently light as discussed 
in Subsec. H.l.ll 

It is lattice data that models must attempt to fit. Many of these mod- 
els resort to the use of interquark potentials — often the above Vqq — 
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in various forms of wave equation. Now the lattice QCD data will pos- 
sibly be able to justify — or rule out — such models. At present these 
models are often simply mimicking techniques that have proven successful 
in Nuclear Physics. Hence my earlier statement that they are "Nuclear 
Physics-inspired" and not "QCD-inspired" as is often claimed. 

Above I said that the results of lattice QCD can be considered on a 
similar footing as experimental data. However, when setting up models, 
in some ways lattice QCD data can sometimes be superior to experimen- 
tal data, since it can be generated in "unphysical worlds".^ Such worlds 
can have the following unphysical features that should, in some cases, also 
be inserted into the corresponding models to test the generality of these 
models: 

(1) The real world of three coloured quarks {i.e. SU(3)) can be replace 
by one with two coloured quarks {i.e. SU(2)). Such a world is eas- 
ier to deal with in lattice QCD — for example, there is essentially 
no distinction between quarks and antiquarks. Also the system 
corresponding to a baryon now consists of only two quarks. How- 
ever, this is not simply an academic exercise, since, in practice, it is 
found that the ratio of many observables are similar in both SU(2) 
and SU(3) — but with a computer effort that is about an order of 



magnitude smaller. An example of this is the ratio R — togb/a/STj 

where ttigb is the glueball mass (see Chapter 2) and bg the string 
tension. It is found for the glueball with the lowest mass (0+"'") that 
R w 3.5 for both SU(2) and SU(3) — see Sec. 2.2.1 in Chapter 2 
and Ref. [221 • In Ref- (HHI this is extended to the general case of 
SU(A^C') for several glueball states and for the 0++ case results in 
i?(0++) = 3.341(76) -t- 1.75/iV§. 
(2) The real world with 3 space coordinates and 1 time coordinate 
(3-1-1) can be replaced by one with 1 or 2 space coordinates and 
1 time coordinate (1-1-1 and 2-|-l). On the lattice the latter are eas- 
ier to study so that results with such high accuracy can be achieved 
that there is little ambiguity in any final conclusions. Also the 
(2-1-1) world has interesting features in its own right and enables 
comparisons to be made between SU(2), SU(3), SU(4),. . . , SU(iVc) 




tin Sec. 3.1 of Chapter 3 these are referred to as "virtual worlds" 
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(3) In the real world space is isotropic, but this need not be so on a lat- 
tice, since the four axes can be treated differently by having unequal 
lattice spacing — in principle we could have ^ ciy ^ at- 
However, the most common choice is Ux = a-y — ciz ^ at. This is 
appropriate for finite temperature systems |65| . where the temper- 
ature is defined to be inversely proportional to the lattice size in 
the t-direction. For high temperatures this would mean a lattice 
that contained fewer steps in the t-direction and so lead to difhcul- 
ties in extracting accurate correlation functions. However, if at is 
made smaller than the three spatial a's, then a given temperature 
is defined by more steps in the ^-direction and better correlation es- 
timates — see Ref. 22j- For studying high-momentum form factors 
such as i? ^ K*^, B irlv or B ^ plv it has been suggested that 
the 2-1-2 anisotropic lattice ax ^ ay ^ a^ ~ at is more suitable [HHI ■ 
Explicit anisotropic forms of the clover action — see Eq. 11.51 — 
and the pure gauge action in Eg. 11.171 can be found from Ref. |67| . 
More recent studies can be found in Refs. |68| . 

(4) In the real world the vacuum (sea) contains gg-pairs, where the q 
are sea-quarks that can have any flavour u, d, s, c, b or t. These 
pairs are being continuously created and annihilated. Lattice QCD 
calculations that take this into account are said to involve dynam- 
ical quarks. In practice, only u, d and, possibly, s sea-quarks 
are included and these two cases are usually referred to as having 
Nf = 2 or 3. However, frequently this effect is neglected to give the 
so-called quenched approximation, which is numerically an order 
of magnitude less demanding on computer resources. In view of 
this last point, there has been much work attempting to show how 
realistic quenched results can be compared with their dynamical 
quark counterparts. The conclusion seems to be that, although 
no formal connection has been established between full QCD and 
the quenched approximation, the similarity of the results ( ~ 10% 
differences) has led to the belief that the effects of quenching are 
generally small, so that quenched QCD provides a reasonable ap- 
proximation to the full theory This has been taken one step 
further in Ref. IJO] , where it is suggested how quenched results can 
be corrected in a systematic way to retrieve the corresponding full 
QCD prediction. 
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However, there are cases where the quenched versus unquenched 
comparison can lead to quaUtative differences. In Ref. [JJ the 
asymptotic potential between two quark clusters is calculated al- 
gebraically in a quenched approximation effective field theory and 
found to have a pure exponential decay and not the usual Yukawa 
form. But it is not clear that this qualitative difference for large in- 
tercluster distances leads to any overall quantitative differences. It 
must be remembered that, at these distances, calculations involving 
only quarks are expected to be incomplete with the introduction of, 
in particular, explicit pion fields being necessary. The same authors 
also study a partially— quenched effective field theory, where the 
masses of the sea-quarks are much larger than the valence quarks 
connected to the external sources [T^j- This they suggest would help 
in the understanding of the NN-potential, when extrapolations are 
made of NN-lattice calculations to realistic quark masses. 
(5) The real world of fixed quark masses can be replaced by one where 
the quark masses take on other values. In many cases, this is of ne- 
cessity, since for light hadrons the use of realistic light quark masses 
of « 10 MeV is computationally too heavy — see subsection II . 1 . II 
However, the variation of quark masses has interesting features in 
its own right. In particular, by carrying out lattice calculations 
with a series of light quark masses, we can extract the differential 
combination J = "lif^^^ry, where My and Mp are the corre- 
sponding vector and pseudoscalar meson masses and rriK* is the 
mass of the K*. This quantity J, which can be shown to be inde- 
pendent of a and the so-called hopping parameter that is related 
to the bare quark masses, serves as a check on the consistency of 
lattice QCD — see Refs. [73] and [TJ. Such analyses can be per- 
formed by varying the valence- and sea-quark masses separately. 
In Subsection. 11.9.41 a similar argument is used for estimating the 
matter sum rule. 

In Ref. |75| the authors have emphasized the importance of the 
quark-mass dependence of the nuclcon-nuclcon interaction by say- 
ing: "While the mg-dependence of the nuclear force is unrelated 
to present day observables, it is a fundamental aspect of nuclear 
physics, and in some sense serves as a benchmark for the devel- 
opment of a perturbative theory of nuclear forces. Having this 
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behaviour under control will be essential to any bridge between 
lattice QCD simulations and nuclear physics in the near future." 
In Ref. |76j the authors are more interested in the reverse situa- 
tion, namely, how to extrapolate nuclear forces calculated in the 
chiral limit to larger pion masses pertinent for the extraction of 
NN-observables from lattice calculations. 

(6) In the real world, hadron-hadron scattering is thought to be de- 
scribed directly in terms of quark-gluon physics at small inter- 
hadron distances, but at larger distances a description in terms of 
meson-exchange is expected to be more appropriate. Both of these 
limits must be included in a single model, if a direct comparison 
with experiment is to be made. However, lattice QCD can concen- 
trate on just the small distance physics and generate "data" that is 
exact not only for that limit but, in principle, for larger interhadron 
distances — until the numerical signal becomes unmeasurable. In 
this way, models can be constructed and compared directly with 
the lattice data — ignoring the effect of explicit meson exchange 
that enters in the real world. However, it should be added that, in 
some cases, lattice calculations based purely on quarks and gluons 
seem to be able to generate effects that resemble meson exchange. 
An example of this is Ref. [77| discussed in Subsec. 11.10.11 

(7) In the real world, space is essentially infinite, whereas quantities cal- 
culated on a lattice are restricted to volumes with L « 1 — 2 fm 
often being comparable to the size of the object under study. This 
leads to results that could depend on L. Usually this is considered a 
negative feature and so lattices must be chosen sufficiently large to 
avoid this problem. However, it was shown by Liischer in Refs. |7H] 
that this volume dependence can be utilized to extract the inter- 
action between hadrons. Recent examples of this idea consider the 
two-nucleon interaction [791 8()| and tttt scattering [HJ. This ap- 
proach is discussed in more detail in Chapter 4 Subsec. 4.1.5. 

(8) In the real world, the systems encountered contain a number of 
valence-quarks and antiquarks. However, on a lattice, systems con- 
sisting of only gluons can be studied — with quarks only entering 
as sea-quarks in the case of dynamical quarks mentioned in item 4 
above. This pure-glue world enables a cleaner study to be made of 
glueballs — see Sec. 2.2 of Chapter 2. 
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These possibilities of a different number of colours, spatial dimensions and 
quark masses greatly expand the scope of lattice QCD and give model 
builders much more data on which to test their models. 

1.3 The Energies of Four Static Quarks (QQQQ) 

1.3.1 Quark descriptions of hadron hadron interactions 

Much of particle and nuclear physics studies the interaction between hadrons 
With increasing complexity in terms of the number of quarks thought to 
be involved, this ranges from meson-meson scattering up to heavy ion col- 
lisions. Clearly, any attempt to describe such processes at the quark level 
must begin with an understanding of meson-meson scattering. Unfortu- 
nately, even with this system there are several major complications pre- 
venting a direct comparison between theory and experiment. Firstly, the 
only mesons for which there is suitable experimental data are the pseu- 
doscalars — the tt, K and B, since beams of these can now be generated 
at the various tt-, K- and B-factories [(tt) PSI (Villigen) and TRIUMF 
(Vancouver); {K) DAPHNE (Frascati) and KEK (Tsukuba); {B) BaBaR 
at SLAC, Belle at KEK, Hera-B at DESY and CLEO III at Cornell]. Of 
course, having beams of mesons does not lead directly to obtaining data on 
meson-meson scattering. This can only be done indirectly, as a final state 
interaction, with the net result that essentially only tttt scattering data (and 
considerably less ttK data) are at present available — and even those are 
very limited. This means that most theoretical attempts to imderstand 
meson-meson scattering concentrate on the tttt system. However, quark 
descriptions of this particular system are then complicated by the fact that 
the pion, being a Goldstone boson, does not have a quark structure as 
simple as a single qq configuration. Models of the pion (see, for exam- 
ple, Weise in Chapter 2 of Ref. j^Q) suggest that it has large multiquark 
components i.e. — X]n'^"(99)"- Furthermore, the total interaction be- 
tween the two pions can not be only due to interquark interactions between 
the constituent quarks, since for large interpion distances it is expected 
that meson exchange — another multiquark mechanism — also plays a 
role i. e. TTTr-scattering involves much more than a discussion of the [qq) (qq) 
system. Having said that, it should be added that these complications have 
not detered the construction of models for TTTr-scattering that are essentially 
nothing more than a {qq)n interacting with a {qq)n through an interquark 
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potential of the form in Eq. 11.61 The references are too numerous to Hst 
here and, furthermore, they often involve physicists who are my friends. In 
my opinion, these models are, as yet, not justified. They are simply hop- 
ing that the success in treating multi-hadron systems in terms of two-body 
potentials will repeat itself. 



1.3.2 The role of lattice QCD 

To make a bridge between quark and hadron descriptions of, say, meson- 
meson scattering needs reliable experimental data. But, as said above, this 
is not available — and this is where lattice QCD enters. The latter is 
based on QCD, which is thought to be the exact theory of quark-gluon 
interactions, and its implementation on a lattice leads (in principle) to ex- 
act results — upto the lattice spacing, lattice finite size and quark mass 
reservations mentioned in Subsec. [1.1.11 Therefore, if we want to study, for 
example, (qq) (qq) systems we simply calculate these on a lattice and we get 
exact results that can now be considered as "data" . Model builders then 
try to understand these data in terms of {qq){qq) states. Such a procedure 
guarantees one of the necessary requirements of bridge building — the need 
to compare like-with-like. In this way, the lattice data generated for the 
(99) (99) system can possibly be modelled with purely {qq){qq) configura- 
tions. In other words the conventional approach for model building 

1 . . 2 

Experimental data — > Hadron description — > Quark description 

is replaced by the alternative 

Lattice data — ^ Quark description Hadron description 

i.e. by concentrating on step 3 we avoid: a) At step 1 the shortage of ex- 
perimental data; b) At step 2 the need to guess the hadron quark structure 
and how a model based on this structure matches on to models more appro- 
priate at larger intcrhadron separation where meson exchange dominates. 
However, there are several problems when attempting to implement this 
second alternative: 

• Step 4 is similar to step 2 each with their uncertainty in the phys- 
ical hadron structure. However, now it is less serious, since step 3 
enables a cleaner description to be made at the quark level. This is 
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in contrast to the conventional approach, where the quark descrip- 
tion is "shielded" from the experimental data by needing to go via 
the hadron description, which could well also involve explicit meson 
exchange. 

• Step 3 is only feasible technically for a very restricted number of 
quark systems — four quarks being essentially the limit at present. 
However, there are a few six quark lattice calculations — mainly 
studies of the NN-interaction for static quarks (see Subsec. 4.5.2 
in Chapter 4) and, in addition, attempts to determine whether or 
not the H dibaryon is bound. In Ref. 1 821 the indications are that 
such a (uuddss) bound state is ruled out. Also recently there have 
been multiquark lattice calculations |83| with uudds quenched con- 
figurations in an attempt to describe the 0+(154O) seen in several 
experiments and thought to be the first observed pentaquark sys- 
tem P] . 

• The best lattice calculations are carried out with dynamical quarks 
— the so-called unquenched formulation. There the possibility 
arises for the creation (and annihilation) of quark-antiquark pairs. 
This is in contrast to the quenched approximation where such pairs 
do not enter. This means that in the unquenched approximation 
the configurations included do not have a fixed number of quarks 
and antiquarks. In model building this effect is often ignored and 
only a fixed type of configuration is used, e.g. {qq){qq) for meson- 
meson interactions. Fortunately, it is found that in most cases of 
interest here the refinement of dynamical fcrmions does not lead 
to significant corrections to the quenched results. But this is only 
known in hindsight and should, if possible, always be checked. 

Ideally, in the above example of meson-meson scattering, we would want to 
perform a lattice calculation with four light quarks i. e. where the quarks in 
{qq){qq) are m, d -quarks with masses of less than 10 MeV. Unfortunately, 
with present day computers, this is not yet possible. Therefore, the problem 
must be simplified — a process that can be done in several stages by making 
more and more of the quarks infinitely heavy i.e. q Q. This makes the 
lattice calculations easier and easier. Examples are as follows: 

• [(Qq)(Qq)] configurations. The energy of this system can be ex- 
pressed as V{R), where R is the distance between the two static 
quarks [Q). Ground and excited state energies in V{R) can then 
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be calculated in the Born-Oppenheimer approximation by assum- 
ing that the Q's have some definite finite mass. A specific example 
would be the B — B system, where mg is the mass of the 6-quark 
i.e. about 5 GeV. This will be discussed in more detail later in 

sec.inni 

• [(Qq)(Qq)] configurations. These are very similar to the above but 
with the added feature that we can now have the qq annihilating 
to give simply a (QQ) configuration. Such a process is a model 
for the string breaking mechanism QQ{R) — s- {Qq){Qq), where for 
some sufficiently large R it becomes energetically favourable for a 
gg-pair to be created from the vacuum. This mechanism must occur 
in nature, but at present it has only been conclusively demonstrated 
in simplified versions of QCD. A specific example of this type of 
configuration would be the J5B-system — to be discussed later in 
Sec. lTTTI 

• [(QQ)(QQ)] configurations. These are the most simple four-quark 
configurations to deal with on the lattice. Unfortunately, they 
are also the most distant from appropriate experimental data — 
the nearest being T(55, 9.5GeV)T(55, 9.5GeV) scattering. However, 
this class of experiments is still far in the future. 

Each of these three possibilities will be discussed separately in the following 
sections. 

The problem we are faced with is that the Lattice and Nuclear Physics 
approaches are extremes that concentrate on two opposite aspects. In the 
Nuclear Physics approach only the quark degrees of freedom are introduced 
with the gluon degrees of freedom entering only implicitly in the interquark 
potential. In contrast, with the [(QQ)(QQ)] configurations — the ones 
studied here in most detail — the lattice calculations treat the quarks in 
the static quenched limit, where they play no dynamical role — all the 
effort being to deal explicitly with the gluon field. 



1.4 (QQ) and [(QQ)(QQ)] Configurations 

The most important observable for QQ states is the interquark potential 
Vqq , which is the topic of Chapter 3 — see also Ref. .SS^ for a more detailed 
account. Here only those aspects of the {QQ)~system will be mentioned 
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that are relevant to the later [{QQ){QQ)] discussion — the main interest 
in this section. 

When lattice calculations on four-quark systems were first attempted, 
it was not clear whether an acceptable signal could be achieved. Therefore, 
various simplifications were made in an attempt to maximize the possibility 
of success |52I53I86I87I88| . The most important of these was the use of 
quarks with only two colours i.e. SU(2). In addition, only configurations 
where the four quarks were at the corners of a square were studied, since 
it was hoped that this degenerate situation would lead to the maximum 
interaction — a feature that in fact turned out to be so and also one that 
had been assumed in the flip-flop model of Ref. (HH! ■ 

1.4.1 Lattice calculations with (QQ) configurations 

Schematically, in analogy to thermodynamics, the expectation value of some 
operator O can be obtained with an action S through the chain 



where, at the first stage, the variables U are continuous. However, for 
QCD, even though we know the exact form of the action S{U), the result- 
ing integrals are singular. As discussed earlier, it is necessary to carry out 
step 1, in which the variables U are discretized into variables Ui that sit 
on the links (i) of a lattice. This now removes the singularities in a con- 
sistent manner and reduces the problem to a numerical approximation for 
a multiple integral. Unfortunately, there are very many t/^'s. For example, 
in QCD the Ui are in fact colour matrices C/f — defined by 8 independent 
real parameters. Furthermore, for a 10^ space-time lattice there are about 
4 X 10* links. In all this amounts to about 320000 integrations to be done. 
If each of these were approximated by 10 points, we would then be approx- 
imating the multiple integral with about lO^^oooo terms — not a feasible 
task. However, many of the link configurations lead to large values of S{Ui) 
and so would have a negligible effect. These configurations can be avoided 
by "importance sampling" , which essentially generates configurations with 
the probability given by the Boltzmann factor exp[—S{Ui)]. This automat- 
ically encodes the feature that the ratio of the values of nearby links are 
exponential in form. Now only a relatively few configurations (often ~ lO's) 



(1.10) 
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need to be generated to get a good estimate of (O). This is depicted as 
step 2 in Eq. 11.101 where the Bohzmann factor has been replaced by the 
N]^J configurations that tend to maximize this factor — see, for example, 
p. 252 in Ref. [T3|. 

The problem is, therefore, reduced to generating these "important" con- 
figurations and to finding an appropriate operator. These are the topics of 
the next two subsections. 

1.4.1.1 Generating lattice configurations 

Earlier, practitioners of lattice QCD generated their own lattice configura- 
tions. However, nowadays there are groups that specialize in this very time 
consuming computer task, e.g. the UKQCD group in Edinburgh. The rest 
of us arc then able to simply perform expectation values of our operators 
in the knowledge that we are using configurations that are well tested and 
essentially independent of each other. This also has the added benefit that 
different groups use the same configurations, i.e. with the same lattice 
parameters, but with different operators. At times, this can lead to useful 
direct comparisons of the lattice results with experiment that do not have to 
be scaled in any way beforehand, i. e. if one group has extracted the lattice 
spacing for their particular problem, then the rest of the lattice community 
can use this same lattice spacing to convert all of their dimensionless lattice 
results directly into physical quantities. 

In spite of this modern trend, it is probably useful to remind the reader 
of some of the problems that arise when generating configurations: 
Equilibration 

If a lattice simulation is ever started from scratch, then the first lattice 
must be simply "guessed" . This could be "cold" , where all the Ui are taken 
to be the same, or "hot" , where the Ui are random numbers. To carry out 
the above "importance sampling" these lattices must first be equilibrated 
or, in the terminology of spins-on-a-lattice, thermalized. This is where 
clever programming enters to generate quickly configurations that are as 
independent of each other as is conveniently possible. For example, in our 
early work |8f)j , where we generated all of our own configurations, the lattice 
was equilibrated by the heat bath method. 
Updating 

Once the lattice is equilibrated, measurements of the operator O can begin. 
However, the lattice must be continuously updated after each measurement 
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to ensure these measurements are made on lattices that are, as far as pos- 
sible and convenient, independent of each other. In Eg. ll.lOl the number 
of such lattices is referred to as Nm- As an example, in Rcf. 53 — after 
equilibration — the lattice was further updated by a combination of three 
over— relsLxation sweeps, each of which can change the configuration a 
lot but leave the energy unchanged. This is followed by one heat bath 
sweep, which changes the configuration less but can change the energy as 
is necessary for ergodicity. Then a measurement is made of the appropriate 
correlation functions. In that particular calculation, the lattice was updated 
6400 times with 1600 measurements being made. The latter were divided 
into 20 blocks of 80 measurements for convenience in carrying out an error 
analysis. In later work |59| 40 updates were made between measurements 
to reduce possible correlations between successive measurements. However, 
it should be added that, in our experience, there has never been any sign 
of a significant correlation due to insufficient updating between measure- 
ments, when dealing with the types of problem we have been studying. 
This is usually checked by calculating the autocorrelation between blocks 
of measurements. In the above, I mention spins-on-a-lattice — the classi- 
cal example being the Ising model — in order to remind the reader that 
lattice QCD and the thermodynamics of spins-on-a-lattice have very much 
in common. Over the years, many ideas were first developed in the spin 
case before attempting to implement them in the more complicated case of 
lattice QCD. 

1.4.1.2 Appropriate operators on a lattice 

Most of the operators evaluated on a lattice take the form of correlations 
between different Euclidean times. The most simple example is the Wilson 
loop involved in extracting Vqq in Eqs. 11.61 and 11.71 Consider a Q and 
a Q are at the lattice sites x and x -f 1 — with the notation that 1 is the 
lower case letter corresponding to the upper case L, which is reserved for 
the spatial size of the lattice. A QQ(l)-sta.te, is then constructed from a 
sequence of connected lattice links Ui, at a fixed Euclidean time ti, as 



Here the colour indices a,b,c, . . . on successive links are coupled to give 
overall a gauge invariant chain of J7's. Also the Q and Q have the same 
colour to ensure the meson is a colour singlet. To be more specific, if the 



*(0g, 1) = 0^(x)[/"''[/^"[/"'* . . . C/^''0J(x + 1). 



(1.11) 
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Q and Q are a distance of two links apart i.e. they are at lattice sites 
{x, 0, 0, ti) and {x + 2, 0, 0, ti), then one possibility is the direct path 



hp'^'ix.x + l;ti)U'"'{x + l,x + 2;ti)<P''nix + 2, ti), (1.12) 



where the two [/-links are simply along the a;-axis. Other choices of less 
direct paths between the Q and Q are possible and, indeed, necessary if 
excited states are of interest. To construct a Wilson loop a similar wave 
function is written down at time t2 = ti + 1 as 

^^{x, t2W'^'{x,x + l;t2)U^'^'{x + l,x + 2;t2WQ{x + 2, ta). (1.13) 

This gives the two horizontal (wavy) sides of the Wilson loop in Fig. 11.11 
To complete the loop we need to insert the propagators of the Q and Q 



Tt 



ti 

X x+1 



Fig. 1.1 The Wilson loop for a QQ state ^[Q(x)Q(x + 1)] propagating from Euclidean 
time ti to t2. 

from ti to t2. For static quarks these are, for t2 — ti — 2, 

(/.^(a;;ii)C/'^'''(a;;ii,ti + l)C/^''''(x; ti + 1, ii + 2)^^' (x; is) 

and 



0|(x + 2; ti)U"-'''{x + 2; ii, ti + I jU^'""' {x + 2; ii + 1, h + 2)^ {x + 2; ts), 

where here the C/'s are all in the T direction and, without loss of generality, 
can be arranged by gauge fixing to be simply unity. The Wilson loop of 
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Eg 1 1.71 then reduces to the overlap 

W{x, x + 2; h, t2) = S^^-'V^ix, x + l, ti)U'"'{x + 1, x + 2, ti)x 

U"'''' (x, x + 1, t2)U^'''' (x + 1, X + 2, ia)] ^ (1.14) 

This is simply a number, which in principle — > exp[—tVQQ{l = 2)] as t — » oo. 
The emergence of an exponential factor from this overlap should not be 
surprising, since the essence of the importance sampling in Subsec. 11.4.11 
was to encode the Boltzmann factor in Ea. ll.lOl into the values of the links. 
Of course, a single overlap W^(x, x + 2; ^1,^2) would not be very informative. 
Only when — for all three orientations x,y,z — this is averaged over the 
whole L^T lattice and the Nm different lattices does one expect a reasonable 
numerical signal for W(2,t) to emerge. In general, this has the form 

L T 

Nm lo=xo,vo-zo ti 

^ exp[-tVQQ{l)] as t^oo, (1.15) 
where the average over the different lattices includes: 

(1) All values of x in the range < x < L, 

(2) All three spatial directions x, y, z 

(3) All values of i\ in the range < ti < T 

and it is repeated for many as-independent-as-possiblc lattices Nm- 

Having introduced the notion of a lattice link Ui representing the gluon 
field, we can combine four links to form a closed loop called an elemen- 
tary plaquette [/□. For example, a plaquette in the xy-plane would be 
constructed from links \J(x x' : y ^ y' ) as 

U^y = U'"^{x ^ X + 1 : y)U''''{x + l:y^y+l)x 

V^ix + l-fx:y + l)[/'"'(x -.y + l^y). (1.16) 

Similarly there can be plaquettes in the xz, . . . , zt-planes. It was in terms 
of these plaquettes that Wilson 47 first expressed the discretized form of 
the QCD action mentioned in Sec. ll.l.Tl which for SU(iVc) has the form 

5w = ^E[l-l^TrC/a]. (1.17) 
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Usually the basic coupling (g) is expressed in terms oi f3 — 2Nc/g'^- 
1.4.1.3 Fuzzing 

The above Wilson loops were constructed from the basic lattice links U^{n), 
where — using a slight change of notation for convenience — the earlier 
colour indices a,b, . . . are omitted and replaced by an index j = 0, n is a 
lattice site and /i is either the x,y or z direction. These basic (j = 0) links 
can now be generalised by "fuzzing, blocking or smearing" (901 91| by means 
of which the basic link is supplemented by a combination of neighbouring 
links. Here I will concentrate on the the fuzzing option, but blocking and 
smearing are similar. Fuzzing is depicted in Fig. 11.21 



If 

2 



1=1 

* 2i 



I 

1 



b) 



Fig. 1.2 Fuzzing. In a) the solid line 12-link is replaced by the dashed line 12-link. 
In b) this dashed line 12-link is itself now replaced by the wavy line 12-link. 



This illustrates the replacement 



(1.18) 



to give links U'^ constructed from the basic links by a single fuzzing 
J = 1. This is followed by 



(1.19) 



to give links constructed from the by a second fuzzing j = 2 and so 
on. In general 



(1.20) 
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Here the are normalisation factors chosen to project the U^{n) into 
SU{Nc) and c is, in principle, a free parameter. However, experience when 
measuring energies 90 91 has shown for SU(2) that c = 4 is a suitable value 
for the present class of problems. The value of c could be optimized, but it 
is found that results are never crucially dependent on the precise value of c. 
Another degree of freedom is the amount of fuzzing (m) . For correlations 
over large distances the greater the fuzzing the better the efficiency of the 
calculation in the sense that the state ^(QQ) in Eq. 11.111 generated by 
connecting together a series of fuzzed links, is expected to have a greater 
overlap with the true ground state wave function. In some cases this has 
been carried to great lengths e.g. in Ref. [^2] m =110 fuzzing iterations 
were used, since there the emphasis was on quark separations upto 24 lattice 
spacings. 

In the calculation of the interquark potential in Eq. 11.151 the fuzzing 
procedure plays the second role of generating different paths {Pi) between 
quarks, for use in the variational approach of Ref. |90l 91| . In this notation 

the basic path from Ea. lTTTl is Pq = [/O.afc^o.bc^^o.cd ^o,za fuzzed 

paths = ui,abJJ^,bcu^,cd JJ^,za ^ rpj^^ variational approach then yields 

the correlation matrix between paths with different fuzzing levels i, j 

=< P,\f'\Pj > . (1.21) 

Here T ~ exp{~aH) is the transfer matrix for a single time step a, with 
the basic QCD Hamiltonian H, the Pij are paths constructed as prod- 
ucts of fuzzed basic links and, as before, t is the number of steps in the 
imaginary time direction. As shown in Ref. [901 91| . a trial wave function 
ip — J2i 0'i\Pi > leads to the eigenvalue equation 

I^l^.a* = A(*)l4^*,ria*.. (1.22) 
For a single path this reduces to 

where Vqq is the potential of the quark system being studied. Unfortu- 
nately, in this single path case i.e. with only Pq as in the previous subsec- 
tion, t needs to be large and this can lead to unacceptable error bars on the 
value of Vqq extracted. However, if — in addition — a few fuzzed paths 
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are included, it is found that t need only be small (t « 5) to get a good 
convergence to Vqq with small error bars. 

A further very important advantage of fuzzing is that not only can the 
lowest eigenvalue be extracted but also higher ones, since a matrix diago- 
nalisation is involved. These higher states correspond to excitations of the 
gluon field. However, with the direct paths from the Q to Q in Eq. 11.121 
these excitations are purely S-wave. To generate non-S-wave gluonic exci- 
tations combinations of indirect paths are needed — see Chapter 2. 

The above fuzzing is an attempt to improve the overlap of the lattice 
wavefunction with the true ground state wavefunction and so only involves 
spatial links (N.B. v ^ i m the summations in Eg. 11.181 i.e. the staples in 
Fig. II. 2| are purely spatial). However, a similar procedure can be applied 
to a time-like link In the notation of Eq. 11.181 the basic time- like 

link U^^Q^n) is now replaced by VF^^qC*^) — ^^^^ average of the six staples 
in the planes {zkx,T), {zLy,T) and (±z,r). The benefit gained from this 
is a reduction in the Noise/Signal ratio i?NS for configurations involving 
a static quark. Without this time-like fuzzing, it has been observed for 
i3-meson correlations, with time extent xq, that i?NS oc exp(a;oA£^), where 
AE = Eo ~ rriTr with Eq being the ground state energy of the B- meson [SJ . 
This leads to very noisy signals as xq increases. However, using W^^Q(n) 
results in almost an order of magnitude reduction in i?NS for a;o ~ 1.5 fm. 
This is shown in Fig. 11.31 

An even greater reduction in Rns can be achieved, if the single staple 
fuzzing in Fig. ll.2l is replaced by the hypercubic blocking of Ref. There 
each level of fuzzing is described by three equations (and so needing three 
free parameters Ci) similar to Eg. 11.181 But these three equations mix links 
that are connected to the original link only i.e. this "fuzzing" only involves 
links in the hypercube defined by the original link. This procedure yields 
a fuzzed link Whyp that is much more local than three fuzzing levels from 
Eq. 11.181 and so preserves short distance spatial structure. Furthermore, 
when this is applied to a time-like link, then it reduces i?NS even further 
than the use of W°^q (n) . The use of hypercubic links is still in its infancy 
and we should expect to see much more of this development in the future 
— see Refs. [HE]- This is also shown in Fig. 11.31 
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Fig. 1.3 Noise to signal ratio (-Rns) of a B— meson correlation function as a function of 
the time extent xq. The solid dots are when using simply (70^(,(n). The open squares 
are for iy^_Q(n) and the open triangles for Whyp I!J3I . 



1.4.2 Lattice calculations with [{QQ){QQ)] configurations 

The previous section outlines the techniques for extracting the two-quark 
potential. This is now generahsed to the [{QQ) (QQ)] case, where the quarks 
are on the corners of a rectangle with sides of length d and r. During the 
Monte Carlo simulation, the correlations Wij (r) in Eq. 11.211 — appropri- 
ate for extracting the two-quark potential VQQ{r) — and the correlations 
Wij{d,r) for the four-quark potential V4((i, r) are evaluated at the same 
time. As an example, in Ref. ^6 for Wij{r) three paths were generated by 
the three different fuzzing levels rn =12, 16 and 20. On the other hand, for 
the Wi'j'{d,r) only the level 20 was kept, with the variational basis {i'j') 
being the two configurations A and B in Figs. 11.41 These lead to the Wilson 
loops in Fig. 11.51 In the case of squares {i.e. r — d), Eg. 11.221 gives the 
potential energy of the ground state as 



Eo ^aV4,{d,d) 



log 



aE^id, d) - 2E{d) 



'-\d,d) + wi^\d,dy 



Wl,{d,d) + Wl2{d,d) 



- 2E{d) 



(1.24) 
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and that of the first excited state 



El = log 



Wl^\d,d) 



wliid,d) 



Wl^id,d) 



2E{d), 



(f.25) 



where E{d) is the result of diagonalising the 3x3 variational basis for Vqq. 
In these equations, t should tend to oo. However, in practice it is found to 
be sufficient to have t « 5 for extracting accurate values for Eq and Ei . 

The results are shown in Fig. 11.61 For comparison this also includes the 
corresponding 1985 results of Ref. [HZl as the two points with large error 
bars at r/a « 1.25 and 2.5 — remembering that the lattice spacing for 
Ref. inZ] is a' « 0.15fm « 1.25a, where a « 0.12 fm from Ref. It is 

seen that the present calculation gives energies that have error bars which 
only become significant for r/a > 7, whereas the 1985 work was unable to 
generate any meaningful numbers beyond r/a' = 2. In fact, for the present 
purpose, neither of the 1985 points is of use, since the r/a' = 1 result could 



t 

d 



(a) 



(c) 



Fig. 1.4 a) Four quarks in a rectangle of sides r and d: 

The two partitions h) A = [QiQ3][Q2Q4] and c) B = [QiQ4][Q2Q3] 




'1: "12 

Fig. 1.5 The 4-Q Wilson loops Wii and W12. 
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Fig. 1.6 The comparison of the lattice data Eq^i (continuous hue) with the / = 1 
model results Eq (dashed line) in Subsec. il.5.11 — see Ref. |87I . Also the two lattice 
data points from Ref. i97i are included. 

well suffer from lattice artefacts — a problem afflicting all calculations of 
configurations involving a single lattice spacing — and the r/a' = 2 result 
has error bars that are too large for the analysis in Subsec. to be carried 
out. However, the following should be added in defence of Ref. 97 . Firstly, 
their calculation was for SU(3) and so it was, at least, an order of magnitude 
more demanding on CPU time — especially in 1985. Secondly, it should 
be remembered that the energies of interest (Eq i) are very small compared 
with the total four-quark energy E4 in Eg. 11.241 For example, at r = d = Aa 
the value of uEq = —0.050(1) is obtained from the difference between aE^ 
and 2aE{d), which is 1.505-1.555 i.e. the error quoted on aEo is less than 
0.1% of aEi. 

One of the most outstanding features of the results is that for r equals d 
the value of aE^ decreases smoothly in magnitude from -0.07 for r — d — a 
to -0.04 for r = d = 6a. However, for the few cases where r ^ d the 
value of Eq is an order of magnitude smaller than the adjacent r = d cases 
e.g. Eq{2,3) w -0.006, whereas £^0(2, 2) and £^0(3, 3) are w -0.055. This 
result is reminiscent of the conclusion found with the flux-tube model 
and the ansatz made in the flip-flop model of Refs. [HHI- In both of these 
models the interaction between the two separate two-quark partitions of 
Fig. 11.41 is very small (in fact zero in the flip-flop model) except when the 
unperturbed energies of the two partitions is the same i.e. at r = d in the 
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case of rectangles. 

1.4.3 Lattice parameters and finite size/scaling check 

Most of the above (and later) calculations involving only static quarks were 
carried out with the parameters in Table 11.21 There the coupling /3 is 
defined after Eq. 11.171 and the lattice spacings a„ were determined by, for 
example, the Sommer method described in Subsec. 11.2.1.21 Each of these 
three sets serve a specific purpose. 

1.4.3.1 Benchmark data 

Set 1 shows the original parameters with which most calculations are per- 
formed {e.g. for the results in Fig. I1.6|l and against which the following 
results with Sets 2 and 3 are compared. These parameters are chosen so 
that — unlike Sets 2 and 3 — the problems of computer memory and 
CPU time do not prevent generating many "independent" measurements 
to ensure good error estimates. 

1.4.3.2 Finite size effect 

With Set 2 the effect of the finite lattice size is checked. However, in the 
present case it was found that the results for the size of squares considered 
{i.e. up to 7 X 7) were unchanged within error bars — see Table 2 in 
Ref. |HH|- This type of check is important for large squares because of 
the spatial periodic boundary conditions at and L. Two quarks can 
be connected by two paths — a direct path on the lattice {e.g. that in 
Eq. 11.12(1 or an indirect path that encircles the boundary. If the two quarks 

Table 1.2 Typical SU(2) lattice parameters used in [{QQ){QQ)\ calculations. Much of 
the notation is explained in the text with M/G being the number of Measurements per 
Geometry. 



Set 


xT 




a„ (fm) 


M/G 


Refs. 


1 


16^ X 32 


2.4 


« 0.12 


720 


[Hni-EHl 


2 


243 X 32 


2.4 


« 0.12 


160 


w 


3 


243 X 32 


2.5 


« 0.082 


660 


88^ 
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are separated by r > L/2, then the indirect path is of length L — r < r 
and so is potentially more important than the direct path — the only one 
of the two that is explicitly treated. In the present case with L = 16, 
such problems should only begin to occur seriously with 8x8 squares. For 
this larger lattice, fewer configurations per square were generated, since the 
storage and computer time increased by about a factor of « (24/16)'^ ~ 3. 

1.4.3.3 Smaller lattice spacing 

In most cases the possible importance of finite size effects can be seen 
beforehand by using simple arguments. We know, for example, that the 
length r of the direct path should be smaller than the length L — r oi 
the indirect path. However, the effect as a —> is not obvious and needs 
checking. In order to isolate this scaling effect from the finite size effect, 
it is convenient to use lattices that have approximately the same physical 
size. Therefore, since the spatial volume for Set 1 is (16 x 0.12)'^ « 7 fm^, 
the value of 03 for a 24^ lattice should be « 7-^/^/24 « 0.08 fm — a lattice 
spacing that corresponds to /3 « 2.5. For the results to show scaling, the 
physical energies {i.e. in, say, MeV) at the same values of r (in fm) should 
be the same. In this case, it means 



where the a„E(Set n, r/a„) are the dimensionless numbers given by the 
lattice calculations and the r/a„ are the number of lattice links between 
the quarks. Of course, in general only r/ai or r/a^ is an integer, so that the 
comparison in Eq. 11.261 must be done by interpolation. This approximate 
equality is sufficiently well satisfied in the present problem (see Table 3 in 
Ref. [HH|)- In Ref. 49 a more complete test of the a ^ limit is carried 
out with the addition of the non-rectangular geometries to be discussed 
in Sec. 11.61 There the f3 values (and lattice sizes) used were (3= 2.35, 2.4 
(16^ X 32), 2.45 (203 x 32), 2.5 (24^ x 32) and 2.55 (26^ x 32) and led to 
4-quark energies that were essentially independent of a over this range of /3 
i.e. scaling was achieved for /3 > 2.35. 

It should be added that the above procedure is called scaling in contrast 
to asymptotic scaling. The latter relates results from different values of a(/3) 
by perturhative arguments and is not expected to apply at the comparatively 
large values of a in Table 11.21 — see Ref. . 



E[^ei 1 (or 2), r/ai] w E[?,et 3, r/ag]. 



(1.26) 
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1.5 Potential Model Description of the Lattice Data 

The assumption often made by those who create models for multi-quark 
systems is that these systems can be treated in terms of two-body poten- 
tials. This is the Nuclear Physics inspired approach that is very successful 
for, say, multi-nucleon systems, where three- and four-body forces are small. 
Of course, the fact that the latter were small was at first simply an assump- 
tion. However, later this was found to be justified by Weinberg [23] using 
various low-energy theorems that force the ttN — > irN interaction — the 
main mechanism for multi-nucleon forces — to essentially vanish in nuclei. 
At present, there seems to be no such simplifying feature for multi-quark 
interactions. In the following, I will first show the consequences of using 

— in the four-quark system — simply two-quark potentials unmodified by 
the presence of the other two quarks. Then this model will be improved by 
including also a direct four-quark interaction. These are examples of what 
was referred to as Effective Potential Theories (EPTs) in Subsec. 11.1.2.21 
However, since the quarks are now static, there is no kinetic energy and so 
the question of which quark wave equation to use does not arise. Similarly, 
the effective quark masses do not enter, since the quantities of interest are 
the binding energies between the two 2-quark clusters. Therefore, of the 
three ingredients usually needed for an EPT of Subsec . [1.1.2.21 onlv freedom 
with the interquark potential remains. This is both a weakness, by being a 
model that is unrealistic and not comparable with any experimental data 

— except for possible future T{bb, 9.5GeV)T(66, 9.5GeV) scattering — and 
a strength, in that only the potential ingredient plays a role. 



1.5.1 Unmodified two-body approach 

Once the quark-quark potential Vqq is known and, if it is assumed to be 
the only interaction between the quarks, then the energy of a multi-quark 
system can be readily calculated — provided the wave function for that 
system is expressed in terms of a sufficient number of basis states. For the 
present situation, the most obvious choices for such states are A and B in 
Figs. II. 41 b"! and c). 

In this extreme two-body approach, since the presence of the gluon 
fields have been explicitly removed — their only effect now being in the 
colour indices of the quarks — the states A and B form a complete but 
non-orthogonal basis. This implicitly assumes in Figs. 11.41 the quark (Q), 
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antiquark {Q) assignment A — [QiQ3][Q2Q4]- In SU(2) the other quark 
assignment [QiQ3\[Q2Qi\ is numerically equivalent and so leads to nothing 
new. In this approach, states are excluded in which the gluon fields are 
excited. However, later this restriction is relaxed in an extension to the 
model in Sec. 11.71 and Appendix A. 3. The energies (i?-) of this static four- 
quark system can be extracted from the eigenvalues of the Hamiltonian 



(V-A,N)*, =0, (1.27) 
where the normalisation matrix 

N^(/ '(M (1-28) 



1/2 1 



and the potential energy matrix 



Wl3 + V2i Vab 
VbA Wi4 + V2Z 

Several points need explaining in these equations. 



(1.29) 



(1) The off-diagonal matrix element N12 = (^l-S) =1/2 shows the 
non-orthogonality of the A, B basis. In general, for S\]{Nc) this 
becomes N12 = 1/Nc and arises simply by recoupling the colour 
components of the quark terms i.e. 



II14I23) = ^^IllsW + V l3^24), (1.30) 

where |ljj) and \A^j) are the SU(Afc') singlet and adjoint represen- 
tations. 

At this stage, the lack of orthogonality could have been avoided by 
simply using the basis A± B. However, later it will be seen that 
the A, B basis is in fact more convenient and suggestive, when the 
gluon fields are reintroduced in a more explicit manner. 
The interquark potential VQQ{ij) = Vij in Eg. 11.61 has been ex- 
tracted as the potential energy between a single static quark (Q) 
and a single static antiquark (Q). In order to evaluate general 
multiquark potential energy matrix elements, a further assumption 
is needed concerning the colour structure of Vij. Here the usual 
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identification 




(1.31) 



will be made, where the are the Pauli spin matrices appropriate 
for SU(2). This choice ensures, for a colour singlet meson- like state 



Strictly speaking, the form in Ea. ll.3l1 is only true in the weak cou- 
pling limit of one-gluon exchange, since this has replaced the local 
gauge invariance — ensured by the series of J7- links of Ea. ll.lll con- 
necting the two quarks — by the global gauge invariance reflected 
by the Tj.Tj factor. 
(3) With the choice of Vij in Eq. 11.311 the off-diagonal potential matrix 
element becomes 



Since the following discussion only involves quark configurations in 
the rectangular geometries of Figs. 11.41 it is convenient to use the 
notation 

Vl3=V24: = Vd, Vi4=V23=Vr and Vi2 = V34 = v^, (1-34) 

where the suffix x refers to the diagonals of the rectangles in Fig. II. 41 

Even though the form of Eq. 11.331 is derived in the one-gluon exchange 
limit, it is now assumed that a more realistic model emerges if the Wy are 
taken to be the complete potential of Eq. 11.61 and not just the one-gluon 
exchange component. This clearly has the correct form when the distance 
between the two two-quark clusters of Fig. 11.41 are far apart, since in this 
case the only interactions are those within the separate clusters — due to 
the cancellation in Eg. 11.331 of vi2, with either U13, 1124 or D14, ^23. 

The 2x2 Hamiltonian of Eq. 11.271 is easily diagonalised to give the 
eigenvalues Ao,i. Since it is the binding energy E'^ of the four-quark system 
that is of interest, and also that was extracted in Subsec. 11.4.21 from the 
Monte Carlo simulation, the internal energy of the meson-like state with 



[zjf, that([zjf|T,.r,|[zjf) = -3 and 



= Vi 



(1.32) 



{A\V\B) ^Vab = Vba 

^\ («13 + 1'24 + Wl4 + W23 - V12 - V34) ■ 



(1.33) 
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the lowest energy (i.e. 2vd for d < r) is now subtracted from the A,j to give 

El = \,-2vd (1.35) 

in analogy with the lattice expressions in Eas. ll.2j| and n~^ 

Therefore, in this simplest version of the two-body approach, the E^^s 
should correspond to the Ei 's from the Monte Carlo simulation — a compar- 
ison which is made in Fig. II. 61 Since the values of Eq for squares (i.e. r = d) 
are much larger than those for neighbouring rectangles, only the results for 
squares are shown in Fig. 11.61 In these cases 

E'o = -^{v, - Vd) and E[ = 2{v, - v^) = -SE'^, (1.36) 

In addition, the corresponding normalised wave functions are of the form 
^o) = 7^\A + B) and i^{E[) ^\A-B). 

Two comments should be made on these results: 

• For the smallest squares and rectangles the agreement between Ei 
and El is best. This is reasonable, since it is expected that at 
small interquark distances perturbation theory is adequate i.e. the 
lowest order gluonic effects are already incorporated correctly into 
the interquark potential. 

• As the squares and rectangles get larger the differences between the 
Ei and the E'^ grow until Eq is more than three times Eq, and E[ 
more than seven times Ei, for the largest squares d « 7a « 0.8 fm. 
Since Eq is too attractive and E[ too repulsive, this suggests that 
the off-diagonal matrix element Vab in Eg. 11.331 is too large. 

The main conclusion to be drawn from the above is that the two-body 
potential of Eq. 11.311 does not give the potential energy of the four-quark 
system — the indication being that the off-diagonal potential energy Vab 
is too large. This is nothing more than the well known van der Waals effect 
of Ref . |1UU| — see Oka and Yazaki in Chapter 6 of Ref. . 

1.5.2 The effect of multiquark interactions 

In the above model it is assumed that all of the gluonic effects are in- 
corporated into the two-body potential Vij. However, this is clearly an 
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oversimplification that is at best only applicable in situations where per- 
turbation theory holds, namely at short distances, as already noted in the 
discussion of Fig. 11.61 In more realistic models, the QCD coupling is suf- 
ficiently strong to constrain the gluon field into flux-tubes connecting the 
quarks in a given meson — as visualised by the wavy lines between the 
quarks in states A and B in Figs. 11.41 Therefore, the overlap of states A 
and B, i.e. N12 = {A\B), is not simply the colour recoupling factor of 1/2, 
but should also involve the lack of overlap of the gluon fields in states A 
and B. This can be incorporated by introducing an entity /, which sim- 
ply multiplies the original N12, and which is an unknown function of the 
position coordinates of the four quarks. With this interpretation of / as a 
gluon field overlap factor, it is also necessary to multiply the off-diagonal 
potential matrix element Vab of Eg. 11.331 bv the same factor /. This fac- 
tor must be the same in N12 — f/2 and V12 = /Vab, since otherwise the 
binding energies E'q ^ would be dependent on the self-energy term c in the 
form of Vij given in Eq. 11.61 — which would be unphysical. The E'q are 
now extracted from Eq. 11.351 after diagonalising 

[V(/)-A,(/)N(/)]vI/, = (1.37) 

with 

N(/) = ( 1 ) and V(/) = ( ) . (1.38) 

V //2 1 / V J^BA VU + V23 ) 

The two equations H1.37|l and H1.38|l are the basis of the following anal- 
ysis. They give a procedure for extending the model of Eas. II. 271 - 11. 291 
which was justified in the weak coupling limit, into the domain beyond 
one-gluon exchange. The off-diagonal potential matrix element performs 
this extension in two ways. Firstly, even though Vab in Eg. 11.331 still has 
the same algebraic structure in terms of the as dictated by the one- 
gluon exchange limit, the 's themselves are the full two-quark potential 
of Eq. 11.61 Secondly, in the off-diagonal correlation Wyi of Fig. 11.51 the 
one-gluon exchange model suggests the presence of the overall multiplica- 
tive factor / due to gluon exchange within the initial and final states A 
and B at euclidean times T = and cxd. In this interpretation, the terms 
in Vab arise during the period of propagation between T = and 00. 

The strategy is now to adjust / to get an exact fit to Eq or E\. In the 
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case of squares, since 



-f 



[vx - Vd) and E[ 



f 



^0 1 + f/2 
the appropriate expressions are 



1-//2 



f{Eo) = 







and f{Ei) 



{Vx - Vd), 



E, 



'Vd + Vx + Ei/2 



45 



(1.39) 



(1.40) 



for fitting Eq and Ei respectively — the results being shown in Fie. II. 71 For 
general rectangles the corresponding equations are somewhat more compli- 
cated. 
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(1.1) (1,2) (2,2) (2.3) (3,3) (3.4) (4,4) (5,5) (6,6) (7,7) 
Fig. 1.7 The values of f{Eo), f(Ei) and / (from Subsec. lLOl . 

Two points should be noted from this figure: 

• All values of / are less than unity as is expected from the interpre- 
tation of / as a gluon field overlap factor. In addition, this idea is 
supported by the fact that the values of / decrease as the quarks 
get further apart. 

• The values of f{Eo) and f{Ei) are rather similar, which suggests 
that a compromise value of / would give a reasonable description 
of both Eq and Ei — as will be seen below. 
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1.5.3 A compromise for the overlap factor f 

The previous subsection suggests that for each geometry a single value of 
/ could give a reasonable description of both Eq and Ei . Here one such 
possibility is given by finding that value of / which minimizes the expression 



where the AEi are the errors quoted for the Ei from the lattice calculation. 
The result is shown by the dashed line in Fig. 11.71 

It is seen that indeed a single value oi f — f suffices to explain reason- 
ably well both energies. This is a non-trivial observation, since it indicates 
that the parametrization suggested in Eqs. 11.371 11.381 contains the most 
important features of the more precise lattice calculation. 

It should be added that the extraction of a compromise value of / is 
not simply a curiosity, since any model that needs different values of / for 
Eq and El would be more difficult to use in practice for more complicated 
multi-quark systems. 

At this point, even though for each geometry a single value of / = / 
gives values of Eq i that are in reasonable agreement with the -Eo.i j it might 
be asked about the remaining small differences. Several possibilities are now 
open: 

• The lattice energies Ei may not be sufficiently accurate due to finite 
lattice size and scaling uncertainties. This was checked in R.ef . |88| 
and found not to be a problem. 

• The parametrization in Eqs. 11.371 11.381 may be inadequate. One 
possibility would be to combine the notion of a gluon overlap factor 
/ with a generalized form for the two-quark potential in Eq. 11.61 
since this could introduce more free parameters. 

• Any model based only on states A and B in Fig. 11.41 is incom- 
plete and other states are necessary in addition. This point will be 
discussed further in Sec. 11.71 

• The step from Eq. 11.271 to Eq. 11.371 was motivated by one gluon 
exchange. However, this could possibly be extended by performing 
a two-gluon exchange calculation to see what new terms arise and 
to then be guided by this in making an improved parametrization 
— the topic of the next subsection. 




2 



2 



(1.41) 
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1.5.4 The effect of two-gluon exchange 

In Ref. it was noted that the two-body / = 1 models discussed in 

Subsec. Il . 5 . II correspond to lowest order perturbation theory in the quark- 
gluon coupling i.e. to order a = 5^/47r in the notation of Ea. IT.17l This was 
extended in Ref. |l()2j . where a perturbative calculation to fourth order in 
the quark-gluon coupling (i.e. to 0(a^)) was made for the potential of the 
QQQQ system. This was performed in the general case of colour SU(Af). 
Considering the quarks to be at points Ri define 

Va = V{Ri3) + V{R2i), Vb = V{Ru) + F(i?23), Vc = V{Ri2) + V{R^i), 

where V{Rpq) is the two-body potential between the different quark and an- 
tiquark pairs a distance Rpq apart. For the two possible energy eigenstates, 
diagonalization yields the following two potentials correct to 0{a^): 

(iV2_2) {Va + Vb) + 2Vc - N^N^{VA-VBf ^-AiVA- Vc) {Vb - Vc) 

2{N-^-l) 

Vi = 

{N^~2) [Va + Vb) + 2Vc + N"^ [Va-Vb)'' + A [Va - Vc) [Vb ~ Vc) 

2 (iV2 _ 1) ■ 

(1.42) 

These potentials are exactly equal to those given by the naive (/ = 1) two- 
body model in Eas. 11.271 - 11.291 This fact that a straightforward two-body 
model is correct also to next-to-leading order in the quark-gluon coupling 
may be surprising in light of the non-abelian nature of QCD. However, 
Ref. |l()2j does go on to show that this two-body model fails at 0[a^) as 
three- and four-body forces appear due to the onset of three-gluon vertex 
effects. In general, their nature seems to be complicated, but for some 
geometries simplifications are possible; e.g. for the four quarks on the 
corners of a regular tetrahedron there will be no contribution from quark 
self-interactions to four-body forces to 0[a^). 

The overall conclusion from Ref. |102| is as follows: "Looking at the 
Monte Carlo lattice calculations for the QQQQ-system in Refs. [871 88| . 
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it is observed that for small interquark distances of a few lattice spacings 
(with a w 0.12 fm) the / = 1 two-body model gives a reasonable ap- 
proximation in the sense that the four-quark potentials calculated from 
Eg. 11.421 using the Monte Carlo two-body potentials are comparable to the 
four-quark potentials from the lattice simulation. The agreement improves 
the smaller the distances get. By comparing the perturbative (i.e. 1/R) 
and non-perturbative {i.e. linear) part in the usual parametrization of the 
QQ-potential in Eg. II. 61 one would expect to start entering the perturbative 
regime when distances get down to about two lattice spacings. However, at 
that stage the approximation provided by the two-body model is already 
very good. The fact that the two-body model is correct to fourth order 
in perturbation theory certainly suggests that it should be a reasonable 
approximation in the perturbative domain. This result supports the belief 
that the results of the lattice simulations for small enough distances indeed 
are correlated to continuum perturbation theory, and thus that continuum 
physics is extracted from the Monte Carlo calculations." 

So far in this section various models have been proposed in an attempt to 
understand the results of the Monte Carlo simulations of lattice QCD. The 
main outcome — summarised in Fig. 11.71 — is the emergence of a function 
/ that depends on the coordinates of the four quarks involved. This shows 
that the usual models based on purely two-quark interactions need to be 
modified considerably — essentially by the factor /, which becomes <C 1 
for large interquark distances. This observation is in itself of much interest, 
but at this stage it is not clear how the effect can be incorporated into 
more realistic situations in which the quarks are not so restricted in their 
geometry. It is the purpose of the next section to tackle this problem by 
first studying how / can be parametrized. 



1.5.5 Parametrizations of the gluon-field overlap factor f 

In Subsecs. 11.5.21 11.5.31 models were introduced in an attempt to under- 
stand the ground state binding energy (Eq) and excited state energy (£^1) 
emerging from a Monte Carlo simulation, in which four quarks were at the 
corners of a rectangle. These models are summarised bv Eos. 11.371 II. 381 For 
each quark configuration, both of the energies i?o,i are described in terms 
of a function / of the four quark positions — Eg. 11.391 As it stands, this 
is not particularly useful, when wishing to extend these ideas to systems 
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containing more than four quarks, unless / can be parametrized in some 
sensible and convenient manner. In the literature, several parametrizations 
have been suggested. In Refs. |l()ll motivated by strong coupling ar- 
guments, the phenomenological form is taken to be 



where bs — 0.0736 is the string tension in the interquark potential of Eg. II. 61 
and S is the minimal area of the surface bounded by the straight lines 
connecting the quarks and antiquarks. The other form, the one proposed 
in HOI, is 



i.e. the cut-down is governed by the average of the six links present in a 
Q^Q^ system. In Eos. II. 4^11 II. 441 the a and k are at present free parameters 
to be determined later. Both of these parametrizations of / accommodate 
the following two extreme models. 

• Weak coupling, which has f — I when all — 0. This / = 1 limit 
was assumed to apply for all values of in Subsec. 11.5.11 

• Strong coupling, which has / = when any — > oo. In this limit 
the flux-tubes between the quarks in the separate mesons — seen 
in Fig. II. 4| — are then very narrow and straight. In this case the 
flux-tube overlap of configurations A and B tends to zero in the 
limit that any r^j — > oo. 

With squares {r = d), for which the most accurate values of / exist, k equals 
3a/4. One measure of how meaningful these parametrizations really are, is 
given by extracting a and k for each quark configuration. The hope would 
then be that, for squares, a (and therefore k) would be independent of 
the separate configurations. Only the few points for non-square rectangles 
(r ^ d) would be able to distinguish between /i and /2. The outcome is 
depicted in Fig. 11.81 for the values of k{EQ), k{Ei) and k corresponding to 
the values of f{Eo)J{Ei) and / in Fig. O This shows that k{EQ) for 
the squares appears to decrease slowly from about 0.7 to about 0.5 as the 
sizes of the squares increase from (2 x 2) to (7x7), whereas for non-squares 
k{Eo) appears to be stable at about 0.7±0.1. On the other hand, k{Ei) 
and k decrease somewhat less and also the square and non-square values 
are consistent with each other. 



/i = exp[~absS], 



(1.43) 




(1.44) 
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Fig. 1.8 The values of fc(i5o); k(E\) and k corresponding to the values of f{Eo), f{Ei) 
and / in Fig. ILTI 

The indications from this are that the single parameter k « 0.5(1) is a 
suitable compromise value that results in a reasonable fit to both Eq and 
El for a series of square and near-square geometries. 

1.5.5.1 A reason for fi = exp[—absS] 
In the limit of large I, Eg. ll.7l becomes 

W{l,t) ^ eyip[~tb'J]:^exp[-bsS] as < ^ oo, (1.45) 

where S is the minimal space-time area enclosed by the loop in Fig. 11.11 
For this QQ case, the meaning of 5* is clear — it is simply a^lt. But for 
the [QQ][QQ] case in Fig. II. 51 the form of the appropriate minimal space- 
time area is less clear. However, in the extreme strong coupling limit the 
area S is the one produced by the minimum number of elementary plaque- 
ttes needed to tile the enclosed area in question. Any fluctuations about 
this space-time surface would need more plaquettes and so be higher or- 
der in the strong coupling model. In this limit, the diagonal loops Wa in 
Fig. II. 51 are again simply two Wilson loops each tiled separately by the min- 
imum number of plaquettes. However, the off-diagonal loops Wij are more 
complicated. A model for this was suggested in Ref. |1()5| and developed 
in Refs. \UYM lOfillOTj . In the notation of Fig. O the Euclidean Green's 
function for a Qa, Qb, Qci Qd system can be thought of as a 2 x 2 matrix 
for a two channel problem with the {AB){CD) and {AD){CB) configura- 
tions. The transition potential between these two configurations may then 
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Fig. 1.9 Tiling the transition between the states [{QaQb){QcQd)] and 



[(QaQdKQcQb)]- 

be extracted from the expression for the Wilson foop of this system, 



where Li and L2 are the minimum flux tube lengths in the basis states 
and bs is the string tension of Subsec. 11.2.1.11 The diagonal terms are 
each simply the product of two propagators — exp[— fe^LiT] for Wu and 
exp[— 65X2^'] for W22- In a language more familiar in nuclear physics, these 
are simply the Green's functions 



Wn = G^s{T)Gcd{T) and W22 = G^B(T)GcB(r). (1.47) 



corresponds to tiling of the off-diagonal loops Wij . This now resembles the 
standard expression for the lowest order transition of a Green's function 



where here i, f denote the initial and final channels [{QaQb){QcQd)] and 
[iQAQD){QcQB)]- Remembering that a^^ J dt we can identify the 




(1.46) 



The mixing term 




(1.48) 




(1.49) 
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transition potential as simply 



Vfi = exp[-bsLiL2]/a exp[-~bsS]/a 



(1.50) 



i.e. the transition potential can be expressed in terms of a spatial area — as 
anticipated by the models in the previous subsections. In more complicated 
geometries this should be the minimum area in coordinate space associated 
with the given boundary conditions. 

1.6 More Complicated [{QQ){QQ)] Geometries 

So far the only [{QQ){QQ)] geometries considered above were squares and 
rectangles. This suggested that the strongest interaction between two sepa- 
rate two-quark clusters occurs when the clusters are degenerate in energy — 
namely for square geometries compared with rectangles. To test this more, 
in Ref. 52 the six different geometries in Fig. II. lOl were studied. Since this 
confirmed that configurations degenerate in energy always gave the largest 
binding — not just in the square versus rectangle case — in Ref. j4H| the 
study concentrated only on those configurations that were near degenerate. 

In Ref. 53 the case of tetrahedral configurations was considered in 
some detail, since in SU(2) this has the three degenerate partitions seen in 
Fig. ll.lUf T) when d = r. Since tetrahedral and linear configurations have 
certain interesting features, they will be discussed separately below. 

However, real life is an average over all spatial configurations, so why 
should we be so interested in such special geometries for four static quarks? 
Firstly, it must be remembered that model (bridge) builders should consider 
all spatial possibilities — a failure with one configuration indicating that 
the proposed model is faulty. Also, in some ways the tetrahedral and linear 
configurations have simplicities and symmetries not present in other cases, 
which could make "Lattice QCD <-> Model" comparisons easier. 

1.6.1 Tetrahedral configurations on a lattice 

Since tetrahedral configurations are so symmetrical, at first sight there 
is no reason to consider only two of the three possible cluster partions 
A + B, B + C oi A + C. In this case, using the notation that the suffices 
1, 2, 3 arc any combination of A, B, C that forms a basis, then the appro- 
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Fig. 1.10 Six four-quark geometries a) Squares (S) and Rectangles (R), b) Tilted rect- 
angles (TR), c) Linear (L), d) Quadrilateral (Q), e) Non-Planar (NP) and f) Tetrahedra 
(T). Their energies are calculated using Lattice QCD in Refs. 1491 521531861871881 and 
analysed in Ref. |108( . 

priate Wilson loop matrix is 



W* = 







wis 






wis 






Wis 



(1.51) 



For regular tetrahedra, the general symmetries W^^ = 



Wis, Wl, = 

WI2, Wli = wis, WI2 = wis ^''"S expected and, in addition, there are 
the equalities Wis — ^12 ^^id VK^g — —Wis- Therefore, in all, there are 
only two independent Wilson loops Wl^ and Wi*2- Here the minus sign 
appearing in the last equation is a reminder that the quarks are in fact 
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fermions even though quarks and antiquarks transform in the same way 
under SU(2). This detail is discussed more in the Appendix to Ref. |52| . 

Unfortunately, the inclusion of all three partitions does not lead to even 
more binding. However, it shows the curious feature that the ground and 
first excited states become degenerate in this highly symmetrical limit, since 
the eigenvalue equation (discussed earlier as Eq. I1.22|l 

W,*,.a* = AWw^*,-ia*, (1.52) 

is easily solved to give for the lowest energy (occurring twice) 

]/{/t^ _|_ yy*2 

Ai 2 = — ; — 1 ; — T (1.53) 

and for the excited state 

In comparison, using only two partitions gives 

W*=fW^n ^1*0 (155) 

where not only is the general symmetry PV^*2 — ^21 expected but also for a 
regular tetrahedron Wl^ = W^22- I'^ this case Ea. ll.52| is again easily solved 
to give for the lowest energy 

i.e. the same as Ai in Ea. ll.53| However, for the energy of the first excited 
state 

which is quite different to the complete result in Eqs. 11.531 11.541 So the 
degeneracy is easily explained as a feature of the more complete 3x3 lattice 
QCD simulation. The effect is depicted in Fig. 11.111 as d ^ r. This also 
shows that, whereas ground state energies can be quite stable, those of 
excited states are more model dependent. 

The dominance of the cluster interaction by degenerate configurations 
was carried to the extreme in the so-called flip-flop model of Ref. |89| . 
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111 
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-0.1 

12 3 4 

r/a 

Fig. 1.11 The binding energies - in units of the lattice spacing — of the four-quark states 
for the tetrahedral geometry of Fig. IlToI fT) for d/a = 3 and r/a = 0, 1, 2, 3, 4. 
Sohd hnes show lattice results: o — the ground-state binding energy Eq. 
X — the first excited state energy Ei . 

Dashed and dotted lines show model results with / = 1 from Egs. 11.271 - 11.2!^ for Eq 
and El respectively — see Ref I53| . 

which makes the ansatz that the 4-quark interaction occurs only in the 
case of exact degeneracy. However, the model is developed for a purely 
linear interquark potential i. e. with only the bs term in Eq. 11.61 which 
means that the interaction only occurs when — in the notation of Fig. 11.41 
— the spatial distances ri3 -I- r24 and ri4 + r23 are equal. 

In several of these geometries it is not clear what is the best "natural" 
partition into two-quark clusters. Therefore, in the lattice simulation all 
three possibilities A, B and C in Fig. 11.101 should be taken into account. 
However, in some cases it is found that one (or more) of the combinations 
A + Bot:A-\-CotB + C\s sufficient to give — for both Eq and Ei — 
similar results to the complete A -\- B -\- C simulation. This is particularly 
true if the configuration with the lowest energy is one of the configurations 
used. The study of tetrahedral and near-tetrahedral geometries on a lattice 
can be summarized as follows: 
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• From earlier work in Rcfs. 52 8 61 871 88| . for the corresponding 
squares [i.e. (d, 0) in Fig. ll.lOl (T)]. where two of the basis states 
are degenerate, the binding energy of the lowest state ranges from - 
0.07 to -0.05 'AS d/a goes from 1 to 5. However, now — even though 
at least two of the basis states are always degenerate {A and B) — 
the ground state binding energy [Eo{d, r)] is always less than that 
of the corresponding square [Eo{d,0)]. For a fixed d, \Eo{d,r)\ de- 
creases as r increases from r — 0. Nothing interesting happens to 
Eq a.t r = d, the point at which all the basis states are degenerate 
in energy. 

• For fixed d, as r increases from to d, the energy of the first excited 
state El decreases until Ei{d,d) = Eo{d,d). For r > d, Ei{d,r) 
increases again. This degeneracy of i?o,i for the tetrahedron is a 
new feature compared with earlier geometries. As will emerge in 
Sec. 11.71 this is a severe constraint on models. 

• The choice of which 2x2 basis to use depends on the particular 
geometry, because one of these two basis states must have the lowest 
unperturbed energy. Since A and B are degenerate in energy, for a 
given d this amounts to using A + B for r < d. On the other hand, 
for r > d it is necessary to use B + C, since C now has the lowest 
unperturbed energy. 

• Except for the tetrahedra, the values of Ei are essentially the same 
in the 2x2 and 3x3 bases. 

• The values of E2{d, r) are always much higher than Ei{d, r). How- 
ever, as will be discussed in the Sec. 11.71 this second excited state 
is dominated by excitations of the gluon field and so is outside the 
scope of the models introduced in that section. 

The above is for colour SU(2). However, in the real world of SU(3) the 
state C — [QiQ2\[Q3Qi\ cannot appear asymptotically as two clusters — 
see Eg. I A. 31 Even so this SU(2) lattice data should be understandable in 
terms of models that are similar to those of SU(3) — see Subsec. 11.2.21 



1.6.2 QCD in two dimensions (1+1) 

The above analyses unavoidably involve numerical results extracted from 
Monte Carlo simulations on a 4-dimensional lattice. However, it is possi- 
ble to study colinear colour sources in a simple approximation for which 
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exact theoretical results are known This is QCD in two dimensions 

(QCD2) - the one spatial direction allowing colinear (and only colinear) 
configurations. For quenched QCD2, the spectrum is known exactly even 
on a lattice. This can be summarised as the requirement that each link is in 
a representation of the colour group [i.e. SU(2) here]. Therefore, the links 
can be in the singlet ground state (J=0), or they can be excited to a funda- 
mental representation (J=l/2) or an adjoint representation (J=l) and so 
on. The energy per unit length for such links is given by bj — ^KJ{J+ 1). 
For 4 colinear quarks, the lowest state [A in Fig. ll.lOl fLl] has energy 2Kd, 
while the first excited state has, in the middle, an adjoint link of length 
(r — d) resulting in an energy 2Kd + 8K{r — d)/3. This exact result applies 
at any /3- value, i.e. strong coupling or weak coupling. 

A comparison with the above /-mixing model of Fas. II .371 II .381 shows 
that there is agreement with the exact QCD2 results provided / = 1 — 
independent of which combination of states is used for the analysis A+B, 
A-l-C or B-l-C. This is also true in the more general case when di3 7^ ^24. If 
the interpretation of / being a gluon-ficld overlap factor is correct, then it 
is easy to understand that here / must be unity, since with only one spatial 
direction the colour flux must overlap fully. Farlier, the main motivation 
for the /-mixing model had been from weak coupling arguments, so that 
this agreement between the mixing model and QCD2 does suggest that the 
model is sensible even at large /3. This adds support to the claim that 
it may be a useful phenomcnological tool when the extension to / 7^ 1 is 
made. 



1.7 Extensions of the 2 X 2 /-Model 

In Sec. II. 51 a model for describing 4-quark interactions was developed as a 
2x2 matrix equation for the two heavy-quark states A — [QiQs] [Q2Q4] 
and B = [QiQ4][Q2Q3] in Fig. O— see Eas. [07I and [QHl There the 
only geometries studied were those in which the quarks were at the corners 
of squares or rectangles. However, when the four quarks are in more general 
geometries, such as those in Fig. ll.lOl the choice of which two partitions out 
of the possible three to use is less clear i.e. C = [QiQ2][Q3Q4] should be 
included. This is completely analogous to the problem that arose in Sub- 
sec. when deciding which configurations to use in the corresponding 
lattice calculations. 
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In Appendix A and in Ref. |1U9| the model of Sec. II. 51 is extended from 
being a 2 x 2 matrix equation for states A, B into a 3 x 3 version, where 
all 3 basic partitions A, B, C are included. Finally a 6 x 6 extension is 
developed, in which interquark excited states are introduced to give three 
additional basis states A*, B* , C* . 

Much of the discussion will concentrate on the problems that arise in 
trying to model the energies of regular tetrahedral configurations. This 
might be considered a minor point to worry about, since such a configura- 
tion is very special. However, the philosophy is that, if any configuration 
cannot be fitted, then the model fails, since then there is no reason to expect 
other configurations not checked explicitly to be fitted. 

A summary of the main points that emerge from Appendix A are as 
follows: 

f) The 2x2 — > 3x3 extension, to some extent, clarifies the reason 
why the earlier 2x2 f-model was, in many cases, quite successful. Also it 
shows that an understanding of the tetrahedron spectrum — in particular 
the degenerate ground state — requires a generalisation of the two-quark 
approach. 

2) The step 3x3 — > 6x6 has the very positive feature for tetrahedra 
of giving a ground state binding energy that initially increases with the size 
of the tetrahedron — a result also seen in the corresponding lattice data. 
This arises naturally, since the energies of the additional states A* , B* , C* 
are each excited by an energy of tt/R with respect to the A, B, C states 
in the 3x3 model. Here R is the interquark distance for that pair of 
quarks containing the excited gluon field. Therefore, as the four-quark 
configuration gets larger spatially, the energies of the A*, B*, C* states 
approach from above the energies of the A, B, C states. The subsequent 
mixing between the two sets of states then manifests itself as an additional 
overall attraction that also grows with the spatial size. 

3) An even more interesting conclusion from the 6x6 extension is that 
it partially justifies the / = 1 model of Subsec. If .5.fl In Appendix A it is 
seen that, by fitting simultaneously the lattice energies for the geometries 
in Fig. If .101 the 6x6 model shows two features: 

• Outside the range where perturbation theory holds (i.e. beyond 
about 0.2 fm) the binding is dominated by the A*, B*, C* config- 
urations. 

• The overlap factor between the A*, B*, C* states (corresponding 
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to / in Eg. ll.38l for the A, B, C states) is essentially unity. 

We, therefore, come to the following scenario for the four-quark interac- 
tion: At the shortest distances, up to about 0.2 fm, perturbation theory 
is reasonable with the binding being given mainly by the A, B, C states 
interacting simply through the two-quark potentials with little effect from 
four-quark potentials. However, for intermediate ranges, from about 0.2 to 
0.5 fm, the four-quark potentials act in such a way as to reduce the effect 
of the A,B,C states so that the binding is dominated by the A*,B*,C* 
states, which now interact amongst themselves again simply through the 
two-quark potentials with little effect from four- quark potentials. This sug- 
gests that models involving only two-quark potentials could be justified — 
provided excited gluon states (such as A* , B* , C* ) are included on the same 
footing as the standard states A, B, C. The above result that excited states 
play an important role in the overall binding is reminiscent of the nucleon- 
nucleon interaction, where nucleon excitations — especially the A(1236) — 
are also responsible for a sizeable part of the attraction. 

1.8 Heavy-Light Mesons (Qq) 

In the previous sections only static (i.e. infinitely heavy) quarks with two 
colours were discussed. Even though this is a far cry from the real world 
of finite mass quarks with three colours, it resulted in several interesting 
conclusions. However, these were of a somewhat academic nature useful 
for creating models, but could not be compared directly with experimental 
data. In this section a compromise situation of two quarks is studied where 
one of the quarks is still static but with the other being light. This is es- 
sentially the "hydrogen atom" of quark physics and is expected to be a rea- 
sonable representation of the heavy-light i?-mesons. This also means that 
in the development of the Effective Potential Theories of Subsec. [1.1.2.21 it 
is possible that the one-body Schrodinger or Dirac equation is applicable. If 
this proves to be so, then it will be a further reason for studying _B-mesons 
in addition to the more basic ones discussed below. 

It should be added that the hydrogen atom analogy also partially holds 
for the interaction — the coulomb potential oc e^/r of the hydrogen atom 
versus the one-gluon exchange cx a/r in the heavy-light meson. However, 
as mentioned in Subsec. 11.1.2.21 for light quarks there are indications that 
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this 1/r attraction gets damped by form factor effects and that much of 
the needed attraction could arise from a short ranged instanton-generated 
interaction. Also, it must not be forgotten that beyond r «0.2~0.3 fm the 
interaction in Ea. ll.fil becomes dominated by the linear confining potential. 

From Eg. 11.21 for lattice calculations with a heavy quark Q, we had the 
condition that a should satisfy a <C Wg^. At present, this rules out direct 
lattice calculations with b quarks of mass 5 GeV, since they would require 
a 0.04 fm. This prompted the weaker condition a ^ '^q^ i^i Eq. 11.31 
One way to partially avoid this problem is to perform lattice calculations 
for lighter quarks with mg 2 GeV, which do not require such a fine 
mesh, and then to extrapolate the results to the h quark mass. However, 
if the results from a siatic-light system are also included in the analysis 
we end up having to interpolate to mg ~ 5 GeV. Since heavy quarks with 
such masses are essentially non-relativistic, the appropriate interpolation 
parameter is 1/toq, so that the actual interpolation is between 1/mg = 
to 0.5 GeV~^. This is usually a more reliable procedure |110| . 

1.8.1 Bottom {B)-mesons 

Bottom mesons are the bound qq states of a 5-antiquark of mass ~5 GeV 
and a lighter quark. These are listed in Table II!^ This century has opened 
with there being renewed interest in S-physics — see the Bottom meson 
summary in the Particle Data listings of Ref. |111| . There are new gen- 
erations of i3-meson experiments at BaBaR (SLAG), Belle, CLEG III and 
Hera-B. These machines have started to accumulate i?-mesons and the long- 
awaited B-factory era has begun. The hope is that these experiments will 
deliver the fundamental constants of the Standard Model and also improve 
our understanding of CP violation. However, having seen that _B-mesons 
are important objects, it must be confessed that their study by lattice QCD 
is very incomplete. This will be the topic of the present section, where we 
concentrate on their energies, and in Sec. 11.91 where we extract density 
distributions. 

1.8.2 Lattice parameters 

In the past few years there have been several detailed measurements of 
i3-meson excited state energies. Some of the parameters used in these 
studies are given in Table [TM The quenched calculations of Refs. |56l 58) 
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Table 1.3 Properties of the Bottom mesons lllll . The state marked with * is from 
Ref. 11121 and those marked with ? do not, at the time of writing, have I{J^) confirmed. 



Meson 


qb 


i{.n 


nL 


Mass(GeV) 


B+ 


ub 




15 


5.2790(5) 




db 




15 


5.2794(5) 


b'" 


db 




25 


5.859(15)* 


B*+ 


ub 




15 


5.325(6) 




ub 


i(0,l,2+) 


IP 


5.698(8) ? 




sb 


o(o-) 


15 


5.3696(24) 


b; 


sb 


o(i-) 


15 


5.4166(35) ? 




sb 


0(0,1,2+) 


IP 


5.853(15) ? 




cb 


o(o-) 


15 


6.4(4) 



Table 1.4 Lattice parameters. The notation is explained in the text. Refs. 1561 581 are 
in the quenched approximation and Refs. 1591 1141 are unquenched. 



Rcf. 


L3 


X T 


/3 


Csw 


K 


a (fm) 


Mps/Mv 


m 


12^ 


X 24 


5.7 


1.57 


0.14077 


« 0.17 


0.65 


m 


16^ 


X 24 


5.7 


1.57 


0.14077 


« 0.17 


0.65 


m 


16-3 


X 24 


5.2 


1.76 


0.1395 


w 0.14 


0.72 


im 


16^ 


X 32 


5.2 


2.02 


0.1350 


w 0.11 


0.70 



are very similar to each other with the latter having a somewhat larger 
lattice. The parameter Csw has been tuned for the clover improved action 
mentioned earlier j5l 6j . The hopping parameter k essentially determines 
the mass ruq of the light quark as being slightly smaller than the accepted 
value of the strange quark mass i.e. rriq = 0.91(2)tos. The ratio of the 
pseudoscalar to vector masses {Mps /My) — i.e. the ratio of the "tt"- and 
p-meson masses generated by those particular configurations — is another 
measure of mq. Since this ratio is much larger than the experimental value 
of 0.18, we again see that iriq ^ rriu^d- Refs. |59l 114| are expected to 
be a distinct improvement, since these are unquenched calculations with 
smaller lattice spacings. Now rriq — 1.28, 1.12ms respectively. However, in 
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practice, it is found that the energy sphttings of the excited states are only 
weakly dependent on mq for the range of values used here. 

1.8.3 Maximal Variance Reduction (MVR) 

One of the reasons why the energies of Qg-states can now be calculated 
reliably is not only the improvement in computer capabilities, but also 
by developments in formalism. It has been demonstrated that light-quark 
propagators can be constructed in an efficient way using the so-called Max- 
imal Variance Reduction (MVR) method. Since this has been explained in 
detail elsewhere, for example Subsec. 4.5.3 of Chapter 4 and in Ref. |5Hj . 
the emphasis here will be mainly on the differences that arise when esti- 
mating on a lattice the two- and three-point correlation functions C2, C3 
needed for measuring spatial charge and matter densities (C3) in addition 
to the energies (C2). In the MVR method the inverse of a positive definite 
matrix A is expressed in the form of a Monte Carlo integration 

j Dcj) 0*0, exp(-i0M0), (1.58) 

where the scalar fields </> are pseudofcrmions located on lattice sites z, j. For 
a given gauge configuration on this lattice, N independent samples of the (p 
fields can be constructed by Monte Carlo techniques resulting in a stochastic 
estimate of as an average of these N samples i.e. Aj-^ — {4>*4>j). The 
N samples of the fields can be calculated separately and stored for use 
in any problem involving light quarks with the same gauge configurations. 
In practice, « 25 is found to be sufficient. 

In LQCD the matrix of interest is the Wilson-Dirac matrix Q — 1~kM, 
where M is a discretized form of the Dirac operator + m) and is the 
mechanism for "hopping " the quarks from one site to another. However, 
Q is not positive definite for those values of the hopping parameter k that 
are of interest. Therefore, we must deal with A = Q^Q, which is positive 
definite. Since M contains only nearest neighbour interactions, A — with at 
most next-to-nearest neighbour interactions — is still sufficiently local for 
effective updating schemes to be implemented. In this case the light-quark 
propagator from site i to site j is expressed as 

Gq = G,, = Qt/ = ((Q,fc0fc)*0,) = (V*0,>. (1.59) 

This is the key element in the following formalism. The Wilson-Dirac ma- 
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trix also leads to an alternative form for the above light-quark propagator 
from site i to near site j 



In practice both forms are used, since — for the same correlation — they 
lead to independent measurements, which can then be averaged to improve 
the overall statistical error. Later, it will be essential to use at some lattice 
sites operators that are purely local. This then restricts us to using at such 
sites only the (j) fields that are located on single lattice sites. In contrast the 

fields, defined above as Qikfl'k, are not purely local, since they contain 
(j) fields on next-to-nearest neighbour sites. 

In the above, the term "Maximal Variance Reduction" comes from the 
technique applied to reduce the statistical noise in Eq. 11.591 The lattice is 
divided into two boxes {0 < t < T/2 and T/2 < t < T) whose boundary 
is kept fixed. Variance of the pseudofermionic fields is then reduced by nu- 
merically solving the equation of motion inside each box. This allows the 
variance of propagators from one box to the other to be greatly reduced. 
However, in the case of a three-point correlation in Subsec. 1 1.91 two prop- 
agators are needed and this is best treated by choosing one of the points 
to be on the boundary of the boxes while the other two are inside their 
own boxes. Furthermore, the field at the boundary must be local to avoid 
the two propagators interfering with each other. This means that only the 
4> fields should be used on the boundary and there they can couple to the 
charge, matter or any other one-body operator. For the points in the boxes, 
the temporal distances from the boundary should be approximately equal 
to give the propagators a similar degree of statistical variance. 

1.8.4 Energies of heavy-light mesons (Qq) 
1.8.4.1 Two-point correlation functions C2 

The basic entities for measuring energies are the two-point correlation func- 
tions C2 . These are depicted in Fig. 11.121 and are seen to be constructed 
from essentially two quantities — the heavy-quark (static-quark) prop- 
agator Gq and the light quark propagator Gq. As discussed in detail in 
Ref. |2HIj when the heavy quark propagates from site (x, t) to site (x, t-l-T), 



(1.60) 
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75 75 75$* 75 $ 75 75*1' 




7b * 75 7b* 7b ** 7b 7b** 



a) b) c) d) 

Fig. 1.12 The four contributions to the two-point correlation function d- 

Gq can be expressed as 

Gq(x, t ; X, i + T) = 1(1 + 74)C/^(x, t, T), (1.61) 

where U^{x,t,T) — YIJSq U4{x,t + i) is the gauge hnk product in the time 
direction. On the other hand, as the hght-quark propagates from site i 
to site j, it can be schematically expressed as one of the two alternatives 
{Gq or G'g) in Eas. ll.St^l and ll.GUI Knowing Gq and Gq, the general form of 
a two-point correlation can be constructed from a heavy quark propagating 
from site (x, t) to site (x', t + T) and a light quark propagating from site 
(x', t + T) back to site (x, t) as 

G2{T) = Tr(rt Gq(x, i; x', t + T)T Gq{-x! , t + T;x, t)) 

=2(Re [U^ [?A*(x, t + T)<P{x, t) + 0* (x, t + T)'tl;{x, t)]]). (1.62) 

Here F is the spin structure of the heavy quark - light quark vertices at t 
and t + T. In this case F = 75 , since we are only interested in pseudoscalar 
mesons such as the _B-meson. For clarity, the Dirac indices have been 
omitted. The four contributions to G2 are depicted in Fig. 11.121 Here the 
a) term uses the light quark propagator Gq in Eq. 11.591 and term b) the 
alternative G'^ in Eg. II. 601 — the two terms in Ea. ll.62| Terms c) and d) are 
the corresponding ones for a heavy antiquark (Q). It is necessary to include 
the Q-terms to ensure C2 is real. It would be sufficient to use only a)-|-c) 
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or b)+d), since both combinations correspond to measuring C2. However, 
since these two combinations are independent measurements of the same 
correlation, keeping both improves the statistics on the final measurement. 

The above has been written down for a single type of gauge field. The 
correlations can now be greatly improved by fuzzing as discussed in Sub- 
sec. 11.4.1.31 This makes the two-point correlation function into the fuzzing 
matrix C2,ij- Since i, j usually take on 2 or 3 values, this means that S-wave 
excited state energies and properties can now be studied in addition to those 
of the ground state. 

1.8.4.2 Analysis of C2 to extract energies 

There are several ways of analysing the correlation functions C2 in order 
to extract the quantities of interest i.e. energies. For a review of these 
methods see Ref. |113| — with more details using the present notation 
being found in Ref. j58j . 

The actual analysis gives not only the energies {rria) but also the eigen- 
vectors (v) for the states of the Qg-system. These values of rUa and v" are 
later fixed when analysing the three-point correlation data C3 to give the 
charge and matter densities x°'^{r) in Sec. 11.91 Each element C2,ij{T) is 
fitted with the form 

Ah 

C2^^J{T) « C2^^J{T) = ^ < exp(-TO„r)v;, (1.63) 

a=l 

where M2 is the number of eigenvalues to be extracted (usually 3 or 4) and 
mi is the ground state energy of the heavy-light meson. The values of TOq, 
and vfj are then determined by minimizing the difference between the C2 
data from the lattice and the parametric form C2. The function actually 
minimized is the usual 



E E 



C2,ij{T) — C2.ij{T) 

AC2,.,(T) 



(1.64) 



where AC2,ij{T) is the statistical error on C2,ij{T) and T2,min, 72, max are 
the minimum and maximum values of T2 used in the fit. The latter depend 
on the lattice size and the future use to which the TOq and vf^ are destined. 
A typical outcome is depicted in Fig. 11.131 from Ref. |114| . This is for a 
dynamical fermion calculation and a corresponding quenched calculation on 
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Fig. 1.13 The energies of Qq states from Ref. 11141 . The solid line uses dynamical 
fermions and the dashed line is a corresponding quenched calculation. The energies of 
the L-wave excited states and the S-wave radial excited state are relative to the ground 
state (IS) — see text for notation. These energies are given both in terms of GeV — right 
axis — and in the more usual Sommer units of Rq Ri 0.5 fm defined in Subsec. 11.2.1.21 
— left axis. 



a 16^ X 24 lattice. In both cases, for numerical reasons, the light quark has 
a mass that is approximately that of a strange quark. Since the heavy quark 
is static, the energies can be labelled by L±, where the coupling of the light 
quark spin to the orbital angular momentum gives j ~ L ± 1/2. The total 
angular momentum (J) is then obtained by coupling j to the heavy quark 
spin giving J = j ± 1/2. However, since the heavy quark spin interaction 
can be neglected, the latter two states are degenerate in energy i.e. the P_ 
state will have — 0+, 1+ and the P+ state will have = l+,2+ etc. 
The _D-waves show the interesting feature that there appears to be little 
or no spin-orbit splitting between the D- and £*+ states — contrary to 
some expectations |115l llfi) that there should be an inversion of the level 
ordering (with lighter than L_) at larger L or for radial excitations. 
This has important implications for phenomenological interpretations of 
the data. We return to this in Subsec. 11.9.3.21 For F-waves, only the 
energy from a spin independent mixed operator F± is shown. The latter 
is expected to approximately correspond to the usual spin-average of the 
F_ and F+ states. 
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Fig. 1.14 The two contributions to the three-point correlation function C3. 

1.9 Charge and Matter Distributions of Heavy-Light Me- 
sons (Qq) 

In many cases, when phenomenological models are constructed to describe 
lattice data, the emphasis is on fitting the energies. However, there are 
other observables that can be measured on a lattice. Of potential value 
when constructing models are lattice data for radial distributions of various 
quantities such as the charge and matter densities, which can be measured 
using three-point correlation functions C3. Here we concentrate on the 
radial distributions of the q in the Qq system, whereas in Ref. |117l 118| the 
much more ambitious task of measuring radial distributions in the tt, p, N 
and A is tackled. The reasons why distributions in few-quark systems have 
received much less attention are two-fold. Firstly, unlike energies, these 
distributions are not directly observable, but only arise in integrated forms 
such as sum rules, form factors, transition rates etc. Secondly, as will be 
seen later, their measurement on a lattice is more difficult and less accurate 
than the corresponding energies. In spite of this, it is of interest to extract 
lattice estimates of various spatial distributions. 

1.9.1 Three-point correlation functions C3 

When the light-quark field is probed by an operator 0(r) at i = as the 
heavy quark propagates from t — — 12 to t = ii, the result is the three-point 
correlation function depicted in Fig. 11.141 This involves two light-quark 
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propagators — one (G^) going from t = —t2 to t — and a second (Gq) 
going from t = to t = ti and has the form 

C3H2, <i, r) = 

Tr(rt Gq(x, -is; X, ti)TGg{x, ii; x + r, 0) e(r)G;(x + r, 0; x, -is)), 

(1.65) 

which can be expressed in terms of the pseudofermion fields 0(x, t) and 
■0(x, t) — similar to that for G2 in Ea. ll.fi2| Here the vertex Q — J4 for the 
charge (vector) distribution and 1 for the matter (scalar) density. Again 
the Q-term in Fig. ^^31 ensures that G3 is real, but now there are only two 
terms compared with the four for C2 in Fig. 11.121 This is because the fields 
connected to the probe O must be local, since the purpose of the probe is to 
measure the charge or matter distribution at a definite point r. Therefore, 
only those light quark propagators that involve the local basic field at r 
can be used, since the tp field contains contributions from (jj fields at next- 
to-nearest neighbour sites and so is non-local. This must also be kept in 
mind when fuzzing is introduced to give a matrix Cz.ij{—i2, ti, r)- Here 
the fuzzing indices i,j refer to the various fuzzing options of the t/j's at the 
Qq vertices. As with the energies extracted from G2, the fuzzing permits 
the measurement of excited state distributions. 



1.9.2 Analysis of C3 

The analysis of the three-point correlation functions 03(0, T = ti +t2,r) 
is performed using a generalisation of the one for G2 in Eq. 11.631 namely, 
fitting G3^ij(0,T, r) with the parametric form 

M3 Ms 

G3,y(e,T,r) = ^^<exp[-m,ii]x"'^(r)exp[-m0(T-ii)]vf (1.66) 

a=l 13=1 

The ttIq and v- vectors are those obtained by minimizing the G2 in Eq. 11.681 
and, for each value of r, the x°'^{r) are varied to ensure a good fit to 
(6, T, r) by the model expression C^i^ij (8, T, r). 
Two forms of x°'^{r) have been used: 

(1) A non-separable (NS) form, where each x°'^{r) is treated as a single 
entity. 
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Fig. 1.15 The ground state charge (C) and matter (M) densities as a function 

of r/a from Ref. [591 . The lines shows a fit to these densities with a sum of two lattice 
exponential functions — see Subsec. 11.9.31 The scaled quenched results of Ref. [581 are 
also shown by filled circles and squares. 

(2) A separable (S) form x°'^{r) — ya{r)yi3{r). It is seen from Eg. II. 651 
that this appears to be a more natural parametrization, since it 
contains the product of two light-quark propagators Gq. Also the 
ijair), to some extent, resemble a wave function for the state a, 
since its square yields a distribution. 

The outcome is shown in Fig. ll.lSl for the ground state and in Fig. ll.ltil for 
the excited states. Several points are of interest in Fig. 11.151 

• At small values of r the two densities are comparable i.e. x}j ~ x\]. 

• As r increases from zero the matter density drops off faster than the 
charge density. A similar difference has also been seen in Ref. |118j , 
where the authors measure these densities for the tt, p, N and A on 
a 16^ X 32 lattice with f3 = 6.0 for both quenched and unquenched 
configurations. 

• The densities calculated with the quenched approximation in 



Ref. [SHI same, within error bars, as those for the full dy- 

namical quark calculation of Ref. 59 . However, as will be discussed 
in Subsec. 11.9.41 the matter sum rule does seem to differ. 
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Fig. 1.16 a) and b): The ratios y2{r) = x^^/x^^ and ysir) = x^^/x^^ for the charge 
distribution, c) and d): These ratios for the matter distribution — see Ref. 1591 . 

• The densities do not have a smooth variation with r, but, as will 
be shown in Subsec. 11.9.31 many of the kinks can be understood 
in terms of latticized forms of standard Yukawa, exponential or 
gaussian functions. 

In Fig. ll.lSl the excited state results are shown for different types of analysis. 
The label "sep" refers to the separable assumption for a;*-' and "non-sep" to 
the non-separable assumption. The two numbers associated with each set 
of data refer to r2,min and Ts^min in Ea. ll.64l and the corresponding equation 
for Eq. II. 661 In this figure the main point of interest is the appearance of 
nodes. The presence of these nodes is very clear and also the number of 
nodes is as expected. The first excited state has a single node at about 
0.3 fm in the charge case and about 0.2 fm for the matter. The second 
excited state seems to have two nodes with one being at about 0.1 fm and 
a second at about 0.4 fm for the charge and 0.3 fm for the matter. It will 
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be a challenge for phenomenological models to explain these data. 

When discussing the use of the separable form x°'^{r) = ya{r)yp{r), 
it was stated that Uair) can possibly be interpreted as a wave function 
for the state a. However, there are other radial distributions associated 
with the Qq system that can also be interpreted as wave functions. These 
are the Bethe-Salpeter wave functions Wa{r) discussed in Ref. •")() . They 
are extracted by assuming the hadronic operators Caa {ri , r2, T) to be of 
the form Wa{ri)wa{r2) Q'^'p{~'^aT), where the sink and source operators 
are of spatial size ri and r2. Qualitatively, the forms of yi{r) and y2{f) 
plotted in Fig. 11.161 and the corresponding ones for wi (r) and W2 {r) are 
found to be similar. However, even though they do bear some similarities, 
it should be added that there are several reasons why these two types of 
wave function should not agree in detail with each other. In particular, 
the [i(;a(r)]^ cannot be identified as a charge or matter distribution. In 
addition, they are found by using an explicit fuzzed path between Q and 
q and so are dependent on the fuzzing prescription, whereas the yaif) are 
defined in a path-independent way. 

The above considers only S-wave distributions for both ground and ex- 
cited states. The extension to other partial waves is now in progress |61| . 
Preliminary results for the P_ state indicate that there the distributions 
are qualitatively of the form expected from a Dirac equation description 
i.e. the charge and matter distributions are not zero for r = 0. 

1.9.3 Fits to the radial forms 

1.9.3.1 Fitting data with Yukawa, exponential and gaussian forms 

In Figs . 1 1 . iKl and IT. 1 61 the results are presented as a series of numbers. How- 
ever, even though these are the actual lattice data, they are, in practice, not 
very convenient to use or interpret. To overcome this it is, therefore, use- 
ful to parametrize the data in some simple form. Furthermore the results 
in Fig. 11.151 do not follow smooth curves but exhibit several kinks. If the 
latter are first ignored, then average fits can be reasonably well achieved 
with simple Yukawa (Y), exponential (E) or gaussian (G) forms giving 
X^/^doi ~ 1-4. The reason for using exponential and Yukawa radial func- 
tions is that they arise naturally as propagators in quantum field theory — 
usually in their momentum space form (g^ -I- m^)~^, where m can be inter- 
preted as the mass of a meson being exchanged between the heavy quark 
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and the point at which the light quark is probed. On the other hand, when 
going away from quantum field theory and attempting to understand the 
radial dependences in terms of wave functions from, for example, the Dirac 
equation, then gaussian forms can arise naturally — see Subsec. 11.9.3.21 
However, this fitting can be greatly improved by using lattice versions (LY, 
LE, LG) of the above Yukawa, exponential or gaussian forms, namely 



LY 



TT 



E 



cos(r.q) 



D + 0.25[a/r 



LY]2 ■ 



(1.67) 



[exp(-r/rL^)] 



LE 2r^^L-^ 



E 



cos(r.q) 



[i:i + 0.25(a/rLE)2]2' 



(1.68) 



[exp[-(r/rLG)2]] 



LG 



aL 



J2cos{r.q)exp[~{r^^/a)^D]. (1.69) 



Here L is the lattice size along one axis and D = X]i=i sin^(a<7i/2), where 

"t/z — L ' ■ • ■ ' L 

These are able to give much of the kink structure in the data — as is 
seen in Fig. 11.171 where in each case two Yukawas, exponentials or gaussians 
are used to give the fits 2LY, 2LE and 2LG with x^/n-doi ~ 1- All three 
of these forms are equally acceptable [SHI- In the 2LE and 2LY fits it is 
of interest to express the range parameters r^^ and r^^ in terms of the 
mass (1/r*) of the particle producing this range. Both types of fit find that 
the longest range corresponds to the exchange of a vector meson of mass 
niQ «1 GeV for the charge density and a scalar meson of mass mg ^1.5 GeV 
for the matter density. These masses of the mesons have been extracted 
in a rather indirect manner. However, in the literature there have been 
direct calculations of the energies of these qq states using the same lattice 
parameters and lattice size as those employed here. In Ref. |119j they 
got = 1.11(1) GeV and in Ref. [HII = 1.66(10) GeV — numbers 
consistent with the above indirect estimates. 

It should be added that the expression LY in Eq. 11.671 is often used in 
the Coulomb limit r'"^ — )■ oo for discretizing the 1/r term of VQQ{r) in 
Eg. 11.61 However, in that limit the condition q 7^ is necessary. 
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The best fits for charge and matter densities 
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Fig. 1.17 Fits to the lattice data in Fig. II. 1.^ with lattice exponential (2LE), Yukawa 
(2LY) and gaussian (2LG) forms from Ref. 59 . 

1.9.3.2 Fitting Qq data with the Dirac equation 

The above fits are purely phenomenological with there being no connection 
between the forms fitting the charge and matter distributions. However, 
there is an alternative approach in the spirit of an Effective Potential Theory 
(EPT), in which the data are fitted by the solutions of the Dirac equation. 
In the present situation using the Dirac equation and not the Schrodinger 
equation does not seem unreasonable on three counts: 

(1) The light-quark propagators are generated by a discretized form of 
the Dirac operator — see Subsec. 11.8.31 

(2) The mass (or any effective mass) of the light quark is ^ 1 GeV, so 
that relativistic effects are expected. 

(3) Figs . 1 1 . 1'Hl and ITTtI show that the charge and matter distributions 
are different — a feature not easy to understand in a non-relativistic 
Schrodinger approach, where the wave functions have only one com- 
ponent (f), so that both the charge and matter distributions would 
be proportional to In comparison the Dirac equation has two 
components. Of course, one could also resort to the argument famil- 
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iar in nuclear physics, when discussing the difference between the 
charge and matter RMS radii of nuclei, namely, that the matter 
distribution is the basic one with the charge distribution obtained 
by folding in the charge radius of an individual nucleon. We return 
to this in the next subsection. 

In the notation that G and F are the large and small components of the 
solution to the Dirac equation then the charge (C) and matter (M) distri- 
butions can be expressed as 



x'^^ir) = G^{r)Gi3{r)+F^{r)Fi3{r) and rr^f (r) = G^{r)Gp{ryFa{r)Fp{r) 



respectively. Attempts are now underway jST; to study to what extent the 
above distributions can indeed be interpreted by these two relationships. 

Since we are now in the realm of EPTs, we need the three ingredients 
discussed in Subsec. 11.1.2.21 — a wave equation (now the Dirac equation), 
a potential Vq^ and an effective quark mass (a free parameter). The main 
problem is the form of Vqq, since it is not at all clear whether the "natural" 
generalisation of the form appropriate for static quarks in Eq. 11.61 to 



is correct. Here the first term is treated as a four-vector simply by the 
analogy between one-gluon exchange and one-photon exchange. However, 
if this term for large r is in fact the so-called Liischer term |122) , a vibrating 
string correction of — tt/ (12r) to the leading string term bgr, then it should 
be a scalar. On the other hand, there are theoretical reasons that partially 
justify the use of the Dirac equation for heavy-light mesons with the form 
for the potential in Eq. I1.7UI For example, in Ref. |121| the authors derive 
a Dirac equation for a heavy-light meson by starting from the QCD La- 
grangian and taking into account both perturbative and nonperturbative 
effects. The Coulomb-like effect is treated rigorously and the confining po- 
tential heuristically. The outcome is that the confining potential is a scalar 
and the Coulomb part is the fourth component of a 4-vector. However, the 
lack of any significant D-wave spin-orbit splitting, as is seen in Fig. ll.KSI 
does suggest that the confining potential can not be purely scalar. This 
follows from the simple argument that, in a heavy(static)-light quark sys- 
tem, a central potential of the form V{r) = a/r + br should give rise to a 



VQ,{r) 



(1.70) 
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spin-orbit potential 

V (r) — 1 dV _ I (1 71) 

Am^r dr Am? r 

Such a potential would lead to inversion of the level ordering (with i+ 
lighter than i_) at larger L or for radial excitations and this is not seen 
in Fig. 11.131 The authors of Ref. in fact, give arguments why the 

interquark confining potential should have the form 

y(r) = a/r + 6r(l +74)/2 (1.72) 

leading to simply 

^so(0 = T^4' (1-73) 

which would rapidly vanish at larger L. In Ref. |123| the electromagnetic 
structure functions of a heavy(static)-light quark system interacting via the 
potential in Eg. 11.721 are calculated. This ambiguity between the vector vs 
scalar structure of Vqq is an ongoing argument — see Ref. ^1 for a list of 
references on this controversy. Also, as mentioned in Subsec. 11.1.2.21 form 
factor effects reduce the role of any one-gluon exchange and that this loss of 
vector-like attraction could be replaced by a very short ranged instanton- 
generated scalar interaction. 

An interesting property of the solutions to the Dirac equation is that, 
for a linearly rising potential, both the large (G) and small {F) components 
decay asymptotically as gaussians. This follows from the observation that 
the coupled Dirac equations for large r can be written as 

-m{r)G{r) = -F'{r) and m{r)F{r) = G'(r), (1.74) 

where m(r) = m -I- cr —> cr as r —> oo giving asymptotically the simple 
harmonic oscillator equation G" -\- {cr)^G = 0. This was the reason for 
considering not only Yukawa and exponential but also gaussian forms in 
the last subsection. 

1.9.3.3 Fitting Qq data with the Schrddinger equation 

In the previous subsection three reasons were given for attempting to fit the 
lattice data with the Dirac equation. However, two of these reasons are far 
from compelling — only the fact that the effective light-quark mass being 
<C 1 GeV seems unavoidable. The comparison with the discretized Dirac 
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operator used in the lattice simulation, as mentioned in Subsec. 11.8.31 is 
nothing more than an analogy without theoretical basis. Secondly, the fact 
that the charge distribution has a longer range than that for the matter is 
reminiscent of the difference between the charge and matter radii of nuclei. 
There this is usually expressed as 



relating the Mean Square Radius of the nuclear charge distribution with 
that for the matter distribution and the proton charge radius. This sug- 
gests that there is only one basic distribution — that of the matter — with 
the charge distribution arising from a finite size correction to the charge of 
the light-quark. This should not be surprising since the effective quark mass 
{iTT'q, effective) is quitc different (larger) compared with that used in the lat- 
tice calculation, which is about that of the strange quark. Here we are now 
saying that, in addition to a mass renormalisation, the light-quark develops 
a charge form factor. A direct application of Ea. ll.75l on the fits to the data 
in Fig. lTTTl results in q RMS charge radii of 0.51(4), 0.49(7) and 0.35(3) fm 
for the Yukawa, exponential and gaussian fits respectively. These are sizes 
consistent with ruq^ effective 500 MeV {i.e. l/rriq^ effective ~ 0.4 fm) and 
are surprisingly large being 2-4 lattice spacings. Also q charge form factors 
of this size qualitatively explain why in Fig. Il.lfil the node in the x^'^/x^^ 
charge distribution is at a larger value of r than that for the matter. 

Since the problem has been reduced to only one basic distribution, we 
are now able to use an Effective Potential Theory (EPT) based on a non- 
relativistic Schrodinger equation. As before, the other two ingredients for an 
EPT are niq, effective (a free parameter) and the potential Vgq presumably 
based on the form appropriate for static quarks in Eq. 11.61 

In Ref. (10 . when analysing the Q^(f lattice data of Ref. 77 , this 
strategy was used first to extract a value of m^. effective by fitting the spin- 
averaged energies of the Qq system — see Fig. 11.181 The outcome was 
JTiq. effective ~ 400 MeV — a valuc consistent with the above size estimate. 
We return to this in Subsec. [1.10.21 The above approach has also been car- 
ried out in Ref. |117| . where nucleon charge correlations are measured and 
then fitted with a non-relativistic Schrodinger equation. First the authors 
extract two entities C{r/\) and C(ry), which they interpret as the nucleon 
wave functions in either the coordinate ta — (7'i2 + ^23 + ^'13) — the so- 
called A-Ansatz — or the coordinate ry = Min(ri£ -f- r2e + r^^), where e 




Large — 



(f"^) matter + {r'^) 



proton 5 



(1.75) 
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Fig. 1.18 Fits to the spin-averaged energies in Fig. ll.l51 for a series of values for mq. 

is the junction at which the three flux tubes meet. They then go on to 
fit this data with Airy functions that decay as exp(— cr^''^). These are the 
wave functions expected from a non-relativistic Schrodinger equation when 
using hnearly rising potentials V(r/^) oc or ^(ry) cx ry. In this way, 
fitting the C(rA,fy) data seemed to shghtly favour the A-Ansatz. How- 
ever, it should be added that this conclusion is far from being universally 
accepted. For example, in Ref. J,24 the authors say: "In particular the 
A-shape configuration debated in the literature is shown to be impossible 
and the well-known Y-shaped baryon is the only possibility." This latter 
result is supported in Refs. |125l 126| . 

1.9.3.4 Fitting Qq data with semirelativistic equations 

The equations considered in the above subsections are two extremes with 
the Dirac equation being fully relativistic and the Schrodinger equation 
completely non-relativistic. However, if — as in the "real" _B-meson — 
the heavy quark is not static, then other possibilities arise when the basic 
Bethe-Salpeter equation is reduced in a systematic way to a Lippmann- 
Schwinger form. As mentioned in Subsec. 11.1.2. 21 this can be carried out in 
a variety of ways, which give rise to the Blankenbecler-Sugar, Gross, Kady- 
shevsky, Thompson, Erkelenz-Holinde and other equations |35| . Unfortu- 
nately, these equations are not so easy to treat because of the unavoidable 
presence of the typical relativistic factors y/mi/Ei, where Ei = y/mf+pi^ 
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is the energy of a single particle of momentum p; . This automatically leads 
to a non-local interaction in coordinate space, since a local interaction re- 
quires a function of the relative momentum q = Pa — Pb- When dealing 
with one-gluon exchange this does not present a problem, since there the 
basic interaction is oc and so the equations can be formulated directly 
in momentum space. However, as seen in Eq. 11.61 a crucial part of the 
interquark interaction is the confining term bsr — an interaction that can 
not be conveniently treated directly in momentum space. In the literature 
there are several attempts to fit directly the meagre experimental i3-meson 
data with these forms. For example, in Refs. |126| and |127| the Gross and 
Blankenbecler-Sugar equations, respectively, are employed to interpret the 
i3-meson spectrum and a series of available transition rates. 

Some of the complications that arise when dealing with the above equa- 
tions can be partially overcome by using instantaneous interactions. Shortly 
after formulating the Bethe-Salpeter (B-S) equation in 1951 jl28| . Salpeter 
in 1952 |129 | replaced the interaction G{q^ = Ql^ — oi) in the B-S equation 
by its three dimensional counterpart G(q^). This resulted in an equation 
for two particles a and b, which could be written in the centre-of-mass 
system as 

[E - Ha{p) - if6(p)]</>(p) =PJ d\74G(q)740(P + q), (1-76) 

where Hi{p) ~ niiP^ +p.Q:* with /3 and a* being the usual Dirac matrices. 
The operator P = A^A^ - A^ A^L is a combination of the projection opera- 
tors A^ = [Ea(p)±Ha{p)]/Ea(p), whcrc Ea[p) = ^Mf+p^. It should be 
added that Ea. ll.76| differs from an earlier one by Breit in 1929 '130' by the 
presence of the P factor. In the nonrelativistic limit both equations reduce 
to the same form, but in general the Breit equation has a more limited 
applicability than Eg. 11.761 — as discussed in Ref. |T29j . 

In Refs. |131| , Eq. 11.761 is further simplified by first considering only 
positive-energy solutions i.e. omit the A° A^ terms in P. This results in 
the reduced Salpeter equation 

[E - ^ml + p^ - VmiT^ ]0(p) - j d\Kir^G[ci)^lK\<jy{p + q). 

(1.77) 

Then the formalism is further restricted to the positive energy components 
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to give the semirelativistic spinless-Salpeter equation 



(1.78) 



where V{x) is the Fourier transform of the G(q). Finally the authors of 
Rcfs. |131| go one more step by studying the equal mass case 



where, instead of fixing f3 at 2, they show that /3 > can simulate the 
effect of several particles all of mass m. They interpret this equation as 
"the generalization of the nonrelativistic Schrodinger Hamiltonian towards 
relativistic kinematics" and study algebraically its properties for a variety 
of forms for V{x). The reason why these authors go through explicitly the 
simplification of Ea. lTTSl is to show that Eg. II. 791 — the obvious extension of 
the Schrodinger equation — can indeed be derived in a systematic manner 
and is not just an educated guess. 

1.9.4 Sum rules 

In addition to measuring C^lr) for various values of r, the correlation where 
r is summed over the whole lattice is also obtained. This leads to the 
charge sum rule as discussed in Ref. [SHI- For the charge distribution, the 
outcome is that J^r^^^i^) ~ -^^^ ~ l-3(l)j which is consistent with the 
earlier quenched result jSHl- The fact that X^^ is not unity, as expected 
in the continuum limit, can be qualitatively understood by introducing a 
renormalisation factor of w 1/1.3 w 0.8 into the basic 74 vertex used to 
measure the charge density. Such a factor of this magnitude is reasonable 
as shown in Ref. |132| . It is also reassuring that the X"^ with a ^ (3 are, 
in general, consistent with zero — as expected in the continuum limit. 

The matter sum rule has a somewhat wider spread of values with 0.9(1) 
being a reasonable compromise — a number that is about twice the esti- 
mate of 0.4(1) for the quenched calculation of Ref. [1H|- Perhaps this is 
an indication that, unlike the corresponding matter radial distributions in 
Figure [1.151 the quenched and unquenched results, due to the effect of dis- 
connected processes, can differ even with the present sea-quark masses of 
about that of the strange quark. Certainly differences should appear with 
very light sea-quarks, since then the contribution from the disconnected 
processes, that only enter for the matter distributions, can become signif- 



(1.79) 
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icant. However, we do not have the data to cross check with Refs. |74j . 
which advocate the existence of such differences for the matter sum rule. 

The fact that the matter sum rule is considerably less than that for 
charge can be qualitatively understood by employing data from different 
hopping parameters (k) and using the identity 

X" = ,1.80) 

where ami is the ground state energy and k the hopping parameter — see 
Subsec. Fl. 8. 21 and Ref. 74 . When the mi's correspond to the cases where 
the light quark is of about one and two strange quark masses, Refs. 
and [Sni give X^^ « 0.34(8) and 0.31(6) respectively — consistent with the 
present value of 0.4(1). These values are also consistent with the following 
simple estimate: If the Qq-meson mass (ami) is taken to be simply the sum 
of the quark masses i.e. amq + aruq and k^^ = 8 + 2amg, then Eq. 11.801 
gives 

Xn^d{amQ+am ) ^^^^^ (1.81) 
d{8 + 2amq) ^ ' 

Another reason for expecting < X}j also follows from a potential 
approach using the Dirac equation as in Subsec. 11.9.3.21 This results in 
~ G? + and Xl] ^ Gj - F^. Here G? and Ff are integrals of the 
large and small components of the solution to the Dirac equation. 



1.10 The B - B System as a [{Qq){Qq)] Configuration. 

In Sees. and IL^ the energies and some radial distributions of a single 
heavy-light meson were studied. In this section the interaction between 
two such mesons is extracted using lattice QCD and the outcome is fit- 
ted with an extension of the /-model developed earlier for the [{QQ){QQ)\ 
system in Sees. 11.51 and 11.71 However, it should be added that the study 
of [{Qq){Qq)\ configurations is much less academic than their [{QQ){QQ)] 
counterparts. Many years ago simple multi-quark systems have been pro- 
posed to exist as bound states (133. 134. 135) . Also four quarks forming 
colour singlets or as bound states of two mesons are candidates for parti- 
cles lying close to meson-antimeson thresholds, such as ao(980), /o(980) 
{KK), /o(1500), /2(1500) {uw, pp), /j(1710) iK*K*), ^(4040) {D*D*), 
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T(10580) {B*B*) Systems involving b quarks are particularly in- 

teresting since they should be more easily bound provided the potential is 
attractive, since the repulsive kinetic energy of the quarks is smaller, while 
the attractive two-body potential remains the same. In so-called deuson 
models the long-range potential between two mesons comes from one- 
pion exchange, suggesting that meson-meson systems are significantly less 
bound than meson-antimeson systems. Other models used for realistic 
four-quark systems include string-flip potential models (see Ref. |1.S8| for 
a review), bag models I .'i 1 . and a model-independent approach In 
fact four-quark states with two heavy quarks have been predicted to be sta- 
ble |14()j . Most models give stability for systems where the heavy quarks 
have the b mass, but long range forces might push the required heavy-to- 
light mass ratio down sufficiently so that ccqq states would be bound as 
well. 

1.10.1 Lattice calculation of the [{Qq){Qq)] system 

The interaction between two (Qq) states is depicted in Fig. ll.l9l as a sum of 
two terms — the uncrossed and crossed diagrams. In the latter diagram the 
q hops from one Q to the other. Exploratory studies of two-meson systems 



Fig. 1.19 The interaction between two (Qq) states — the so-called uncrossed and 
crossed diagrams. The solid lines represent Q's and the wavy ones q's. 
NB Compared with earlier figures such as Figs. ll.l'^ and [r.l4l the Euclidean time direction 
is now horizontal. 

have been made for the cross diagram only for SU(3) colour |141| and for 
both diagrams in Refs. [HH, [HSj for SU(2), SU(3) colour respectively. 
This topic is also discussed in Sec. 2.4 of Chapter 2 and Subsec. 4.4.3 of 
Chapter 4. 

In Refs. |77l 144| quenched lattices are used with SU(3) colour and static 
heavy quarks with light quarks of approximately the strange quark mass. 
Also, when the two Q's are at the same point {i.e. R = 0) and by using the 
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SU(3) colour relationship 3 (8) 3 = 3 ® 6, then the two Q's behave like a Q. 
This equivalence implies that the Iq — 1, Sq = 1 state will have the same 
light quark structure as the S;, baryon and the Iq = 0, Sq = state as 
that of the Af, baryon. The other two allowed states at i? = correspond 
to a static sextet source. When the results for i? = are compared with 
the known spectrum of baryons with one heavy quark (A^ and Sfc), good 
agreement is found. Note that this link to baryons at small separation R 
cannot be explored using SU(2) of colour. 

Below, the mass of one quark in each meson is taken to be very heavy — 
the prototype being the _B-meson. The static limit is then the leading term 
in the heavy quark effective theory for a heavy quark of zero velocity and 
there will be corrections of higher orders in 1 /mg , where mq is the heavy 
quark mass. In the limit of a static heavy quark, the heavy quark spin is 
uncoupled since the relevant magnetic moment vanishes which implies that 
the pseudoscalar i3-meson and the vector i?*-meson will be degenerate. 
This is a reasonable approximation since they are split by only 46 MeV 
experimentally, which is less than 1% of the mass of the mesons — see 
Table [T31 Since it is often convenient to treat these two mesonic states as 
if they were degenerate, here they are described collectively as the ;B-meson. 
Because of the insensitivity to the heavy quark spin, it is then appropriate 
to classify these degenerate S-meson states by the light quark spin. The 
system of two heavy-light mesons at spatial separation R will be referred to 
as the BB system. With both heavy-light mesons static, this BB system is 
then described by the two independent spin states of the two light quarks 
in the two mesons i.e. Sq = and 1. Thus there are four possible states 
and it is necessary to classify the interaction in terms of these spin states. 

This situation is very similar to that of the hydrogen molecule in the 
Born-Oppenheimer approximation — with, however, the additional possi- 
bility that the two "electrons" can have different properties. Another simi- 
larity is with the potential between quarks which has a central component 
and also scalar and tensor spin-dependent contributions. 

Each yS- meson will have a light quark flavour assignment. For the BB 
system, it will be appropriate to classify these states according to their sym- 
metry under interchange of the light quark flavours. For identical flavours 
{e.g. ss or uu), we have symmetry under interchange, whereas for non- 
identical flavours {e.g. su or du), we may have either symmetry or antisym- 
metry. For two light quarks, it is convenient to classify the states according 
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to isospin as Iq — 1 (with uu, ud + du and dd) or 1^ — (with ud — du). 
To ensure overall symmetry of the wave function under interchange and, 
assuming symmetry for spatial interchange, the flavour, total light quark 
spin {Sq) and total heavy quark spin {Sh) must be combined to achieve this. 
Thus in the limit of an isotropic spatial wave function i.e. L = 0, there will 
be four different ground state levels of the BB system labelled by (Ig, Sq) 
in the following discussion. These are shown in Table 4.2 of Chapter 4. 

To check for possible finite size effects, the numerical analysis was car- 
ried out on quenched lattices of sizes 12^^ x 24 and 16^ x 24, at /3 = 5.7, 
corresponding to a w 0.18 fm. The bare mass of the light quark was near 
that of the strange mass and light quark propagators were generated using 
the Maximal Variance Reduction method in Subsec. 11.8.81 This enabled 
measurements of the strength of the interaction to be made out to separa- 
tions of i? « 8, which corresponds roughly to 1.4 fm. 

As seen in Figs. 2.5, 2.6 of Chapter 2 and Figs. 4.21 of Chapter 4, attrac- 
tion between two B mesons is found at small values of R for (J^, Sq)—{Ofi) 
and (1,1) and at moderate R (^0.5 fm) for (1,0) and (0,1). For very heavy 
quarks, this will imply binding of the BB molecules with these quantum 
numbers and L = — see Sec. 2.4 of Chapter 2 for more details. 

It is also possible to extract from the lattice data quantities that can be 
identified as tt- and p-exchange between the two S-mesons. This is shown 
in Fig. 11.201 For the 7r-exchange at large R the potential is expected to be 
of the form 

2^^2 -MR, 

ViR)^n.f.L_—. (L82) 

where g/ f-n is the pion coupling to quarks |18f)j . From the lattice studies of 
BB*!: coupling |145| the value of = 0.42(8) is predicted and is the pion 
decay constant (132 MeV). Because the comparison with the lattice results 
is for light quarks, the pion mass is taken to be Ma = 0.53 i.e. a "pion" 
of mass w 580 MeV — a value appropriate for this lattice. Therefore, in 
Fig. 11.201 the solid curve being compared with that data containing one- 
pion-exchange is from Eq. 11.821 and is seen to have the correct features — 
giving support, at these large values of i?, for the deuson model of Ref. |13(ij . 
However, the corresponding p-exchange comparison with the dash-dotted 
line in Fig. II. 201 is less informative, since the normalisation is ad hoc. 

Even though the agreement with one-pion-exchange is seen to be very 
good, the authors of Ref. [77] are quick to point out that their result is 
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0.01 — 



R 

Fig. 1.20 The ratio of the crossed-diagram contributions to the spin-averaged uncrossed 
contribution to the BB correlation. The meson exchange expressions, exp{— M R) / R, are 
compared with these results — using pion exchange with M = 580 MeV (continuous line) 
and rho exchange with M = 890 MeV (dash-dotted line). Note that the pion exchange 
expression is normalised as described in the text, whereas the rho exchange contribution 
has an ad hoc normalisation. Here R is in units of a Si 0.18 fm 1771 . 



"better than should be expected" and give arguments why "This imphes 
that we should not take our estimate of the magnitude of one pion exchange 
as more than a rough guide at the i?- values we are able to measure." It is 
possible that here we are seeing an effect reminiscent of the Chiral Bag 
Model described by Myhrer in Chapter 4 of Ref. [HE!- In such models 
there are pion fields outside some radius « 0.5 fm and only quarks 
inside rf,. Calculations of one-pion-exchange potentials (OPEP) between, 
say, two nucleons then naturally give the usual form of OPEP for r > rt- 
However, for some smaller values of r the interaction seems to be simply a 
continuation of this "usual form of OPEP" i.e. there is a precocious onset 
of OPEP in a region where there are no pions. Perhaps it is this that is 
being seen in Ref. jZZj and Fig. 11.201 It should be added that in some 
studies of interacting clusters, the OPEP is expected to emerge with an 
exponential dependence and not the Yukawa form in Eg. 11.821 This has been 
demonstrated in Ref. j71l , where the authors show that this unconventional 
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form of OPEP is due to the use of the quenched approximation. To avoid 
this problem with the long range part of the interaction, in a later paper 
|187j these authors study the AqAq potential, which does not contain a 
one-pion or one-eta contribution. 

1.10.2 Extension of the f -model to the [{Qq){Qq)] system 

In Sec. 11.51 a model was developed for understanding the lattice energies of 
four static quarks (5(ri)(5(r2)Q(r3)(5(r4) in terms of two-quark potentials 
and is summarised by Eqs. 11.371 11.381 This f-model, although very simple, 
contains the same basic assumptions made in the more elaborate many- 
body models that incorporate kinetic energy e.g. the Resonating Group 
Method described by Oka and Yazaki in Chapter 6 of Ref. It is, 

therefore, reasonable that this simplified f-model can to some extent check 
the validity of these more elaborate counterparts. In Ref. [HD] the model 
in Sec. II. 51 was extended as below to study the interaction between two Qq 
mesons. This resulted in a non-relativistic Schrodinger-like equation 



— for more details see Appendix B. 

So far this model has only been developed for spin independent poten- 
tials, which means that it should only be compared with the spin-averaged 
results of Ref. (J2| shown in Fig. 4.21 in Chapter 4. The outcome from 
Refs. |fi()l 146j is shown in Fig. 11.211 where it is seen by the solid line that 
the use of only two-quark potentials (i.e. in the weak coupling limit / = 1 
OT kf =0) results in a considerable overbinding at i? = 0.18 fm. The dashed 
line shows the effect of using a form factor with kf = 0.6. Admittedly this 
is less convincing than the earlier four static quark case, since the conclu- 
sion depends essentially on only the two data points corresponding to the 
two Q's being 1 and 2 lattice spacings apart. However, the fact that kf 
is consistent with k « 0.5(1) — see Fig. 11.81 — the corresponding param- 
eter needed in Subsec. ri.5.51 for [{QQ){QQ)] configurations in squares or 
near-squares lends support to the general approach of this model. 

In Refs. |148j SB-scattering has been treated in the weak coupling limit 
of fc^ = — a limit that appears to be ruled out by the comparison in 
Fig. 11.211 However, it is possible that this limit can, to some extent, be 
salvaged if the model is extended by including states with excited glue — 
as in Sec. 11.71 and Appendix A. 3 for the [{QQ){QQ)] system. 



|K'(i?) + V'(i?) - £;(4, R)N'{R)\^j; = 



(1.83) 
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Fig. 1.21 Comparison between the spin independent part of the Q^q^ binding en- 
ergies calculated on a lattice [771 (solid circles — quenched approximation with 
a = 0.170 fm) 11471 (solid squares - with dynamical fermions and a = 0.142 fm). The 
crosses, with the solid line to guide the eye, show the model in the weak coupling limit 
{kf = 0). The dotted line shows the result with kf =0.6 and m5=400 MeV. The dy- 
namical fermion data was not used in any fit. They are simply included to show that 
it is qualitatively consistent with the quenched data but with considerably smaller error 
bars. 

1.11 The B - B System as a [{Qq){Qq)] Configuration 

The B-factories discussed in Subsec. 11.3.11 are not able to study directly 
BB reactions. However, the related BB system is accessible as a final state 
in the decay of the T(4S', 10580 MeV), whose main branching (> 96%) is 
into BB. 

At first sight it may be thought that the [{Qq){Qq)] and [(Qg)((5^)] 
configurations have similar properties. However, this is not so, since now 
the q and q can annihilate each other. This means that there is a coupling 
between {QQ) and [{Qq){Qq)] states i.e. the string (fiux) connecting the Q 
and Q in the {QQ) state can break into two mesons {Qq) and {Qq). This 
becomes very clear if the mechanism for the interaction between two {Qq) 
states, shown in the Fig. 11.191 is compared with the corresponding ones 
for the present \{Qq){Qq)\ case on the first row of Fig. 11.221 — only the 
uncrossed contributions look similar. The diagram b) represents the two 
step process involving the annihilation and creation of a qq pair 
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[iQq)iQq)]^{QQ)^[iQq)iQ<l)]- 

This breaking of a long flux tube between two static quarks into a quark- 
antiquark pair is one of the most fundamental phenomena in QCD. Because 
of its highly non-perturbative nature it has defied analytical calculation, 
while its large scale, e.g. when compared to the sizes of composite particles 
in the theory, has caused difficulties in standard nonpcrturbativc methods. 
Thus string breaking has remained a widely publicized feature of the strong 
interaction that has never, apart from rough models, been reproduced from 
the theory. 

String breaking can occur in hadronic decays of QQ mesons and is 
especially relevant when this meson is lying close to a meson-antimeson 
{Qq){Qq) threshold. Its effect should be seen most directly by measuring 
the QQ potential, since the onset of string breaking would change the form 
of the standard QQ potential in Eg. II. 61 so that the confining term — br for 
all r — would become br only for r < Tc and have the constant value hvc for 
r > Tc- Unfortunately, this direct approach of trying to see the flattening in 
the static QQ potential at large separation has only had limited success. An 
examDle |149) of the usual outcome is seen in Fig. 11.231 There the potential 
V{QQ) continues to rise linearly way past the value of r/rp « 2.4, where 
the breaking into two mesons should occur — denoted by the horizontal 
dotted lines. There is clearly no sign of the expected flattening. The 



a) 




b) 




d) 



Fig. 1.22 The [{Qq){Qg)] interaction and string breaking. 

The a) uncrossed and b) crossed [{Qq){Qq)] — > [iQq){Qq)] contributions to the inter- 
action; c) The (QQ) (QQ) Wilson loop and d) the [{Qq){Qq}] (QQ) off-diagonal 
correlation. 
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Fig. 1.23 The static potential V{QQ) obtained from simply the Wilson loop of Fig. ll.l| 
N.B. There is no sign of the expected flattening at r/rg fs: 2.4, where it becomes ener- 
getically favourable to create two mesons [l49l . 

failure of this Wilson loop method seems to be mainly due to the poor 
overlap of the operator (s) with the [{Qq){Qq)] state [1501 151| . In three- 
dimensional SU(2) with staggered fermions an improved action approach 
has been claimed to be successful with just Wilson loops |152| . However, 
in full QCD with fermions, effective operators for [iQq){Qq)] systems are 
hard to implement; part of the problem is the exhausting computational 
effort required to get sufficient statistics for light quark propagators with 
conventional techniques for fermion matrix inversion. In Ref. 153 this 
problem was to a large extent overcome by applying the Maximal Variance 
Reduction method of Subsec. [1.8.31 using SU(3) on a 16"^ x 24 lattice with 
the Wilson gauge action plus the Sheikholeslami-Wohlert quark action and 
P = 5.2 {i.e. a « 0.14 fm) with two degenerate flavours of both valence- 
and sea-quarks. 

Since Fig. 11.221 b) is a two-step process it is natural to consider the 
process as a coupled channels problem between the [{Qq)iQq)] and (QQ) 
configurations. This leads to the two new processes in Fig. ll.22l d^ and c) — 
the off-diagonal term [{Qq){Qq)]< — *{QQ) and the corresponding diagonal 
term {QQ)< — *{QQ), which is nothing more than the Wilson loop discussed 
in Subsec. UXTJ 

In Fig. 11.241 the results from a calculation using just the most fuzzed 
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basis states for both (QQ) and [{Qq){Qq)] (a 2 x 2 matrix) are shown. Un- 
fortunately, the statistics are not sufficient to give accurate plateaux for the 
energies. However, the authors of Ref. 1!t^ are able to extract a quantity 
they call the mixing matrix element x, which can be interpreted as an in- 
direct measure of possible string breaking. They find that x — 48(6) MeV. 
At Tc one would expect the ground and excited state energies to be sepa- 
rated by 2x. However, in Fig. ll.24l a larger separation is observed, which is 
presumably again due to insufficient statistics. In Ref. |I54) . utilizing the 



3.0 



A- - AGround 

- Vlst excited 
& - O Wilson 




Fig. 1.24 Ground and excited state from a variational calculation including QQ and 
QqQq operators 11531 . The highest fuzzed basis state for both is used here. The ground 
state of the Wilson loop and 2m,Qq arc also shown. 



SU(3) colour relationship 3 ® 3 = 1 © 8, the QQ at i? = behave as a sin- 
glet(vacuum), so that the remanent qq component with {Sqq,Iqq) = (0, 1) 
can be compared successfully to a pion with non-zero momentum. 

In the above, the failure of the QQ correlations alone to give a flattening 
of V{QQ) was considered a negative feature. However, in the Conclusion 
it will be seen that this "failure" could be useful for constructing models. 

So far no one has attempted to understand this BB data with the 
extended /-model of Subsec. ri.lU.2l for BB states. Such an extension would 
also have to incorporate qq creation and annihilation using some model such 
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as the so-called Quark Pair Creation or ^Pq model |155| . When this ^Pq 
model was combined with the harmonic oscillator flux tube model of Isgur 
and Paton (31, it proved successful for describing flux tube breaking and 
formation |15f)j . 



1.12 Conclusions and the Future 

In this chapter there has been an attempt to bring together two distinct 
lines of research: 

(1) Lattice QCD was applied to the multiquark systems [{QQ){QQ)], 
iQq), [{Qq)iQq)] and [iQq)iQq)] in Sees. imi IlO OniandlLTTl 

(2) Effective Potential Theories (EPTs) — based on interquark 
potentials with four-quark form factors included — were developed 
in Secs. ll.5tJ1.7l and ll.l().2l to interpret phenomenologically this lat- 
tice data. 

However, most of the comparison between these two lines has been devoted 
to constructing EPTs that give some phenomenological interpretation of the 
Lattice QCD energies. For the [{QQ){QQ)] system the latter concentrated 
on the energies of the six geometries in Fig. 11.101 and this amounted in 
all to about 100 pieces of data. In spite of four static quarks being a 
very unphysical system, the Lattice QCD <-> EPT comparison showed that 
a 4-quark potential seemed to be needed and that simply using a sum 
of 2-quark potentials failed by generating far too much binding for large 
interquark distances. On the other hand, for small interquark distances the 
use of only 2-quark potentials was sufficient — the so-called weak coupling 
limit. However, in Sect. ll.7l it was found that, if the model was extended to 
explicitly include gluonic excited states, then much of the attraction came 
from these excited states and that the need for a 4-quark potential was 
greatly reduced. This observation for the interaction between four static 
quarks suggests that the usual approach of simply using 2-quark potentials 
in multiquark problems could, to some extent, be justified provided gluonic 
excited states are explicitly included. 

In the more physical case of the [{Qq){Qq)] system the EPTs are com- 
plicated by the presence of the light quark kinetic energy and mass. Even 
so, the Lattice QCD ^ EPT comparison in Fig. 11.211 still shows the same 
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effect that simply using a sum of 2-quark potentials fails, if the basis states 
only contain the gluon field in its ground state. 

This apparent need for a 4-quark potential can be viewed as a form 
factor — a familiar and successful technique in, for example, parametriza- 
tions of the NN-potential |2ZI- However, there the form factor is needed to 
regularize the potential at small values of r, where the potential can have 
l/r" singularities. It can, therefore, be simply incorporated as a short- 
ranged vertex correction modelled from meson fields. In contrast, the form 
factor needed in the present model is introduced to eliminate problems — 
essentially the van der Waals effect — at large values of r. This form factor 
is, therefore, modelling a long-ranged effect which could be excited gluon 
states, which become more effective the larger the system. 

Clearly for a better "Lattice QCD ^ EPT" comparison more data is 
needed from a given quark system. This suggests that comparisons should 
be made using other quantities in addition to the few available lattice en- 
ergies. The most obvious candidates are radial distributions of the light 
quark, since — being a function of distance — these introduce many new 
pieces of data in a way that is more systematic than the earlier choice of 
simply the six "convenient" geometries of Fig. 11.101 In addition, there can 
be several types of radial distribution: 

• The charge (vector) density where the probe in Eg. 11.651 is Q — J4. 

• The matter (scalar) density where the probe is = 1 . 

• The pseudovector density with the probe (7^75). This is needed 
for the B*Btt coupling jl45| and was exploited in Subsec. Fl.lO.ll 

• In addition to operators that probe the radial distributions of the 
light quark(s), it is also possible to map out the structure of the 
properties of the underlying gluon field. This is achieved by us- 
ing different orientations of the elementary plaquette defined in 
Eq. 11.161 Using purely spatial plaquettes the radial distributions 
of the various components of the colour magnetic field can be 
extracted e.g. Un^ B^, whereas those plaquettes with a eu- 
clidean time superfix give the spatial distributions of the vari- 
ous components of the colour electric field e.g. [/□' E^. In 
Ref. II57 158 159_ these distributions were calculated for two- and 
four-static SU(2) quark systems, where the latter were restricted to 
the corners of squares with sides upto 8 lattice spacings (i.e.« 1 fm). 
The lattice was 20^ x 32 with (3 = 2.4 (a « 0.12 fm). There 
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Fig. 1.25 Comparison of the neutron-proton radial correlations in the deuteron, ^He 
and *He HM . 



is a wealth of information in these flux profiles and they present 
a formidable challenge to models that attempt to describe them 
|54l 55| . However, the few comparisons that have so far been made 
are very encouraging |157| . 

In principle, all of these distributions can be measured for the light quarks 
in the 4-quark systems [{Qq){Qq)] and [{Qq){Qq)] discussed in Sees. ILlHl 
and 11.111 However, this is probably too ambitious at the present time 
for meaningful Lattice QCD <-> EPT comparisons. First, the radial dis- 
tributions in the 2-quark (Qq) system should be better understood — the 
topic of Sec. 11.91 It is possible that the radial distributions of the q in 
the [{Qq){Qq)] and [{Qq){Qq)] systems are related to that in the simpler 
(Qq) case, since for the corresponding situation in light nuclei the neutron- 
proton radial correlations in both ^He and '^He are very similar to that 
in the deuteron — see Fig. 11.251 This shows that, although the 3- and 
4-nucleon calculations needed to extract such correlations are very compli- 
cated, they result in some simplicities. Possibly comparisons of correlations 
within few quark systems could lead to similar simplifications and so en- 
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hance our understanding of such systems. The work in Rcf. "EP for the 
QQq system is a step in this direction. -I- 

In Sec. ll.9l the charge and matter distributions of the q in the Qq system 
( w i?-meson) were measured and attempts made to fit these distributions 
with simple functions in Subsec . fl .9 . 3 . H and also by using the Dirac equation 
in Subsec. [1.9.3.21 The latter can be viewed in two ways: 

(1) As simply an alternative form of parametrization of the two distri- 
butions with no physical interpretation of the parameters needed. 

(2) As the construction of an EPT, in which the parameters do have 
a physical interpretation — albeit phenomenological — on which 
extensions to multiquark systems can be based. 

This leads us to the main problem of how to set up an EPT in order to 
understand multiquark systems. In the past the key word has been "mimic- 
ing" , in which models based on potentials have been mimicing the successful 
models of multi-nuclcon systems. Unfortunately, since such models aim di- 
rectly at the experimental data, they may be attempting to describe some 
mechanisms that are outside the scope of the model. This I have called 
the Nuclear-Physics-Inspired- Approach (NPIA) and it would correspond 
to Option 1 in the comparison 



1 2 

NPIA < — > Experimental Data versus Lattice data < — > QCDIA 

However, I believe it is more reasonable to try to create QCD-Inspired- 
Approaches (QCDIA) as in Option 2. Here the models attempt to mimic 
directly details of the lattice data that have been obtained under "controlled 
conditions" and so, hopefully, do not contain undesirable processes not 
included in the model. In this way the data have a better chance of deciding 
the form of the model. In contrast, with the NPIA the form of the model 
tends to be decided beforehand with the experimental data only leading 
to a tuning of the parameters. An example of this was demonstrated by 
Fig. 11.231 There it was seen that using only {QQ) correlations did not 
lead to a flattening of the interquark potential V{QQ) — even though 
such a flattening should arise with the onset of {Qq) and {Qq) mesons 
being created. Therefore, for an understanding of V{QQ) this defect was 

■fSome recent lattice calculations on doubly-charmed baryons have concentrated on their 
masses lltill . 
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considered to be a negative feature, since this implementation of LQCD 
clearly did not agree with experiment. However, for model building this 
could be envisaged as a positive feature, since LQCD is now generating 
exact data from the QCD Lagrangian that can be interpreted by a model 
that is simply a linearly rising potential between a single Q and a single Q. 
This intermediate stage model could, hopefully, be more easily extended to 
the real life situation, when the more complicated correlations in Fig. 11.221 
are treated by LQCD. This intermediate stage could be considered as yet 
another example of the "unphysical worlds" discussed in Sec. 11.2.21 in the 
context of models with different numbers of colours, spatial dimensions, 
quark masses etc.. This recalls my earlier work in ^^0 |lf)2l Ifii-ij . There the 
"real" ^^O needed to be described by "4-particle 2-hole" states in addition 
to the usual "2-particle 0-hole" states. This is analogous to the matrix of 
correlations in Fig. 11.221 It would have been very useful if there had been 
an "experimental phase" of ^^O that only needed 2-particle 0-hole states 
for its description. In that case the model for "real" ^^O would have had 
less freedom. 

So far the QCDIA has only been attempted for the [{QqjiQq)] con- 
figuration in Sec. 11.101 There the (Qq) lattice data in Fig. 11.181 was first 
fitted with a Schrodinger equation to give an effective light quark mass 
iTT-q, offoctivo ~ 400 McV. Usiug a variational method, this Schrodinger ap- 
proach was then easily extended to the Q^q^ system using the same in- 
terquark potential and m^^ offoctivc determined earlier from the Qq data. In 
principle, this could be extended to any system that can be described in 
terms of interacting quark clusters. 

The problem with the above QCDIA is that it usually results in an 
rriq^ effective ^ 1 GcV Suggesting the need for a more relativistic approach. 
This can be attempted at different levels: 

(1) At the one extreme we can use directly the Dirac equation to de- 
scribe, say, the Qq system as in Subsec. 11.9.3.21 This had some 
success but we are then confronted with the problem of extending 
the comparison to multi-quark systems — a second step (2) that 
is not directly possible for systems containing more than one light 
quark i.e. 

_ 1 2 ^'^ 

Lattice data for Qq — > 1-quark Dirac Eq. — ^ Multi-quark case 



(2) Use some Semirelativistic Schrodinger-like Equation (SRSE) as de- 
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scribed in Subsec. ri.9.3.41 This is most easily formulated in momen- 
tum space with the kinetic energy Enr = niq^e+P^ /'^niq^e being re- 



chosen to fit some piece(s) of experimental or lattice data. How- 
ever, care must be taken to treat the potential terms to the same 
semirelativistic degree by inserting appropriate factors of mq^e/Ep 
as in the Blankenbecler-Sugar equation 35 . The second step (2) 
to multiquark systems is then probably possible i.e. 

_ 1 2 

Lattice data for Qq — > SRSE — > Multi-quark case 

(3) Thirdly, a compromise model emerges. First the Qq lattice data 
is fitted by a Dirac equation. This equation then generates other 
observables that are interpreted in terms of a SRSE, which can be 
extended to multi-quark systems. This is similar to the philosophy 
of Bhaduri and Brack ,,30) . who show how the Schrodinger equa- 
tion with a quark effective mass of mq^ effective ~ 500 MeV is able 
to explain some of the results — energies and magnetic moments 
— of a Dirac equation for a zero mass quark i.e. they make the 
comparison 

_ 1 2 3 

Lattice data for Qq — > Dirac Eq. — > SRSE — > Multi-quark 

So what are the "Bridges from Lattice QCD to Nuclear Physics" as adver- 
tized in the title of this chapter? So far, the only clear examples involve 
mainly static (or heavy) quarks (Q) as in the extraction of the string energy 
and the lattice spacing from Vqq with a non-relativistic Schrodinger equa- 
tion fSec. ll.2|l and the energies of the various Q"* geometries in Sees. I1.3H1.7I 
The introduction of light quarks as in the (Qq), [{Qq){Qq)] and [{Qq){Qq)] 
systems in Sees. II .81 - 11 .lll leads to major complications and the only partial 
success is the Schrodinger description of the (Qq) + [iQq)iQq)] energies in 
Sec. ll.lOl The real test of whether any "Bridge" exists is only now becom- 
ing possible with the advent of the lattice data for radial distributions, since 
such distributions are also at the centre of much of Nuclear Physics. Of 
course, the final outcome could be that there is no useful bridge for treating 
multiquark systems in the way we treat multinucleon systems. This would 
mean that Lattice QCD would only ever be able to directly address prob- 
lems involving a few quarks — perhaps upto the six quarks needed for the 
nucleon-nucleon interaction — but not show a general way for how to deal 
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with multiquark systems. This would result in the two worlds of QCD and 
Nuclear Physics having little direct connection with each other. 

However, before such a pessimistic view is adopted, we should remember 
that Nuclear Physics earlier had another two-world structure that lasted 
for many years. Until the 1960's there were essentially two models for 
the nucleus — the collective liquid-drop-like model that rarely mentioned 
the nucleon-nucleon potential and, in contrast, the shell model based on 
this potential. But with the advent of Brueckner theory it was shown 
how effective interactions in many-nucleon systems could be constructed 
from the basic nucleon-nucleon potential. This was followed up by the 
generation of collective effects as interacting particle- hole states |164| . In 
this way a bridge was made between the basic nucleon-nucleon potential 
and collectivity — but it took many years. 
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Appendix A: Extensions of the /-Model from 2 X 2 to 6 X 6 

A.l: The 3x3 extension of the unmodified two-body approach 
of Subsec. \1.5.1\ 

Because the colour group considered so far is SU(2), there is no distinction 
between the group properties of quarks {Q) and antiquarks (Q) . Four such 
quarks can then be partitioned as pairs in three different ways 

A^{QiQ3){Q2Qi), B ^ {QiQi){Q2Qz) and C = {QiQ2){QiQi), 

(A.l) 

where each (QiQj) is a colour singlet. However, these three basis states are 
not orthogonal to each other. Also, remembering the fact that the quarks 
are indeed fcrmions gives, in the weak coupling limit, the condition in the 
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Appendix of Ref. 

\A) + \B) + \C) = 0. (A.2) 

Since = {B\B) ~ {C\C) — 1, we get — in this limit — the equalities 

{A\B) = {B\C) = {A\C) = -1/2. 

In SU(3) the partitioning problem is different since in that case the three 
partitions are 

A = {QiQ3)iQ2Q4), B = {Q^Q^){Q2Q^) and C = [{Q^Q2)''{Q3Qa)\ 

(A.3) 

where state C is expressed in terms of either colour antitriplet {d = 3) or 
sextet {d = 6) states and so cannot appear asymptotically as two clusters. 

If all three basis states in SU(2) are included, then the matrix below is 
singular for the obvious reason that |j4) + + |C) =0 i.e. 

/ 1 -1/2 -1/2 \ / 1 1/2 1/2 \ 

detN = det -1/2 1 -1/2 = dot 1/2 1 -1/2 = 0. 
V -1/2 -1/2 I I V 1/2 -1/2 1 / 

(A.4) 

Earlier this was interpreted to mean that it was unnecessary to include all 
three states and so the symmetry was broken by keeping the two states 
with the lowest energy, let us say , A and B. A similar effect also occurred 
in the lattice simulations. There it was found that the energy of the lowest 
state was always the same in both a 2 x 2 and 3x3 description, providing 
A ox B had the lowest energy of the three possible partitions. In addition 
the energy of the second state was, in most cases, more or less the same - 
the largest difference occurring with the tetrahedral geometry. 

A.2: The 3x3 extension of the f -model of Subsec. \1.5.S!\ 

The /-model of Subsec. fl. 5. 21 by incorporating multiquark interactions, had 
the good feature that, when fitting the data {Eq,Ei) for a given square, 
only a single / was necessary to get a reasonable fit to both energies — 
see Fig. 11.71 Of course, / was dependent on the size of the square, but a 
reasonable parametrization was 

/(/a) = exp(-6sfc/S') (Version la), (A. 5) 

where S is the area of the square and kf « 0.5. Earlier in Ea. ll.43l a slightly 
different notation was used. The original hope was that, with kf determined 
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from the squares and nearby rectangles, the model would automatically also 
fit other geometries with S being the "appropriate" area contained by the 
four quarks. When the four quarks lie in a plane, the definition of S is 
clear. However, in non-planar cases the situation is more complicated. One 
possibility is to simply take S to be the average of the sum of the four 
triangular areas defined by the positions of the four quarks i.e. the faces of 
the tetrahedon. For example, in the notation of Eq. lA.ll the appropriate 
area S{AB) for / is 



where S(ijk) is the area of the triangle with corners at i,j and k. For 
planar geometries this simply reduces to the expected area, but for non- 
planar cases this is only an approximation to S{AB) - a more correct area 
being one that is not necessarily a combination of planar areas but of curved 
surfaces with minimum areas. These possibilities are discussed in Ref . I165j . 
It would be feasible to incorporate this refinement here, since only a few 
(~ 50) such areas are needed for the geometries in Fig. 11.101 However, 
for a general situation, in which the positions of the quarks are integrated 
over, it would become impractical to use the exact value of S{AB), since 
the expression for the minimum area itself involves a double integration. 
In contrast, the area used in Eq. IA.6I is an algebraic expression and is, 
therefore, more readily evaluated for any geometry. The above model will 
be referred to as Version la. 

This model la has only the one free parameter kf in /(la) of Eg. IA.5I 
Another possibility with additional parameters /o, kp is 



where P is the perimeter bounding S. This form has been used in Refs. |lf)5j . 
However, as shown in Ref. |49j . this reduces in the continuum limit to the 
same as Version la — the differences at /? = 2.4 being mainly due to lattice 
artefacts. 

Unfortunately, in the 2x2 version both of these models have the feature 
that, for regular tetrahedra, they are unable to reproduce a degenerate 
ground state with a non-zero energy, since the two eigenvalues are 



S{AB) = 0.5[S'(431) + 5(432) + S'(123) + S'(124)], 



(A.6) 




(Version lb) , 



(A.7) 



Eo 



//2 



[Vcc - Vaa] and Ei 



.//2 



[Vcc ~ Vaa], (A.8) 



1 + //2 



1-//2 
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where, in the notation of Fig. II.IUI Vaa — ^13 + W24 and Vcc = + 1^23 
and so for regular tetrahedra Vcc = Vaa, giving Eq = Ei = 0. 

Prior to the work on tetrahedra the geometries considered had, at most, 
two of the three possible partitions being degenerate in energy {e.g. for 
squares) — see Subsec. 11.5.21 In these cases, it is found that the lattice 
energies Eq and Ei are essentially the same for the three-basis-state calcu- 
lation {A + B + C) and those two-basis-state calculations {A + B, A + C 
and effectively B + C) which involve the basis state with the lowest un- 
perturbed energy. This is one of the reasons why the 2x2 version of the 
/-model in Eq. 11.381 was quite successful for a qualitative understanding 
of these cases. However, for tetrahedra and the neighbouring geometries 
calculated in Subsec. 11.61 it now seems plausible to extend the /-model to 
the corresponding 3x3 version of 

[V(/)-A,(/)N(/)]vI/, =0, (A.9) 

in which 

/ 1 //2 f/2 \ 
N(/) - //2 1 -f"/2 and (A.IO) 
V f'/2 -r/2 1 J 



Vl3 + V24 fVAB I'Vac 

V(./)=| fVsA v^i + v23 -f'VBc \, (A.ll) 

f'VcA -f'VcB V12 + V34 

where the negative sign in the BC matrix elements is of the same origin as 
the one in Eas. 11.511 and lA.4l 

This apparently leads to the need for two more factors /', /" defined by 



{A\C) = -r/2 and {B\C) = -/72. (A.12) 

However, with the parametrizations of / as in Eqs. IA.5I or I A. 71 and the 
definition of S as in Eq. \EM it is seen that f ^ f" = f, since S is 
simply proportional to the area of the faces of the tetrahedron defined by 
the four quark positions and is independent of the state combination used. 
Therefore, the 3x3 model has for all 4-quark geometries a form where the 
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N and V matrices are 

/ 1 //2 f/2 \ 

N(/) =[ f/2 1 -f/2 and (A.13) 

V f/2 ~f/2 1 / 

(VAA fVAB fVAC \ 

fVBA VBB -fVBC ■ (A. 14) 

fVcA -fVcB VCC J 

This extension from 2 x 2 to 3 x 3 has the good feature that aU three 
basis states are now treated on an equal footing. This is convenient when 
considering some general four-quark geometry, since it is then not necessary 
to choose some favoured 2x2 basis, which could well change as the geometry 
develops from one form to another. In the weak coupling limit [i.e. f —> 1) 
the 3x3 matrix in Eq. IA.13I becomes the singular matrix in Eq. I A. 41 
However, in this limit, each of the 2x2 matrices corresponding to the three 
possible partitions A+B, A+C and B+C gives the same results. Away 
from weak coupling the 3x3 matrix is no longer singular, but now the 
three possible 2x2 partitions do not necessarily give the same results. In 
addition to this general problem as / — > 1, there are also the following more 
specific unpleasant features: 

(1) For regular tetrahedra all eigenvalues are zero as in the 2x2 models. 
The reason for this is clear. There is only one energy scale in the 
model, since all the Vij are the same. Therefore, there can not be 
any excitations. 

(2) For a linear geometry, since the "appropriate" area as defined by 
Eg . I A . 61 vanishes . we get f = 1 i.e. we are back to the weak coupling 
limit and a singular matrix. 

(3) For squares the model gives Ei — —Eq, whereas the predictions of 
the 2x2 version in Eos. lA.Sl seem to be nearer the lattice data. 

(4) The differences between using the various combinations of the three 
partitions are often considerably larger than in the corresponding 
lattice calculation. 

The most glaring problem is the fact that the 3x3 model for the tetrahe- 
dron gives three degenerate states with zero energy, because there is only 
one scale in the model. It is, therefore, necessary to introduce a second 
energy scale. However, any improvements in the model have very limited 
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choices, since there are only two different matrix elements involved — the 
diagonal ones all equal to —E and the ofF-diagonal ones all equal to ±0.5 fE. 
Therefore, the most general modifications are to change the diagonal ma- 
trix elements to di — E and the off-diagonal ones to ±0.5/((i2 — E). This 
results in the eigenvalues 



At first sight it may appear that there is sufficient information to now 
extract the new parameters di,2, since / can be estimated using the pa- 
rameters (assumed to be universal) from other geometries — thus leaving 
two equations for i?o,i and E2 and the two unknowns rfi,2- However, as 
said before, this is too much to demand from the /-model, since in the 
lattice calculation the third basis state in the complete A-l-B-l-C basis gen- 
erally plays a minor role in determining the values of -Bo.i and, therefore, 
the third eigenvalue is presumably dominated by an excitation of the gluon 
field. Furthermore, it is of interest to see that a similar feature now arises 
with E2 in Eq. IA.15I since this third state is removed in the weak coupling 
limit i.e. — > 00 as / — > 1. However, in its present form, the /-model is 
only expressed in terms of the lowest energy gluon configurations, since the 
gluon field is not explicitly in its formulation, but only appears implicitly 
in the form of the two-quark potentials and the /-factors. But already at 
this stage we see from the behaviour of E2 that the effect of excited gluon 
states seem to be playing a role — the topic of the next subsection when the 
model is further extended to a 6 x 6 version. In view of this, no quantitative 
attempt should be made to identify the second excited state emerging from 
the lattice calculation as E2 in this 3x3 version of the /-model. 

In Ref. |1(J2) it was shown that the two-state model of Eas. n'.38l with the 
overlap factor / = 1 agreed with perturbation theory up to fourth order in 
the quark-gluon coupling [i.e. to 0(q;^)] and gave Eq i—Q for tetrahedra. 
Therefore, the non-zero lattice results for small tetrahedra must be of O(a^) 
at least. Another aspect of this special situation for tetrahedra is also seen 
— when extracting or interpreting the value of Ei — by the need for the 
third basis state both in the lattice calculation and in the /-model, since 
in comparison with Eq. IA.15I the two-basis-state version gives 



Eq — El 



di + 0.5/d2 
1 + 0.5/ 



and E2 



di - fd2 
1-/ 



(A.15) 



di + 0.5/d2 
1 + 0.5/ 



and El 



di - 0.5/d2 
1-0.5/ 



(A.16) 
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i.e. both the two- and three- basis-state models have the same ground state, 
but the former does not show the Eg = Ei degeneracy. 

The expressions in Eq. IA.15I are not particularly useful unless there 
is a model for the parameters (ii.2- However, since it is not the purpose 
at this stage to make a comprehensive study of models covering all the 
4-quark geometries considered in earlier works |52l 8()l 871 88| . only a few 
general remarks will be made here for the tetrahedron geometry. Models 
for the di^2 need extensions of the potential in Eq. 11.311 so that for the 
tetrahedron two energy scales arise. Here several ways of achieving this 
goal are suggested: 

The effect of an isoscalar two-quark potential. 

As discussed in Ref. 87 , an isoscalar potential Wij can be introduced into 
Vij — still ensuring Vij = Vij for a colour singlet two-quark system — by 
extending the form in Eq. 11.61 to 



In this case, di = d2 = '^w, since all of the Wij are now equal to w and 
results in Eq = Ei = E2 = 4w. Therefore, here w takes on values that 
range from -0.0035 to -0.0070 i.e. they have values much smaller than the 
corresponding Vij — v in Eq. 11.311 A similar feature was found in Ref. |87| , 
when the form in Eq. IA.17I was introduced to improve the model fit for 
squares and rectangles. However, as shown in Ref. jl()2j . in perturbation 
theory all terms of 0{a^) are included in the two state model of Eg. 11.381 
with / = 1. Therefore, in the weak coupling limit Wij must be of 0{a^) at 
least. 

The effect of a three- or four-body potential. 

The / factor is itself a four-body operator. However, it is conceivable 
that additional multiquark effects arise. Some perturbative possibilities 
are discussed in Ref. |102| . There it is shown that all three-quark terms 
arising from three gluon vertices always vanish, but that the four-gluon 
vertex can contribute to 2-, 3- and 4-quark terms at O(a^). However, in 
the tetrahedral case {r — d), cancellations result in this particular 4-quark 
term also vanishing. 



V^j =--T,-Tj {v^j - Wy ) + Wy. 



(A.17) 



The effect of non-interacting three gluon exchange processes. 

These are also discussed qualitatively in Ref. |102j and contribute at O(a^) 



Februa ry 8, 2008 



23:53 WorldScientific/ws-b8-5x6-0 



chap terSrv 



Appendix A 



103 



to 2-, 3- and 4-quark potentials. 

The effect of two quark potentials where the gluon field is excited. 

The first excited state of the two-quark potential V{r) is approx- 

imately given by V*{r) « V{r) + ir/r — see for example Refs. |90l 91| . 
Therefore, if a fom'th state, based on such an excited state, is introduced 
into the model, it will give attraction in the ground state, since it is higher 
in energy than the three degenerate basis states so far considered. Further- 
more, as the size of the tetrahedron increases this fourth state will approach 
the other three states, so that the attraction felt in the ground state will 
increase — a trend seen in the tetrahedron results for i^o.i in Fig. 11.111 
This possibility will be discussed more in the 6x6 extension below. 

The above isoscalar potential option now offers a reason for Eq — Ei ^ 0. 
But, unfortunately, E2 is still equal to i?o,i since di = d2- However, there 
is no reason to expect any three or four body forces to also be purely 
isoscalars. In this case, their contributions to di and d2 could be different 
and through the presence of the (1 — /) factor in Ea. lA.lSl anv estimates of 
E2 could be very model dependent. 

A. 3: The 6x6 extension of the f -model of Subsec. \1.5.^ 

The above models both have trouble in describing regular tetrahedra. In 
Refs. |1()81 109) an attempt is made to overcome this problem. An interest- 
ing feature of the regular tetrahedron data is that the lowest state becomes 
more bound as the tetrahedron increases in size with the magnitude of 
Eq increasing from -0.0202(8) to -0.028(3) as the d^ cube containing the 
tetrahedron increases from d = 2 to d = 4. This is opposite to what hap- 
pens with squares, where the magnitude of Eq decreases from —0.0572(4) 
to —0.047(3) as d increases from 2 to 5. This indicates that there could be 
coupling to some higher state(s) that becomes more effective as the size in- 
creases and suggests that these higher states contain gluon excitation with 
respect to the A, B, C configurations. Therefore, the 3x3 model in the 
previous subsection is further extended to a 6 x 6 model by adding three 
more states A*,B*,C*, where in analogy with Eq. lA.ll 



A* ^ {QiQ3)eAQ2Q4)e,, etc.. 



(A.18) 
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Here {Q1Q3) denotes a state where the gluon field is excited to the lowest 
state with the symmetry of the representation of the lattice symmetry 
group D^h- Because it is an odd parity excitation, A* ,B* , C* must contain 
two such states in order to have the same parity as A, B, C. The excitation 
energy of an £"„ state over its ground state (Aig) counterpart is ~ ir/r 
for two quarks a distance r apart. As r increases this excitation energy 
decreases making the effect of the A* ,B* ,C* states more important, leading 
to the effect mentioned above. Here we have assumed that these states arise 
from a combination of excited states with However, it is possible that 
they involve other excitations, e.g. 

A* = iQiQ3)A'JQ2Q4) etc., (A.19) 

where the A[g state is a gluonic excitation with the same quantum numbers 
as the ground state (Aig). For this case the following formalism would be 
essentially the same. Another possibility, which is not considered here, is 
that the relevant excitations are flux configurations where all four quarks, 
instead of two, are involved in forming a colour singlet. In the strong cou- 
pling approximation such states would reduce to two-body singlets due to 
Casimir scaling of the string tensions, namely, the string tension for a higher 
representation would be more than double the value of the fundamental 
string tension, thus preventing junctions of two strings in the fundamental 
and one in the higher representation. This would happen both in SU(2) 
and SU(3), the only exception being the unexcited C state in SU(3), which 
would involve an antitriplet string — see Eq. I A. 31 

For the regular tetrahedral case, in addition to / = /' = /", there are 
now several new matrix elements that need to be discussed for N(f) and 
V(f) (Eqs.EHSmi): 

a) With the inclusion of the A*,B*,C* states and the antisymmetry 
condition \A*) + \B*) + \C*) = analogous to Eg. IA.2I there are now two 
more gluon overlap functions /"'"^ defined as 

(A*\B*) = (A*\C*) - {B*\C*) = -r/2 and 

{A*\B) = {A*\C) = .. etc. .. = -/72. (A.20) 

Here it is assumed that f'^ '^ are both dependent on S as defined in Ea. lA.6l 
Since involves only the excited states, it is reasonable to expect it has a 
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form similar to / in Eg. lA.Sl z.e. 

r = eM-bskcS). (A.21) 

b) By orthogonality {A\A*) = {B\B*) = {C\C*) = 

c) In the weak coupling limit, from the \A*) + \B*) + \C*) = condition, 

we expect {A\B*) = {B\C*) = = at small distances. To take this into 

account /° is parametrized as 

r = (/r + bsf.^S) exp{-b,kaS). (A.22) 

If all three parameters fi,f2,ka are varied, it is found that /f is al- 
ways consistent with zero — as expected from the above condition that 
{A\B*) ~ .... = 0. Therefore, usually is fixed at zero. 

d) For the potential matrix V(f) the diagonal matrix elements, after 
the lowest energy Vdd amongst the basis states is removed, are 



{A*\V - Vdd\A*) = vis + y*24 - Vdd, etc., 



where Vdd = min[VAA = vis + ^24, Vbb = vu + V23, Vcc = i'i2 + W34] and 
v*^ Vij + 7r/r is the potential of the excited state, which is a quantity 
also measured on the lattice along with the four-quark energies. 

In the special case of regular tetrahedra Vdd = Vaa = Vbb = Vcc and 
V reduces to the form 



Vaa 


-fVAA/2 


-fVAA/2 





-rva/2 


-rva/2 


-fVAA/2 


Vaa 


-fVAA/2 


-rva/2 





-rva/2 


-.fVAA/2 


-fVAA/2 


Vaa 


~rva/2 


~rva/2 








-rK/2 


~.rva/2 


Vb 


~rvj2 


~rv,/2 


~rva/2 





-rva/2 


~rvj2 


Vb 


-rvj2 


~rva/2 


-rK/2 





-rvc/2 


-rvc/2 


Vb 



(A.23) 

where Va , Vb , Vc can be expressed in terms of Vaa and v*{ij) plus some 
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fine tuning parameters. As with all geometries 



N 



1 

-//2 
-//2 


-//2 
1 

-//2 


-//2 
-//2 
1 




-P/2 
-/V2 


-/V2 


-r/2 


-r/2 
-r/2 






-/V2 
-/V2 




-/V2 


-/V2 

-/V2 



1 

-/72 
-/72 


-r/2 
1 

-r/2 


-r/2 
-r/2 
1 



(A.24) 



The full 6x6 matrix [V — + Vaa)N] now factorizes into three 2x2 
matrices, two of which are identical - giving the observed degeneracy. These 
two matrices have determinants of the form 



-E{l + f/2) 
-r{E~Va)/2 



-r{E~Va)/2 

-E{i + r/2) + Vb + rvc/2 



0, 



whereas the third 2x2 matrix has the determinant 



-E{1 

riE- 



-/) 

Va) 



r{E~Va) 

-E{i - n + Vb - 



= 0. 



(A.25) 



(A.26) 



In this case the problem reduces to solving two quadratic equations for E. 
However, away from the regular tetrahedron the complete 6x6 matrix needs 
to be treated directly. 

By fitting simultaneously the energies Eq and Ei from the lattice re- 
sults for the geometries in Fig. 11.101 an interquark potential model can be 
constructed that is able to explain, on the average, these energies. The 
full model utilized 6 basis states A, B, C, A* ,B* , C* and in its most general 
form has eight parameters. However, in practice, only 3 of these (fc/, fca, r) 
need be considered as completely free when fitting the data. 

The parameters that are of most interest are those connected with the 
ranges of the various interactions, namely, kf and ka- Here "range" is 
defined as r/_a,c = \/^/bskf^a,c- In the 2x2 version, where ka is effectively 
infinite, we get fc/(Ia)=0.57(l) i.e. rf{la) — 5.0 in lattice units of 0.12 fm. 
However, when the excited states A* ,B* , C* are introduced, the interaction 
between the basic states A,B,C decreases by raising kf to 1.51 giving 
r/ = 3.1. But at the same time this loss of binding by the direct interaction 
between A, B, C is compensated by their coupling to the A*,B*, C* states. 
This coupling in Eg. IA.22l is found to have about the same range = 5.1 
as r/(Ia) in Eg. IA.5| whereas the direct interaction between the A*,B*,C* 
■states seems, in all fits, to be satisfied with simply a two-quark description 
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without any four- quark correction (i.e. k^^O) in Eg. \A.2l\ The observation 
that r / (la) w ra suggests that the energy density has a range dictated by the 
longest range available — namely r^. Therefore, when the A* ,B* , C* states 
are not explicitly present, as in Model la, the only available range ?'/(Ia) 
has to simulate the role of r^. In the binding energies the contributions 
from the A*^B*,C* states rapidly dominate over those from the A,B,C 
states. For example, with squares of side R, the A, B, C states contribute 
only 85, 40, 10% to the binding energy for R— 2,4,6 respectively. Of course, 
at the largest distances (~ 0.7 fm) the quenched approximation is expected 
to break down and the role of quark-pair creation to become important. 



Appendix B: Extension of the /-Model to [{Qq){Qq)] Systems 

When only two of the four quarks are static the corresponding matrices 
for Q(ri)Q(r2)^(r3)q(r4) can be expressed in a similar form but where 
the matrix elements are now integrals over the positions of the two light 
antiquarks. In the notation of Fig. II. 41 we consider basis state A to be the 
one realised as two separate heavy-light mesons — [Q1Q3] and [Q294] — 
when the distance R = ri — r2 between the two heavy quarks becomes 
large. In this state the convenient coordinates are then Si = ra — ri and 
S2 = Y4, — Y2, whereas for the other partition B the convenient coordinates 
are ti = ra — r2 = Si + R and t2 r4 — ri = S2 — R. 

To describe this system in terms of an Effective Potential Theory the 
three ingredients quoted in Subsec. 11.1.2.21 are needed — see Ref. for 
more details of the specific example now to be described: 

A wave equation. 

For simplicity, the system is considered to be non-relativistic resulting in a 
Schrodinger-like equation 



i.e. a Resonating Group equation as discussed by Oka and Yazaki in Chap- 
ter 6 of Ref. and also Ref. |27j . This is a generalisation of Eq. 11.371 
to non-static quarks and can be solved using a variational wave function 



|K'(i?) -f- V'{R) - E{A, i?)N'(i?)|V = 



(B.l) 
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taken to have the form |140| 



^{n, /) = /i/2(ri, r2, rs, v^) ^ exp(-XM,X), 



(B.2) 



i=l 



where X = (si, S2, R) and each matrix Mi has the form 



^ I ftj Ci 
M, = - b, di e. 



(B.3) 



Since the present problem considers the masses of the hght quarks to be 
equal, it is sufficient to use a simplified form of Mi with bi — 0, di = and 
Ci — Ci. This is not necessary, but it is expected to be the dominant term 
in such a symmetric case. Already for = 2, this wave function is indeed 
adequate for giving sufficiently accurate four-quark binding energies. Even 
this choice involves five free parameters (ai, ci, a2, C2, 172) in the variation 
- with gi being fixed at unity to set the overall normalisation. In what 
follows the positions of the light quarks are integrated over leaving matrix 
elements that are functions of R. In order to achieve this in any practical 
way it is necessary to have a form for /(ri, r2, ra, that has a simple 
spatial dependence. Here the very symmetrical form in Eq. 11.441 defined 
by a single parameter fc/, is used i.e. 



(N.B. Here kj = fc/6, where k was defined in Eq. 11.441 ') It should be 
emphasised that this form of / is purely for numerical simplicity leading 
to analytical expressions for all matrix elements. As in the static case kf 
is a free parameter, which should be adjusted to fit the four-quark lattice 
energies. 

The wave function in Ea. IB. 21 is used for both states A and B. This is 
an approximation that appears to work well for the Q^(f system, since A 
and B are similar in structure for the R values of interest here. 

An interquark potential. 

This enters in three different contexts: 

1) As w(13), ^(24), w(14), v(23) in the Vqq potential. This is taken to be 




(B.4) 
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of the standard form in Eq. II. 61 namely 

ay(2,r) ^- "'^"^^^^^ +0.1649(36)r/a + 0.629(25), (B.5) 

r / a 

which gives a string energy of (445 MeV)^ for a = 0.18 fm. Tliis was 
obtained by fitting Vqq generated on a 16'^ x 24 lattice. Here the emphasis 
was to get a good fit over the important range of r ~ (2 — 4)a and is in 
contrast to the potential in Ref. which was designed to extract the 
string tension at large values of r. 

2) As w(34) for the Vqq potential. Here it is assumed to also be of the form 
in Eq. lR5l 

3) As w(12) for the Vqq potential. This was calculated from the same gauge 
configurations as the four-quark energies. In this case there was no need to 
fit Vqq with a function of R, since it is only ever needed at discrete values 
of i? - the ones for which the four-quark energies are calculated. 

An effective quark mass niq. 

In this case niq can be determined beforehand by carrying out an EPT 
analysis of the correponding two-body energies in Fig. 11.151 Using the po- 
tential in Eq. IB. 51 a value of niq « 400 MeV is able to give a good fit to the 
spin-averaged energies with L = 0, 1, 2 and 3. However, the results are 
not strongly dependent on niq — see Fig. 11.181 A disturbing feature of this 
result is that ruq <C 1 GeV indicating the need for a relativistic approach 
— as discussed in Sec. 11.121 

In Eg. IB. II the normalisation matrix N'(i?, fc/) — a generalisation of 
N(/) in Eg. 11.381 to non-static quarks in SU(3) — can now be written as 

where, after integrating over Si and S2, N{R, kf) can be expressed as a 
sum of terms of the form 



7r3 



(aX)3/- 



■ exp 



y2 



(B.7) 



where a = 0.5(ai -f- aj) + 3kf, c = 0.5(cj ± Cj) + 2fcy, d = 0.5(cj ± Cj) — 2kf, 
g = 0.5{gi + gj) + Akf, X ~ a — k'j / a, Y — c + kfd/a and Z ~ g — d^/a. 
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Since two of the quarks are not static there is now also a kinetic energy 
matrix in Eq. IB. II namely, 

K'(R k\^f K^iR^O) + K,{R,0) \[K^{R,kf)^Ki{R,kf)] \ 

>^{^,^f) { ^[Ks{R,kf) + K,{R,kf)] Ks{R,0) + K,{R,0) J' 

(B.8) 

where, for example. 



KsiR, kf) = J d?sidPs2ip*{kf) 



2mqdr^ 



i^ikf). (B.9) 



Again these integrals can be expressed in forms similar to that in Eq. IB. 71 
Finally, the potential matrix — a generalisation of V(/) in Eg. 11.381 to 
non-static quarks — has the form 

^^'''''^-^-[{V,s,kj) M14),0) + M23),0) j- (^-^0) 

Here 



{VAB,kf) = {VQ,)-{V,,)-{VQQ), 



where 



(Vq,) =i [(«(13), kf) + {v{2A),kf) + {v{U),kf) + {v{23), kf}] , 
{V,,)=l{v{3A),kf) and {Vqq) ^ ^N{R,kf)V{2, R). (B.ll) 



Here N{R,kf ) is defined in Eg. IB. 61 V{2,R) is the potential between the 
two heavy quarks and, for example. 



{v{13),kf) = J d'sid's2r{kf)V{s^)i;{kf). (B.12) 

For potentials of the form in Eq. 11.61 these integrals can be expressed in 
terms of Error functions. The energy E{A, R, kf) of the two heavy-light 
meson system is then obtained by diagonalising Eg. IB. II Since this equation 
is a 2 X 2 determinant, a prediction can also be made for an excited state 
£'*(4, R,kf) and the corresponding binding energy B*{A, R). 
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